作业帮 > 综合 > 作业

请高手指教Matlab

来源:学生作业帮 编辑:搜搜考试网作业帮 分类:综合作业 时间:2024/06/28 07:29:39
请高手指教Matlab
function dx=c7chua(x,a,b,alpha,beta,gamma)
dx=[alpha*(x(2)-x(1)-nlfun(x,a,b));
x(1)-x(2)+x(3);-beta*x(2)-gamma*x(3)];
function f0=nlfun(x,a,b)
f0=b*x+(a-b)*(abs(x+1)-abs(x-1))/2;
>>a=-120/7; b=-75/7; alpha=15; beta=20; gamma=0.5;
f=odeset; f.RelTol=1e-6;
[t,x]=ode45('c7chua',[0,100],[-2.121304;-0.06617;2.88109],f,a,b,alpha,beta,gamma);
plot3(x(:,1),x(:,2),x(:,3))
grid
Input argument "alpha" is undefined.
Error in ==> c7chua at 2 dx=[alpha*(x(2)-x(1)-nlfun(x,a,b));
请高手指教Matlab
这应该是蔡氏电路
定义函数错误,nlfun参数x是个三维向量,而f0仅需要x(1),可改为如下
function dx=c7chua(t,x,a,b,alpha,beta,gamma)    f0=b*x(1)+(a-b)*(abs(x(1)+1)-abs(x(1)-1))/2;    dx=[alpha*(x(2)-x(1)-f0);    x(1)-x(2)+x(3);-beta*x(2)-gamma*x(3)];
 a=-120/7; b=-75/7; alpha=15; beta=20; gamma=0.5;f=odeset; f.RelTol=1e-6;[t,x]=ode45(@c7chua,[0,100],[-2.121304;-0.06617;2.88109],f,a,b,alpha,beta,gamma);plot3(x(:,1),x(:,2),x(:,3))grid
参数好像有问题图应是这样的而你的却得出