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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】求达人帮我解一阶常微分方程组
21/1返回列表
查看: 495  |  回复: 1
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fan_keke

铜虫 (初入文坛)

[交流] 【求助】求达人帮我解一阶常微分方程组

求达人帮我解一阶常微分方程组,万分的感谢！

dy1/dt= - i*( - u*y1+J*y2+K*y3);
dy2/dt= - i*(J*y1+J*y3);
dy3/dt= - i*(u*y3+J*y2+K*y1);
其中 u, J, K 均为常数，i 为虚数，y1, y2, y3 满足|y1|^2+|y2|^2+|y3|^2=1,
初始条件为:y1(0)=0, y2(0)=0, y3(0)=1, 求 y1(t), y2(t), y3(t) 的解析表达式。

[ Last edited by fan_keke on 2010-10-11 at 08:47 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-10-10 23:26:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijinfeng042

木虫 (小有名气)

Matlab
★ ★
nono2009(金币+2):鼓励应助。 2010-10-15 06:24:17
fan_keke(金币+10):我还在想办法，谢了 2010-10-26 08:47:28
引用回帖:
Originally posted by fan_keke at 2010-10-10 23:26:04:
求达人帮我解一阶常微分方程组,万分的感谢！

dy1/dt= - i*( - u*y1+J*y2+K*y3);
dy2/dt= - i*(J*y1+J*y3);
dy3/dt= - i*(u*y3+J*y2+K*y1);
其中 u, J, K 均为常数，i 为虚数，y1, y2, y3 满足|y1|^2+|y2|^2 ...

一般是没有表达式的
CODE:
function caogao
Y0=[0 0 1];%初值
tspan=0:10;%求解区间
u=1;J=2;K=4; %已知的值
[t,Y]=ode45(@func,tspan,Y0,[],u,J,K);
Y2=sum(abs(Y),2);%条件|y1|^2+|y2|^2+|y3|^2=1
disp(['满足条件的解为 ' num2str(Y(Y2-1<1e-9,:))]) %找出满足的解
function dydt=func(t,y,u,J,K)
dydt=zeros(3,1);
dydt(1)= - 1i*( - u*y(1)+J*y(2)+K*y(3));
dydt(2)= - 1i*(J*y(1)+J*y(2));
dydt(3)= - 1i*(u*y(3)+J*y(2)+K*y(1));
赞一下(1人) 回复此楼
工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
2楼2010-10-14 08:07:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fan_keke 的主题更新
21/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 XJXMuRjD7cK7 2026-11-13 3/150 2026-11-13 23:37 by 6PrO1XyZ5URF
[论文投稿] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 gRND4NVAZ8VG 2026-11-11 5/250 2026-11-13 11:15 by Equinoxhua
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:13 by Equinoxhua
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 by Equinoxhua
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +5 C6CRn66RYYXS 2026-11-10 6/300 2026-11-13 11:04 by Equinoxhua
[考博] 27博士 +5 锦鲤送福至 2026-11-10 5/250 2026-11-13 06:22 by wlf1982
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +4 Z7qfCsdIiD1T 2026-11-11 5/250 2026-11-12 11:44 by HajghKOMynmz
[基金申请] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 Z7qfCsdIiD1T 2026-11-10 9/450 2026-11-12 11:22 by HajghKOMynmz
[考博] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +4 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:11 by HajghKOMynmz
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 10:38 by HajghKOMynmz
[硕博家园] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +8 cqv1zKnzd7hr 2026-11-09 12/600 2026-11-12 09:32 by HajghKOMynmz
[基金申请] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-12 09:21 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 czAiolzR0IJC 2026-11-10 5/250 2026-11-12 06:36 by 9fOTqj074tIg
[硕博家园] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-12 06:14 by 9fOTqj074tIg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 05:41 by 9fOTqj074tIg
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 YHiw1gou1Eda 2026-11-10 7/350 2026-11-12 05:39 by 9fOTqj074tIg
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 4/200 2026-11-12 02:23 by 9fOTqj074tIg
[文学芳草园] 厉娜 +3 张丹枫1991 2026-11-09 5/250 2026-11-09 16:53 by 张丹枫1991
信息提示
关闭
请填处理意见
关闭