版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

Math露珠

木虫 (知名作家)

应助: 10 (幼儿园)
金币: 7730.3
帖子: 8880
在线: 886小时
虫号: 4010226

[交流] Error using plot Vectors must be the same length.如下程序如何更改？

Error using plot Vectors must be the same length.如下程序如何更改？

hold on %figureaxes
>> q1=0.94;q2=0.94;q3=0.94; %x,y,z
>> h=0.05;T=250;  %hT
>> a=0.45;b=0.75; %
>> p1=16.5;p2=17.5;  %pp1p2
>> m1=1;m2=1;  %
>> x0=0.1;y0=0.001;z0=0.001;  %
>> x0=1;y0=1;z0=2; %
>> p=p1; %p%----  -----
>> x(1)=x0+h^q1*y0;
>> y(1)=y0+h^q2*(x0-b*y0-x0*z0);
>> z(1)=z0+h^q3*(-a*z0+x0*x0);
>> t=m1*h;n=1;  %
>> % ------  --------
>> while t<=T  %
n=n+1; %
M1=0;M2=0;M3=0; %
  % xyz
for k=1:n-1
      M1=M1+(-1)^(k)*gamma(q1+1)/(gamma(k+1)*gamma(q1-k+1))*(x(n-k)-x0);
      M2=M2+(-1)^(k)*gamma(q2+1)/(gamma(k+1)*gamma(q2-k+1))*(y(n-k)-y0);
      M3=M3+(-1)^(k)*gamma(q3+1)/(gamma(k+1)*gamma(q3-k+1))*(z(n-k)-z0);
end
%5m1*hm2*h
if p==p1
   p=p2;t=t+m2*h;
else
      p=p1;t=t+m1*h;
end
%
x(n)=x0-M1+h^q1*y(n-1);
y(n)=y0-M2+h^q2*(x(n)-b*y(n-1)-x(n)*z(n-1));  %pp1p2
z(n)=z0-M3+h^q3*(-a*z(n-1)+x(n)*x(n));
end  %
>> t=0:h:T; %
>> figure(13)  %x(t)
>> plot(t,x,'K');xlabel('t'),ylabel('x');

错误提示：
Error using plot
Vectors must be the same length.

提供帮助的虫子50金币奖励，谢谢大家

回复此楼

» 猜你喜欢

售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有4人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有7人回复
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急已经有6人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有8人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有7人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有8人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有5人回复
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有8人回复
一个有机合成实验室都需要哪些设备？已经有9人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有7人回复

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2017-10-24 03:03:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

简单回复

M塞北2楼

2017-10-24 05:09 回复

Math露珠(金币+1): 谢谢参与

发自小木虫Android客户端

dsctg3楼

2017-10-24 08:40 回复