版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

suixin123456

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 342.4
红花: 47
帖子: 96
在线: 23.2小时
虫号: 1578539
注册: 2012-01-15
性别: MM
专业: 粒子物理学和场论

[求助] 求助：用Matlab解微分方程，命令怎么编辑？有例子已有1人参与

如图，a=0.1,b=0.086,g=0.045,c=0.01,d=0.8,e=-0.069+0.109*i(复数）
当t=0时，y1=0.0005-0.00013*i,y2=-6.28-15*i
要求用matlab求此微分方程组的解，形式为exp()+exp()+复数
给我一个可运行的完整命令，谢谢了！！！

图片1.png

回复此楼

» 猜你喜欢

比利时鲁文大学固体物理系招收CSC博士已经有0人回复
🌟 比利时新鲁汶大学（UCLouvain）诚邀CSC博士加入Pascal Gehring教授团队已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有239人回复
第一性原理计算方向2026级博士申请 PRB*1，四级484 已经有1人回复
求助VISSIM破解版软件已经有0人回复
求2026年在台湾举行的物理和材料领域国际学术会议信息已经有0人回复
求国际会议网站已经有1人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

精诚所至，金石为开

1楼2014-02-18 16:44:54

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

回帖支持 ( 显示支持度最高的前 50 名 )

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

感谢参与，应助指数 +1

程序如下：

CODE:

function re_dsolve

a=0.1;

b=0.086;

g=0.045;

c=0.01;

d=0.8;

e=-0.069+0.109*i;

[y1 y2]=dsolve('Dy1 = (-b+a*i)*y1+g*y2', 'Dy2 =(-d/2+c*i)*y2-g*y1+e*i')

结果的一部分见附图1。

附图1.jpg

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

3楼2014-02-18 18:24:54

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

月只蓝

主管区长 (职业作家)

专家经验: +1059
计算强帖: 8
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
suixin123456: 金币+20, ★★★★★最佳答案 2014-02-23 16:18:18

引用回帖:

4楼: Originally posted by suixin123456 at 2014-02-20 11:31:37
结果的形式倒是对，不过abcd都赋值了为什么结果中还有符号？而且初始值都给了，还有c1c2？我就想要最后结果，不知道怎么弄，之前我也能算到带C1C2的，但是我需要用这个算几十组数据的最后结果...

已加上初值条件，并把符号解转换为数值解，代码如下，结果的一部分见附图1。

CODE:

function test3333

clear all;clc

syms a b c d e g

[y1 y2]=dsolve('Dy1 = (-b+a*i)*y1+g*y2', 'Dy2 =(-d/2+c*i)*y2-g*y1+e*i',...

    'y1(0)=0.0005-0.00013*i','y2(0)=-6.28-15*i');

a=0.1;

b=0.086;

g=0.045;

c=0.01;

d=0.8;

e=-0.069+0.109*i;

y1=subs(y1,{'a ' 'b' 'c' 'd' 'e' 'g'},{a b c d e g})

y2=subs(y2,{'a ' 'b' 'c' 'd' 'e' 'g'},{a b c d e g})

附图1.jpg

赞一下(1人)

回复此楼

MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。

5楼2014-02-20 13:54:27

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

普通回帖

smutao

禁虫 (著名写手)

感谢参与，应助指数 +1

本帖内容被屏蔽

2楼2014-02-18 18:04:55

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

suixin123456

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 342.4
红花: 47
帖子: 96
在线: 23.2小时
虫号: 1578539
注册: 2012-01-15
性别: MM
专业: 粒子物理学和场论

引用回帖:

3楼: Originally posted by 月只蓝 at 2014-02-18 18:24:54
程序如下：

function re_dsolve

a=0.1;
b=0.086;
g=0.045;
c=0.01;
d=0.8;
e=-0.069+0.109*i;

=dsolve('Dy1 = (-b+a*i)*y1+g*y2', 'Dy2 =(-d/2+c*i)*y2-g*y1+e*i')

结果的一部分见附图1。

...

结果的形式倒是对，不过abcd都赋值了为什么结果中还有符号？而且初始值都给了，还有c1c2？我就想要最后结果，不知道怎么弄，之前我也能算到带C1C2的，但是我需要用这个算几十组数据的最后结果

赞一下

回复此楼

精诚所至，金石为开

4楼2014-02-20 11:31:37

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

suixin123456

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 342.4
红花: 47
帖子: 96
在线: 23.2小时
虫号: 1578539
注册: 2012-01-15
性别: MM
专业: 粒子物理学和场论

引用回帖:

5楼: Originally posted by 月只蓝 at 2014-02-20 13:54:27
已加上初值条件，并把符号解转换为数值解，代码如下，结果的一部分见附图1。
function test3333
clear all;clc
syms a b c d e g

=dsolve('Dy1 = (-b+a*i)*y1+g*y2', 'Dy2 =(-d/2+c*i)*y2-g*y1+e*i',...
...

你这个太强大了，照你这个编辑，调用的时候怎么输入？

赞一下

回复此楼

精诚所至，金石为开

6楼2014-02-23 11:08:53

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