24小时热门版块排行榜    

查看: 373  |  回复: 1

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 ]
回复此楼
已阅   回复此楼   关注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));

工作了,偶尔会上来~可以关注新浪微博 @云是风的梦_Matlab
2楼2010-10-14 08:07:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fan_keke 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856化工专硕求调剂 +8 董boxing 2026-03-01 8/400 2026-03-01 16:49 by caszguilin
[考研] 材料学硕318求调剂 +8 February_Feb 2026-03-01 10/500 2026-03-01 16:48 by sunny81
[考研] 285求调剂 +8 满头大汗的学生 2026-02-28 8/400 2026-03-01 16:47 by caszguilin
[基金申请] 刚录用,没有期刊号,但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 295求调剂 +6 19171856320 2026-02-28 6/300 2026-03-01 15:52 by jxstnuZYX
[考研] 307求调剂 +5 wyyyqx 2026-03-01 5/250 2026-03-01 15:21 by Fff-1
[考研] 281求调剂 +4 2026计算机_诚心 2026-03-01 4/200 2026-03-01 15:21 by chemisry
[考研] 材料工程274求调剂 +3 Lilithan 2026-03-01 3/150 2026-03-01 14:58 by ms629
[考研] 求调剂 +6 repeatt?t 2026-02-28 6/300 2026-03-01 14:37 by Sakura绘
[考研] 291分工科求调剂 +8 science饿饿 2026-03-01 9/450 2026-03-01 14:22 by Ducount.Y
[考研] 材料学调剂 +8 提神豆沙包 2026-02-28 10/500 2026-03-01 14:14 by peike
[考研] 290求调剂 +8 材料专硕调剂; 2026-02-28 9/450 2026-03-01 12:46 by 闭眼看蓝天
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +3 虫虫虫虫虫7 2026-03-01 6/300 2026-03-01 11:50 by gaoxiaoniuma
[考研] 调剂 +3 简木ChuFront 2026-02-28 3/150 2026-03-01 11:46 by 王伟要上岸啊
[考博] 博士自荐 +4 kkluvs 2026-02-28 4/200 2026-03-01 10:19 by 馥安馥安
[硕博家园] 2025届双非化工硕士毕业,申博 +3 更多的是 2026-02-27 4/200 2026-03-01 10:04 by ztg729
[考研] 298求调剂 +5 axyz3 2026-02-28 5/250 2026-03-01 06:45 by 刘兵
[基金申请] 面上模板改不了页边距吧? +5 ieewxg 2026-02-25 6/300 2026-03-01 00:10 by addressing
[考研] 276求调剂 +3 路lyh123 2026-02-28 4/200 2026-02-28 19:45 by 路lyh123
信息提示
请填处理意见