如果elseif else语句不工作

19日视图(30天)
Yakoob说道
Yakoob说道 2023年4月25日16:16
评论道: 明星黾 2023年4月25日17:47
我有这段代码
r0 = 0.05;
k1 = 0.5;
k2 = 0.5;
K = 1;
μ= 0.5;
ρ= 0.5;
α= 0.1;
b = 0.05;
ε= 0.25;
σ= 0.05;
η= 0.05;
π= 0.05;
ω= 0.001;
信谊M M Msol正
P (M) =ε* M / (1 + M);
N (M) =μ+(ρ* M / (1 + M));
λ(M) = 1 / (1 + N);
r (M) = r0 * (1 + k1 * N * (1-k2 * N));
F = (m / b) * r (m) *(η+ P (m)) * (1 + N (m)) *(1 -(((ασ* * m *(η+ P (m)) * (1 + N (m)) + b * (m +α)*(π* M-omega0))) / (b *ασ* * K)))——(m /σ)*(π* M-omega0) = = 0;
(num,窝)= numden (lh (F));
m_num = linspace (0.01, 0.9, 50);
u = 1:元素个数(m_num)
Msol (u) = vpa(解决(潜艇(num, m, m_num (u))));
结束
R0 = (b。* k *λ(omega0. / pi))。/ (m_num。*(埃塔+ P (omega0. / pi)));
如果R0 < 1
l = 0;
y = 0;
v = 0;
其他的
l = ((b。* k *λ(omega0. / pi)。*(π。* Msol-omega0))。/(ασ*。* R0。*(埃塔+ P (omega0. / pi))));
y =((π。* Msol - ((pi +₀) /ε)。* (((b。* k *λ(omega0. / pi))。/ (m_num * R0)) eta))。/(σ));
v = ((m_num * R0。*(埃塔+ P (omega0. / pi)) *(π。* Msol-omega0))。/(σ。* k *λ(omega0. / pi)。*(埃塔+ P (Msol))));
结束
次要情节(1,1)
p1 =情节(R0、l);
ylabel (“l ^ *’)
xlim ([0 5]);
ylim ([-0.15 - 0.15]);
网格;
次要情节(1、2)
p2 =情节(R0, y);
ylabel (“y ^ *’)
xlim ([0 5]);
ylim ([-0.15 - 0.15]);
网格;
次要情节(3、1,3)
p3 =情节(R0, v);
包含(“R_0”),ylabel (“v ^ *”)
xlim ([0 5]);
ylim ([-0.15 - 0.15]);
网格;
集(p1, {“线宽”},{1.5});
集(p2, {“线宽”},{1.5});
集(p3, {“线宽”},{1.5});
elseif,我试图使用如果其他语法l, y v 0 R0 < 1和积极的其他地方。然而,它看起来像MATLAB并不承认这部分。任何帮助将不胜感激!

接受的答案

明星黾
明星黾 2023年4月25日16:24
编辑:明星黾 2023年4月25日16:54
使用一个 逻辑 向量来选择您要绘制的元素,然后用这个作为一个参数在每个下标 情节 电话-
r0 = 0.05;
k1 = 0.5;
k2 = 0.5;
K = 1;
μ= 0.5;
ρ= 0.5;
α= 0.1;
b = 0.05;
ε= 0.25;
σ= 0.05;
η= 0.05;
π= 0.05;
ω= 0.001;
信谊M M Msol正
P (M) =ε* M / (1 + M);
N (M) =μ+(ρ* M / (1 + M));
λ(M) = 1 / (1 + N);
r (M) = r0 * (1 + k1 * N * (1-k2 * N));
F = (m / b) * r (m) *(η+ P (m)) * (1 + N (m)) *(1 -(((ασ* * m *(η+ P (m)) * (1 + N (m)) + b * (m +α)*(π* M-omega0))) / (b *ασ* * K)))——(m /σ)*(π* M-omega0) = = 0;
(num,窝)= numden (lh (F));
m_num = linspace (0.01, 0.9, 50);
u = 1:元素个数(m_num)
Msol (u) = vpa(解决(潜艇(num, m, m_num (u))));
结束
R0 = (b。* k *λ(omega0. / pi))。/ (m_num。*(埃塔+ P (omega0. / pi)));
l = ((b。* k *λ(omega0. / pi)。*(π。* Msol-omega0))。/(ασ*。* R0。*(埃塔+ P (omega0. / pi))));
y =((π。* Msol - ((pi +₀) /ε)。* (((b。* k *λ(omega0. / pi))。/ (m_num * R0)) eta))。/(σ));
v = ((m_num * R0。*(埃塔+ P (omega0. / pi)) *(π。* Msol-omega0))。/(σ。* k *λ(omega0. / pi)。*(埃塔+ P (Msol))));
Lv = R0 > = 1;%逻辑向量
次要情节(1,1)
p1 =情节(R0 (Lv), l (Lv));
持有
情节(1.02 [0],[0 0),“- - -”,“颜色”p1.Color,“线宽”,1)
持有
ylabel (“l ^ *’)
xlim ([0 5]);
ylim ([-0.15 - 0.15]);
网格;
次要情节(1、2)
p2 =情节(R0 (Lv), y (Lv));
持有
情节(1.02 [0],[0 0),“- - -”,“颜色”p2.Color,“线宽”,1)
持有
ylabel (“y ^ *’)
xlim ([0 5]);
ylim ([-0.15 - 0.15]);
网格;
次要情节(3、1,3)
p3 =情节(R0 (Lv)、v (Lv));
持有
情节(1.02 [0],[0 0),“- - -”,“颜色”p3.Color,“线宽”,1)
持有
包含(“R_0”),ylabel (“v ^ *”)
xlim ([0 5]);
ylim ([-0.15 - 0.15]);
网格;
集(p1, {“线宽”},{1.5});
集(p2, {“线宽”},{1.5});
集(p3, {“线宽”},{1.5});
编辑- (2023年4月25日16:54)
添加常数台词 1.02 [0] 在每一个情节。
4评论
明星黾
明星黾 2023年4月25日17:47
像往常一样,这是我的荣幸!

登录置评。

更多的答案(0)

标签

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!

翻译的