| 查看: 733 | 回复: 6 | |||
| 当前主题已经存档。 | |||
semibeauty木虫 (正式写手)
|
[交流]
【求助】偏微分方程求解【已完成】
|
||
|
下面两个偏微分方程，要用matlab求解，解用图形表示，但总运行不出来，故求助（小虫对这些就是文盲，是我们老师让我去咨询高人，有什么低级错误，大家多担待） 方程1方程2分别对应附件1和2. 方程1 m-file: function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); [T,Y] = ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') 方程2 m-file: function dy = vdp1000(t,y) dy = zeros(2,1); % a column vector dy(1) = y(2); dy(2) = 1000*(1 - y(1)^2)*y(2) - y(1); options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]); [T,Y] = ode15s('vdp1000',[0 3000],[2 0],options); plot(T,Y(:,1),'-o') 上面是我把我们老师发给我的贴上。她说是帮助文件的两个例子，我们要解的方程和这个类似，只是未知数多了一点，现在两个示例都运行不出来，以前曾经运行出来过。 [ Last edited by nono2009 on 2009-9-25 at 07:16 ] |
回复此楼» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有7人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
sitonyl
铁杆木虫 (文坛精英)
- 应助: 33 (小学生)
- 金币: 12935.9
- 红花: 3
- 帖子: 10448
- 在线: 338.8小时
- 虫号: 722261
- 注册: 2009-03-14
- 性别: GG
- 专业: 导航、制导与传感技术
2楼2009-05-20 09:59:59
3楼2009-05-20 14:36:36
4楼2009-05-20 22:27:01
semibeauty
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3602.8
- 散金: 35
- 红花: 1
- 帖子: 653
- 在线: 53.7小时
- 虫号: 651761
- 注册: 2008-11-11
- 专业: 质谱分析
5楼2009-05-21 09:25:18
★ ★ ★ ★ ★ ★ ★ ★ ★
woshilsh(金币+4,VIP+0):thanks, 5-28 15:22
semibeauty(金币+5,VIP+0):谢谢你的详细解答，谢谢 5-31 20:57
woshilsh(金币+4,VIP+0):thanks, 5-28 15:22
semibeauty(金币+5,VIP+0):谢谢你的详细解答，谢谢 5-31 20:57
|
这个问题很好办。有两种方法解决：(针对第一个例子） 1.首先建立函数文件。保存成rigid.m，内容如下： function dy = rigid(t,y) dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); 然后在命令行里面调用 [T,Y] = ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') 2.第二种方法：编写函数文件，保存为fangcheng1.m，内容如下： function fangcheng1 [T,Y] = ode45('rigid',[0 12],[0 1 1]); plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.') function dy = rigid(t,y) （作为一个子函数，被ode45调用） dy = zeros(3,1); % a column vector dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3) = -0.51 * y(1) * y(2); 然后在命令行里输入fangcheng1运行即可。 |
6楼2009-05-28 14:58:41
7楼2009-05-28 15:04:28