24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

wzp0098

木虫 (著名写手)

应助: 2 (幼儿园)
金币: 2524.4
散金: 39
红花: 13
帖子: 1218
在线: 226.6小时
虫号: 853822
注册: 2009-09-22
性别: GG
专业: 光学

[求助] 用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以上，很难精确求解解析解）
请大侠指教！

回复此楼

» 猜你喜欢

招博士已经有4人回复
救命帖已经有7人回复
限项规定已经有6人回复
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向已经有3人回复
英文综述是否需要润色及查重已经有5人回复
为什么nbs上溴没有产物点出现呢已经有9人回复
最失望的一年已经有18人回复

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

求助，有关于matlab作图的问题已经有7人回复
matlab画图问题！！！急！希望会的大神来看看！！已经有11人回复
如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！已经有11人回复
matlab画图拷到word后，想快速更改坐标轴标识已经有3人回复
matlab作图问题已经有9人回复
请问用matlab怎样画下面这幅图? 已经有10人回复
matlab画图问题: 怎样用小图标或自定义图形代替+ x . o等做图已经有10人回复
用matlab怎么绘制z=xy的图像已经有5人回复
有达人愿意共享自己论文写作时的matlab绘图程序吗已经有10人回复
matlab作图问题已经有5人回复
matlab 画图求助已经有10人回复
MATLAB绘图-直线不平滑已经有6人回复
matlab中遇到的画图小问题，程序没问题。标题要长@@@@@@@@@@@@@@@@@@@@@@ 已经有8人回复
matlab（unix版）画图，内存溢出已经有4人回复
如何用Matlab画矩阵的散点图？已经有3人回复
matlab画图求助已经有11人回复
继续求助Matlab绘图！已经有8人回复
关于Matlab作图的问题已经有14人回复
Matlab中如何调用M文件作图已经有9人回复
求助matlab三维画图中：Inner matrix dimensions must agree. 已经有4人回复
MATLAB 画图小问题已经有11人回复
求助：谁会用matlab画图？已经有13人回复
【求助完成】matlab画图surf 已经有5人回复
【求助】matlab公式输入及画图已经有6人回复

1楼 2012-08-14 11:15:07

已阅关注TA 给TA发消息送TA红花 TA的回帖

wzp0098

木虫 (著名写手)

应助: 2 (幼儿园)
金币: 2524.4
散金: 39
红花: 13
帖子: 1218
在线: 226.6小时
虫号: 853822
注册: 2009-09-22
性别: GG
专业: 光学

引用回帖:

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的回帖

查看全部 7 个回答

csgt0

荣誉版主 (著名写手)

彩色挂图

专家经验: +2
应助: 367 (硕士)
贵宾: 0.457
金币: 7807
散金: 8897
红花: 20
沙发: 7
帖子: 2809
在线: 515.8小时
虫号: 449046
注册: 2007-11-02
专业: 金属非晶态、准晶和纳米晶
管辖: 计算模拟

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与，应助指数 +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