Question

怎么给用MATLAB画微分方程的图像？

函数是D2y+3siny=0;y是关于时间的函数，我想画t=[0-5]的图像。在线等程序！！！

Answer 1

把下面的程序放到一个m文件里，直接运行

function [T,Y] = stest

tspan = [0,5];

Dy0 = 0;

y0 = 30*pi/180;

[T,Y] = ode45(@f,tspan,[y0,Dy0]);

plot(T,Y);

xlabel('t/s');

ylabel('y,dy')

legend('y','dy')

function dy = f(t,y)

dy = zeros(2,1);

dy(1) = y(2);

dy(2) = -3*sin(y(1));

Answer 2

给出的微分方程的初值条件为多少？

追问

t=0,y=30°；t=0时，Dy=0;谢谢。

追答

y0=[pi/6,0]

h=0.1;

a=0;

b=5;

[t,y] = runge_kutta(@myfun,y0,h,a,b)

plot(t,y),grid on

xlabel('t/s');

ylabel('y,dy')

legend('y','dy')

 myfun.m

function f = myfun(t,y)

f=[y(2);-3*sin(y(1))];

追问

你好，我怎么运行不出来这个程序呢?

追答

把runge_kutta.m和myfun.m放在同目录下，在此目录下执行
y0=[pi/6,0]
h=0.1; a=0; b=5;
[t,y] = runge_kutta(@myfun,y0,h,a,b)
。。。。。。