Chua电路的状态方程如下：
参数：a=10,b=15,m0=-1/7,m1=2/7.
程序和simulink仿真图我放在下边，仿真出的图形与之前用function编写的图形大致一样但是曲线曲折，不平滑，希望有前辈能解答一下，谢谢！

这是主函数：
[t,x] = sim('simulink_chua',500);
n = length(t);
n1 = round(n/2);
figure(1);
plot(t(n1:n,1),x(n1:n,1));
xlabel('t','fontsize',20,'fontname','times new roman','FontAngle','italic');
ylabel('x','fontsize',20,'fontname','times new roman','FontAngle','italic');
figure(2);
plot(x(n1:n,1),x(n1:n,2));
xlabel('x','fontsize',20,'fontname','times new roman','FontAngle','italic');
ylabel('z','fontsize',20,'fontname','times new roman','FontAngle','italic');

这是将状态方程转化成的simulink：

这是f(x)子模块：

贴吧不能添加附件，若需要源代码程序，私聊我给您发一下

奥对了，这是之前同function编出来的图案：

偶可以试试