用matlab 解这个简单微分方程的数值解，急
y'=-y+t+1,y(0)=1第一步要求先解解析解，再解数值解，还要画出解析的图形曲线和数值解图形，求所有的详细代码和图形。说详细点，新手不懂。...
y'=-y+t+1,y(0)=1
第一步要求先解解析解，再解数值解，还要画出解析的图形曲线和数值解图形，求所有的详细代码和图形。说详细点，新手不懂。 展开
第一步要求先解解析解，再解数值解，还要画出解析的图形曲线和数值解图形，求所有的详细代码和图形。说详细点，新手不懂。 展开
展开全部
clc
clear all
close all
%主程序：
yy=dsolve('Dy=-y+t+1','y(0)=1','t')
tt=linspace(0,1,1001); % t是[0,0.001,0.002,...,1]
for i=1:1001
t=tt(i);
f1(i)=eval(yy);
end
figure(1)
plot(tt,f1,'r')
title('Ana')
%解析解画完
%数值解
[ttt,f2] = ode45(@rigid,tt,1);
figure(2)
plot(tt,f2,'b')
title('Num')
===================================
%存成rigid.m
function dy=rigid(t,y)
dy=-y+t+1;
===================================
一共两个文件main.m, rigid.m, 不明白可追问
clear all
close all
%主程序：
yy=dsolve('Dy=-y+t+1','y(0)=1','t')
tt=linspace(0,1,1001); % t是[0,0.001,0.002,...,1]
for i=1:1001
t=tt(i);
f1(i)=eval(yy);
end
figure(1)
plot(tt,f1,'r')
title('Ana')
%解析解画完
%数值解
[ttt,f2] = ode45(@rigid,tt,1);
figure(2)
plot(tt,f2,'b')
title('Num')
===================================
%存成rigid.m
function dy=rigid(t,y)
dy=-y+t+1;
===================================
一共两个文件main.m, rigid.m, 不明白可追问
本回答由提问者推荐
已赞过
已踩过<
评论
收起
你对这个回答的评价是？
富港检测技术（东莞）有限公司_
2024-04-02 广告
正弦振动多用于找出产品设计或包装设计的脆弱点。看在哪一个具体频率点响应最大（共振点）；正弦振动在任一瞬间只包含一种频率的振动，而随机振动在任一瞬间包含频谱范围内的各种频率的振动。由于随机振动包含频谱内所有的频率，所以样品上的共振点会同时激发...
点击进入详情页
本回答由富港检测技术（东莞）有限公司_提供
推荐律师服务：
若未解决您的问题，请您详细描述您的问题，通过百度律临进行免费专业咨询
