インデックスが配列要素数を超えています。インデックスは 0 を超えてはなりません。
43 visualizzazioni (ultimi 30 giorni)
Mostra commenti meno recenti
Masaaki Yanagihara
il 28 Mar 2024
Commentato: Dyuman Joshi
il 28 Mar 2024
syms s
for kq = -10:1:10
eqnqroop = 1 + (kq*(-((4483*s^3)/5000 + (1915281*s^2)/5000000 - (551896897*s)/5000000000)/(s^4 + (52*s^3)/125 ...
+ (764541*s^2)/1000000 + (4137177603*s)/500000000 - 427810011/250000000))) == 0;
q1 = solve(eqnqroop,s); %%求解
kq
plot(real(q1(1)),imag(q1(1)),'-s','MarkerSize',5,'MarkerEdgeColor','red', 'MarkerFaceColor',[1 .6 .6])
hold on
end
kqを-10~10まで1刻みで変えながら,eqnqroopという方程式の解を複素平面にプロットする,ということをしたいのですが,「インデックスが配列要素数を超えています。インデックスは 0 を超えてはなりません。」とでてきてしまい,-10~0までで計算がストップしてしまいます.
解決方法をご教授いただけると幸いです.
0 Commenti
Risposta accettata
Dyuman Joshi
il 28 Mar 2024
Modificato: Dyuman Joshi
il 28 Mar 2024
When kq is 0, the equation becomes 1==0, which is obviously not true, thus there is no solution for it. For that iteration, q1 is empty and you can't use indexing for it.
You can skip that value or plot a value manually for it. I have implemented the 1st option in my code below.
Note that the solve() returns the real solutions first if there are any. As the equation is a cubic polynomial in s, there will atleast be one real solution, so use the 2nd index for plotting.
syms s
figure
hold on
%Modify the loop indices to skip the value 0
for kq = setdiff(-10:10, 0)
eqnqroop = 1 + (kq*(-((4483*s^3)/5000 + (1915281*s^2)/5000000 - (551896897*s)/5000000000)/(s^4 + (52*s^3)/125 ...
+ (764541*s^2)/1000000 + (4137177603*s)/500000000 - 427810011/250000000))) == 0;
q1 = solve(eqnqroop,s);
%updated indexing
plot(real(q1(2)),imag(q1(2)),'-s','MarkerSize',5,'MarkerEdgeColor','red', 'MarkerFaceColor',[1 .6 .6])
end
2 Commenti
Dyuman Joshi
il 28 Mar 2024
You're welcome!
If my answer solved your problem, please consider accepting the answer.
Più risposte (0)
Vedere anche
Categorie
Scopri di più su 一般の演算 in Help Center e File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!