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

函数是D2y+3siny=0;y是关于时间的函数,我想画t=[0-5]的图像。在线等程序!!!... 函数是D2y+3siny=0;y是关于时间的函数,我想画t=[0-5]的图像。在线等程序!!! 展开
 我来答
slassddd
推荐于2017-12-15 · TA获得超过499个赞
知道小有建树答主
回答量:425
采纳率:0%
帮助的人:237万
展开全部

把下面的程序放到一个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));

本回答被提问者和网友采纳
已赞过 已踩过<
你对这个回答的评价是?
评论 收起
lhmhz
高粉答主

2015-06-27 · 专注matlab等在各领域中的应用。
lhmhz
采纳数:7263 获赞数:17004

向TA提问 私信TA
展开全部
给出的微分方程的初值条件为多少?
更多追问追答
追问
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))];

已赞过 已踩过<
你对这个回答的评价是?
评论 收起
收起 1条折叠回答
推荐律师服务: 若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询

为你推荐:

下载百度知道APP,抢鲜体验
使用百度知道APP,立即抢鲜体验。你的手机镜头里或许有别人想知道的答案。
扫描二维码下载
×

类别

我们会通过消息、邮箱等方式尽快将举报结果通知您。

说明

0/200

提交
取消

辅 助

模 式