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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4325)
>基金申请 (3992)
>休闲灌水 (1361)
>虫友互识 (268)
>论文投稿 (244)
>导师招生 (230)
>有奖起名 (165)
>文献求助 (149)
>硕博家园 (90)
>健康生活 (87)
>考博 (81)
>留学DIY (37)
>考研 (31)
>公派出国 (28)
>博后之家 (25)
>绿色求助(高悬赏) (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金币奖励，谢谢大家

回复此楼

» 猜你喜欢

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

查看全部散金贴

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 回复