matlab 求解微分方程

显示全部楼层 · 2011-6-4 17:20:04

数值解ode45求解微分方程x“+0.06x'=0.35sin(wt+0.07)
初始条件x‘(0)=0,x(0)=0.03,t范围是[0,500],要求输出图像

千问 · 2011-6-4 17:20:04

编个M文件，命名成rigid：%x(1)指原微分方程中的x，x(2)指x'function dx = rigid(t,x)dx=zeros(2,1);%两行dx(1)=x(2);dx(2)=-0.06*x(1)+0.35*sin(2*t+0.07); %令w=2,必须是个给定常数在命令窗口里输入命令：[T,Y]=ode45(@rigid,[0,500],[0.03,0]);plot(T,Y(:,1))就输出x的图像...