24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 651  |  回复: 6

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以上,很难精确求解解析解)
请大侠指教!
回复此楼

» 猜你喜欢

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

已阅   关注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

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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 华南理工0703化学,总分336求调剂 +10 手机用户 2026-04-03 11/550 2026-04-06 22:20 by michael2011
[考研] 308求调剂 +13 倘若起风了呢 2026-04-05 13/650 2026-04-06 14:20 by 蒋皓禹
[考研] 290求调剂085701 +11 1314捧花 2026-04-02 11/550 2026-04-06 13:31 by zllcz
[考研] 085410人工智能 初试316分 求调剂 +7 残星拂曙 2026-03-31 7/350 2026-04-06 10:09 by 蓝云思雨
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +8 Mr. Z 2026-04-05 8/400 2026-04-06 09:24 by dongzh2009
[考研] 生物与医药调剂 +4 十七sa 2026-04-05 4/200 2026-04-05 20:05 by lys0704
[考研] 计算机11408,286分求调剂 +7 木子念晞 2026-04-05 7/350 2026-04-05 19:02 by chy09050039
[考研] 考研生物学考A区211,初试322,科目生化和生物综合,求调剂 +6 。。。54 2026-04-03 6/300 2026-04-05 14:54 by JOKER0401
[考研] 085602调剂 初试总分335 +12 19123253302 2026-04-04 12/600 2026-04-05 08:08 by 544594351
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 22:08 by 啵啵啵0119
[考研] 333求调剂 +12 wfh030413@ 2026-04-03 13/650 2026-04-04 21:02 by jj987
[考研] 334求调剂 +8 曾仰之 2026-04-03 8/400 2026-04-04 11:16 by w_xuqing
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 281求调剂 +10 aaawhy 2026-04-03 10/500 2026-04-03 21:42 by lbsjt
[考研] 化工调剂303分,过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 293求调剂 +5 末未mm 2026-04-02 6/300 2026-04-03 15:20 by 王保杰33
[考研] 312 化工或制药调剂 +8 小小墨123 2026-04-02 9/450 2026-04-03 09:12 by zhouxiaoyu
[考研] 321求调剂 一志愿 浙江工业大学生物医药 +5 嘿嘿HC 2026-04-01 6/300 2026-04-02 15:23 by sophie2180
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 调剂申请 +8 张张张张zy 2026-03-31 9/450 2026-04-01 08:29 by zjbkx
信息提示
请填处理意见