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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 用matlab画图
71/1返回列表
查看: 694  |  回复: 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返回列表
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜？ +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果？ +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
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥？ +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 帮忙看看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
[基金申请] 什么时候出结果，有咨询渠道？？？ +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 确定了，国自然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
信息提示
关闭
请填处理意见
关闭