我想把这段代码,但我的值显示为零,图形显示。如何修复它

22日视图(30天)
clc;明确;
%输入参数
δ= 50;
K1 = 10 ^ 4;
Ko = 0.1;
n = 3;
Oa = 10;
Pa = 100;
mu_1 = 10 ^ 3;
K2 = 5 * 10 ^ 4;
K3 = 10 ^ 4;
γ= 75;
K = 100;
%输入初始条件
M1 (1) = 10;
平方米(1)= 0;
M3 (1) = 0;
%输入时间
t (1) = 0;
h = 0.01;%的时间间隔
dt = 0: h: 100;%时间跨度
%输入空数组
t = 0(长度(dt), 1);% t空数组
M1 = 0(长度(dt), 1);%为M1空数组
M2 = 0(长度(dt), 1);%对M2空数组
M3 = 0(长度(dt), 1);%对M3空数组
K = 0(长度(n - 1), 1);
j = 1:长度(dt)
t (j + 1) = t (j) + h * j;
M11公路= M1 (j) + 1. / (1 + exp (dt));
结束
i = 1:长度(n - 1)
K (i + 1) = K (i) + n *我;
结束
j = 1:长度(dt (j))
M1 (j + 1) = M1 (j) + t (j) *(δ* M1 (j) * (1 - (M1 (j) /γ)]2 * K1 * M1 (j) * M1 (j) M1 (j) * (K (i + 1)。* M1 (j)) - (Oa-n) * K3 * M1 (j) * M3 (j)——(Pa-Oa) * Ko * M1 (j) * 10 - (mu_1 * M1 (j)))
结束
M1 = 10001×1
0 0 0 0 0 0 0 0 0 0
情节(dt, M1 (j + 1),“r”,“线宽”3)
包含(“时间”)
ylabel (“M1”)
4评论

登录置评。

接受的答案

KALYAN ACHARJYA
KALYAN ACHARJYA 2023年5月13日13:26
编辑:KALYAN ACHARJYA 2023年5月13日13:32
M1必须向量,然而在代码中情节M1 (j + 1)是一个标量(单值)
情节(dt, M1,“r”,“线宽”3);
理想的情节,请修改代码。
2的评论
KALYAN ACHARJYA
KALYAN ACHARJYA 2023年5月13日14:34
是的,你必须检查代码。请再交叉校验是否这样做是正确的。

登录置评。

更多的答案(0)

类别

找到更多的在MATLAB移动基础帮助中心文件交换

标签

社区寻宝

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

开始狩猎!