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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 用matlab画图
71/1返回列表
查看: 693  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wzp0098

木虫 (著名写手)

[求助] 用matlab画图

如下的方程组，
I*(K2-K3)*Y1-I*Y2*Y3-I*C*Y3*Exp[-I*Q]=0
(-r-I*K2)*Y2-I*C-I*Y3*Y1=0
(-r-I*K1)*Y3-I*8-I*Y2*Y1*Exp[I*Q]=0
Y1，Y2，Y3为变量；I为虚数单位；
K1=1, K2=3, r=1，C=4, Q=pi/3，
我想利用for循环语句画出Y1实部随K3（取值-20到20）变化的图（请用数值画图方法，因为有时候变量很多10以上，很难精确求解解析解）
请大侠指教！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-08-14 11:15:07
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +1
wzp0098: 金币+60, ★有帮助, 谢谢帮助！ 2012-08-14 14:46:34
xzhdty: 金币+2, 谢谢 2012-08-14 17:53:40
CODE:
function wzp

for k3=-20:20
    y0=ones(3,1);
    Y(:,k3+21)=fsolve(@(y)fun(y,k3),y0);
end
Y1=Y(1,:)
plot(-20:20,Y1)
end

function F=fun(y,k3)
k1=1;
k2=3;
r=1;
C=4;
Q=pi/3;
F=[(k2-k3)*i*y(1)-i*y(2)*y(3)-i*C*y(3)*exp(-i*Q);
    (-r-i*k2)*y(2)-i*C-i*y(3)*y(1);
    (-r-i*k1)*y(3)-i*8-i*y(2)*y(1)*exp(i*Q);];
end
赞一下(1人) 回复此楼
showmethemoney
2楼2012-08-14 13:33:15
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wzp0098

木虫 (著名写手)
引用回帖:
2楼: Originally posted by csgt0 at 2012-08-14 13:33:15
function wzp

for k3=-20:20
    y0=ones(3,1);
    Y(:,k3+21)=fsolve(@(y)fun(y,k3),y0);
end
Y1=Y(1,
plot(-20:20,Y1)
end

function F=fun(y,k3)
k1=1;
k2=3;
r=1;
C=4;
Q=pi/3;
F=;
end ...

向您请教一个问题，方程 Y(:,k3+21)=fsolve(@(y)fun(y,k3),y0);
中 k3+21是什么含义？
如果我想以步长为0.1画图，怎么控制？
赞一下 回复此楼
3楼2012-08-14 14:46:13
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wzp0098: 金币+20, ★有帮助, 非常谢谢！我发现正你说那样，变量多了图不对，能否有办法处理呢？ 2012-08-15 09:11:14
引用回帖:
3楼: Originally posted by wzp0098 at 2012-08-14 14:46:13
向您请教一个问题，方程 Y(:,k3+21)=fsolve(@(y)fun(y,k3),y0);
中 k3+21是什么含义？
如果我想以步长为0.1画图，怎么控制？...

如果以0.1画图，那么就用for k3=-20:0.1:20，但是如果太小可能解方程很费时间，另外如果参数过多也可能解的不对。

k3+21主要是作为Y的列数，我将每次的结果都保存下来，k3=-20时就是1，一直k3=20时到41。如果间隔0.1，你可以单独用个变量来计数，每个循环+1。
赞一下 回复此楼
showmethemoney
4楼2012-08-14 17:46:40
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wzp0098

木虫 (著名写手)
好像步长没法控制
赞一下 回复此楼
5楼2012-08-15 09:15:32
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
wzp0098: 金币+20, ★有帮助, 谢谢帮助！ 2012-08-15 09:27:04
引用回帖:
4楼: Originally posted by csgt0 at 2012-08-14 17:46:40
如果以0.1画图，那么就用for k3=-20:0.1:20，但是如果太小可能解方程很费时间，另外如果参数过多也可能解的不对。

k3+21主要是作为Y的列数，我将每次的结果都保存下来，k3=-20时就是1，一直k3=20时到41。如果间 ...

比较困难，你可以估计下大概与k3的关系而确定可能的解，然后设置合适的初值
赞一下 回复此楼
showmethemoney
6楼2012-08-15 09:22:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

wzp0098

木虫 (著名写手)
function wzp

for k3=-20:0.1:20
    y0=ones(3,1);
    Y(:,k3+21)=fsolve(@(y)fun(y,k3),y0);
end
Y1=Y(1,
plot(-20:0.1:20,Y1)
end

function F=fun(y,k3)
k1=1;
k2=3;
r=1;
C=4;
Q=pi/3;
F=[(k2-k3)*i*y(1)-i*y(2)*y(3)-i*C*y(3)*exp(-i*Q);
    (-r-i*k2)*y(2)-i*C-i*y(3)*y(1);
    (-r-i*k1)*y(3)-i*8-i*y(2)*y(1)*exp(i*Q);];
end
步长按你说法没办法控制，因为我涉及的变量很多（8个以上），没有办法精确求解解析解。
赞一下 回复此楼
7楼2012-08-15 09:26:37
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wzp0098 的主题更新
71/1返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode +8 cratir 2026-08-14 12/600 2026-08-15 18:08 by zyfgau
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +7 Ldrop2023 2026-08-13 7/350 2026-08-15 15:57 by mzhh2000
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 22/1100 2026-08-14 23:45 by Noways
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 奇怪，两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
信息提示
关闭
请填处理意见
关闭