以下のコードを考えてみましょう (MatLab):
w = 0 : 0.0001 : 9.4978;
a = [1 11 46 95 109 74 24];
b = [-1 3 4 3 1];
mu = 1;
a0 = a(7) ;a1 = a(6) ;a2 = a(5); a3 = a(4) ; a4 = a(3) ; a5 = a(2); a6 = a(1);
b0 = b(5);b1 = b(4);b2 = b(3) ; b3 = b(2); b4 = b(1) ;
De = -a6*w.^6 + a4*w.^4 - a2*w.^2 + a0;
Do = a5*w.^4 - a3*w.^2 + a1;
Ne = b4*w.^4 - b2*w.^2 + b0;
No = -b3*w.^2 + b1;
T = 0.01;
e = real((1i*w).^mu);
f = imag((1i*w).^mu);
A = Ne.*cos(T*w) + w.*No.*sin(T*w);
B = e.*(Ne.*cos(T*w) + w.*No.*sin(T*w)) - f.*(w.*No.*cos(T*w) - Ne.*sin(T*w));
C = w.*No.*cos(T*w) - Ne.*sin(T*w);
D = e.*(w.*No.*cos(T*w) - Ne.*sin(T*w)) + f.*(Ne.*cos(T*w) + w.*No.*sin(T*w));
Kp = (-De.*D + w.*Do.*B)./(f.*(Ne.^2 + w.^2.*No.^2));
Kd = (-w.*Do.*A + De.*C)./(f.*(Ne.^2 + w.^2.*No.^2));
figure
plot(Kp,Kd)
line([-24 -24],[-2.24 9.813])
コードを実行すると、次の図が得られます。
曲線の指定された部分(赤い部分、wは[0.6342、0.9985]に属します)
に接線を引きたい:
それを行った後、私の目的は、接線によって生成されるすべての可能な領域間のこの線と曲線によって定義される内向きの半平面の
最大領域を見つけることです(このように):
別の点に別の接線がある別の例は次のとおりです。
そして、最初の領域が 2 番目の領域よりも大きいと結論付けることができます。このアプローチは、赤い部分のすべてのポイントに対して行う必要があります。
MatLabでどうすればできますか?
私の質問が明確であることを願っています。任意のアイデアをいただければ幸いです。