我写了以下(写在我的问题末尾)一段没有错误的代码,但我认为,在运行它时,它有一个覆盖问题。在程序中,有两种情况我想画图;首先,用 ezplot 编写的曲线图,其次,我想绘制回归线的绘图回归。
当我跳过代码plotregression(C_i, D_i)
时,显示所有五个逻辑函数的图表没有问题(实际上这里的一位用户向我展示了暂停代码以帮助执行此操作),但是当我合并时plotregression(C_i, D_i)
,会发生两件事:
它向我展示了所有的回归线,但与所有回归线都在同一个图中相反,它不断改变回归线的回归系数不同。如果您运行代码,您实际上可以看到这种情况。
的效果
plotregression(C_i, D_i)
消失了;它不再绘制五个逻辑函数的图形。
我有两个问题:
如果我想得到两个数字,一个显示所有 5 条逻辑曲线,另一个显示所有 5 条回归曲线,我怎样才能最小限度地修改程序以完成工作?
我怎样才能停止写回归曲线?我使用了'hold on-hold off'来避免逻辑曲线相同,但它不适用于回归曲线吗?
这是代码:
syms t;
hold on;
for i=1:5;
P_i=0.009;
r_i=abs(sin(i.^-1));
y_i(t)= P_i*exp(r_i*t)/(1+P_i*(exp(r_i*t)-1));
t_1= 1+rand; t_2= 16+rand; t_3=31+rand;
time_points=[1, t_1; 1, t_2; 1, t_3];
biomarker_values= double([y_i(t_1);y_i(t_2);y_i(t_3)]);
X=vertcat(X,time_points);
Z=blkdiag(Z,time_points);
Y=vertcat(Y,biomarker_values);
G=vertcat(G,[i,i,i]');
ezplot(y_i,[-50,100]);
C_i=time_points(:,2)
D_i=biomarker_values
plotregression(C_i, D_i)
结尾;
hold off;