24小时热门版块排行榜    

查看: 339  |  回复: 3
当前主题已经存档。

freecn

[交流] 【求助】关于matlab微分方程解的存储和调用问题

假设我已经用伦格-库塔法算出一组微分方程组的解(3个微分方程)
tspan=[0,20] %起始和终了的时间
ek4(@dy,[0,0,1],1000) %ek4是伦格-库塔子程序,[0,0,1]是初值
plot(t,y(:,1),'-')
并且能够画出其中任意一个解随时间的变化关系
问题:怎么存储任意一个解随时间的变化关系,存储成ASCII的?
好像,save t, y(:,1) /-ascii 不管用

第二个问题:
如果要算 sin(t)*y(t),怎么办?也就是说要把时间配套起来(sin(t)和y(t)的时间)
tstep=(tspan(2)-tspan(1))/1000
for n=1,1,1000
z=sin(tstep*n)
z1=z*y(:,1)
end

可是程序报错说维数不匹配
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

★ ★
wangen994(金币+2,VIP+0):感谢参与讨论 7-16 08:56
tstep=(tspan(2)-tspan(1))/1000
n=1:1:1001
z=sin(tstep*n)
z1=z*y(:,1)
2楼2009-07-15 22:34:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

change0618

铁杆木虫 (著名写手)

方丈大师

save t y(:,1) -ascii
3楼2009-07-15 22:36:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freecn

好像不行
4楼2009-07-16 07:15:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 freecn 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见