用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/06 18:52:55
![用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来](/uploads/image/z/3765918-30-8.jpg?t=%E7%94%A8matlab%E6%B1%82%E8%BF%99%E4%B8%AA%E6%96%B9%E7%A8%8B%E7%9A%84%E6%95%B0%E5%80%BC%E8%A7%A3%EF%BC%9Ay%27%27-%281-y%5E2%29y%27%2By%3D0%2Cy%281%29%3D2%2Cy%282%29%3D0%2C%5B0%2C20%5D%2C%E9%BA%BB%E7%83%A6%E6%8A%8A%E7%A8%8B%E5%BA%8F%E7%BC%96%E5%87%BA%E6%9D%A5%E7%94%A8matlab%E6%B1%82%E8%BF%99%E4%B8%AA%E6%96%B9%E7%A8%8B%E7%9A%84%E6%95%B0%E5%80%BC%E8%A7%A3%EF%BC%9Ay%27%27-%281-y%5E2%29y%27%2By%3D0%2Cy%281%29%3D2%2Cy%282%29%3D0%2C%5B0%2C20%5D%2C%E9%BA%BB%E7%83%A6%E6%8A%8A%E7%A8%8B%E5%BA%8F%E7%BC%96%E5%87%BA%E6%9D%A5)
用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来
用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来
用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,
y(1)=2,y(2)=0,[0,20],麻烦把程序编出来
用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来用matlab求这个方程的数值解:y''-(1-y^2)y'+y=0,y(1)=2,y(2)=0,[0,20],麻烦把程序编出来
这是你们老师出的题目吧,这个例子在好多书上都出现过.
具体如下:
function dydt = vdp1(t,y)
dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];
end
保存函数文件
[t,y] = ode45(@vdp1,[0 20],[2; 0]);
%This example uses @ to pass vdp1 as a function handle to ode45.The
resulting output is a column vector of time points t and a solution array y.Each row in y corresponds to a time returned in the corresponding row of t.
作图:
You can simply use the plot command to view
the solver output:
plot(t,y(:,1),'-',t,y(:,2),'--')
xlabel('time t');
ylabel('solution y');
legend('y_1','y_2')