plotyの同じ数字をプロットする

Question

おはよう、 私はこの数字をMatlabの同じコードでプロットしようとしています。しかし、それは全く機能しません。私は単に同じ数字を得ることはできませんが、代わりにplotyyを使用します。 私がplotyyでそれを必要とするのは、同じ関数をプロットして定数を掛けて左手の数字の意味を読者に理解させるためです。右手側は「ゲイン」と呼ばれ、したがって、読者は少なくとも特定のパワーに対してどのくらいの利得が達成されるかを知るでしょう。ここ は、コードは次のとおりです。

[![r=1:0.001:3; 
R=50; 
kappa=0.5; 
L1=1000; 
c1=R/(8*kappa*L1); 
Frmin=(r+r.^-1).*(r-1)./(r+1); 
Frmax=(r+r.^-1).*(atan(r)-atan(1./r)); 
P1min=c1.*Frmin;% the 1st function for the left hand side of plotyy 
P1max=c1.*Frmax;% the 2nd function for the left hand side of plotyy 
YYmin=P1min*20; % the 1st function for the right hand side of plotyy 
YYmasx=P1max*20;% the 2nd function for the right hand side of plotyy 
figure; 
ha1 = area(\[1 3\], \[0 M.*c1\]); 
hold on 
plot(r,P1min,'r','Linewidth',8) 
fill(\[r fliplr(r)\],\[P1min,fliplr(P1max)\],'y') 
plot(r,P1max,'k','Linewidth',5); 
hold off 
legend('\Delta >0','\Delta=0', '\Delta<0','Pmax') 
%,'900','800','700','600','500') 
xlabel('r') 
ylabel('P_p (W)') 
title('P_p vs r for L= 1000 m')][1]][1] 

だから私は、私はこれでいくつかの助けを得ることができると思います。私は試行錯誤の試行をいくつか見せてもらうことができますが、同じ色を使って同じ図をプロットするのではなく、plotyを使って私の目標を明確にすることはできません。ありがとう、長いメッセージをお詫び申し上げます。

Answer 1

これはあなたが探しているものですか？あなたのコードの

編集（これが働いている）：黒とコントロールのフォントへ

r=1:0.001:3; 
R=50; 
kappa=0.5; 
L1=1000; 
c1=R/(8*kappa*L1); 
Frmin=(r+r.^-1).*(r-1)./(r+1); 
Frmax=(r+r.^-1).*(atan(r)-atan(1./r)); 
P1min=c1.*Frmin;% the 1st function for the left hand side of plotyy 
P1max=c1.*Frmax;% the 2nd function for the left hand side of plotyy 
YYmin=P1min*20; % the 1st function for the right hand side of plotyy 
YYmasx=P1max*20;% the 2nd function for the right hand side of plotyy 
figure; 
M = 1; 
ha1 = area([1 3], [0 M.*c1]); 
hold on 
plot(r,P1min,'r','Linewidth',8) 
fill([r fliplr(r)],[P1min,fliplr(P1max)],'y') 
plot(r,P1max,'k','Linewidth',5); 

あなたの要因とplotyy

myFactor1 = 5; 
[ax,~,~] = plotyy(r,P1min,r,myFactor1*P1min); 

設定が色を軸

set(ax,{'ycolor'},{'k';'k'},{'FontSize'},{8;15}) 

残りコード、編集、現在作業中

legend('\Delta >0','\Delta=0', '\Delta<0','Pmax') 
xlabel('r') 
ylabel('P_p (W)') 
title('P_p vs r for L= 1000 m')