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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » Matlab画隐函数，ezplot
41/1返回列表
查看: 2134  |  回复: 3
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

hanye0517

木虫 (小有名气)

[求助] Matlab画隐函数，ezplot

用Matlab画粒子相图，准备用ezplot，当把公式'c_1*omega_k^2/2+c_2*(sin(phi)-phi*cos(phi_s))-H_phi=0'中的参数c_1，c_2,phi_s,H_phi经由计算得出的数值重新写成方程ezplot（‘2.7778*y^2/2+19.5695*(sin(x)-x*cos(-5*2*pi/360))-0.0043’）[直接用这一句就可以]时，可以得到相图.
但是，在直接带入，亦即用下面的程序时，画不出来。请各位指教！
主要有四个问题：
1、当变量不是x-y而是phi-omega_k时，怎么画图，貌似ezplot只能画xy平面。
2、为什么ezplot内的方程参数为数值时可以画图，经由计算得出参数带入就画不了呢？
3、我想解这个隐式方程的话，看到有教程说@（x），这是什么意思？
4、作相图经常要用到隐函数，如何模拟画图，有没有其他的方法呢？
谢谢

clear
clear all
T=1;%渡越时间因子
e=1.6*10^-19;%C
m0=9.1*10^-31;%kg
W0=0.511*10^6*e;%J
Wk=1*10^6*e;% 动能
c=2.99792458*10^8;
gamma_s=(W0+Wk)/W0;
beta_s=sqrt(1-(1/gamma_s)^2);
v=c;
f=2856*10^6;                                                   %Hz
lambda=v/f;
E0=10*10^6;                                                     % V
c_1=2*pi/(beta_s^3*gamma_s^3*lambda);
c_2=e*E0*T/W0;
phi_s=-5*2*pi/360;
H_phi=c_2*(sin(-phi_s)-(-phi_s)*cos(phi_s));
phi=(-360:360)*2*pi/360;
figure
ezplot('c_1*omega_k^2/2+c_2*(sin(phi)-phi*cos(phi_s))-H_phi')
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
道之所在，虽千万人吾往矣。
1楼2013-02-06 17:38:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

脱皮的洋葱

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
hanye0517: 金币+10, ★★★★★最佳答案 2013-02-08 15:12:02
csgt0: 金币+2, 谢谢 2013-02-16 09:19:12
1，虽然 matlab Error 里说 “c_1*omega_k^2/2+c_2*(sin(phi)-phi*cos(phi_s))-H_phi cannot be plotted in the xy-plane.” 这并不是楼主理解的那个意思。。xy 平面其实就是说二维平面， x 轴可以表示 phi 值，y 轴可以表示 omega_k 值。。。

2，这个问题我想是不是 ezplot 里单引号下的函数变成了字符串 所以无法识别字母所代表的数值

3，fun = @(x)x^2+2 是定义一个自变量为 x 的函数，以后再调用时可以直接用名字 fun 来调用 x^2+2

4，隐函数画图还有其他方法的 比如利用 contour 函数

5，楼主给出了 phi 值的变化范围，但为给出 omega_k 的取值范围，是无法画图的。。。。
赞一下 回复此楼
2楼2013-02-08 01:20:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hanye0517

木虫 (小有名气)
引用回帖:
2楼: Originally posted by 脱皮的洋葱 at 2013-02-08 01:20:57
1，虽然 matlab Error 里说 “c_1*omega_k^2/2+c_2*(sin(phi)-phi*cos(phi_s))-H_phi cannot be plotted in the xy-plane.” 这并不是楼主理解的那个意思。。xy 平面其实就是说二维平面， x 轴可以表示 phi 值，y 轴 ...

谢谢你，回答的很详细，已经找到问题了，就是没有定义变量xy，在clear后加syms x y 就好，想在问一下能不能把横纵坐标改成phi和omega_k?
赞一下 回复此楼
道之所在，虽千万人吾往矣。
3楼2013-02-08 15:11:47
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

脱皮的洋葱

木虫 (小有名气)
引用回帖:
3楼: Originally posted by hanye0517 at 2013-02-08 15:11:47
谢谢你，回答的很详细，已经找到问题了，就是没有定义变量xy，在clear后加syms x y 就好，想在问一下能不能把横纵坐标改成phi和omega_k?...

楼主可以试试

xlabel('\phi')
ylabel('\omega_k')
赞一下 回复此楼
4楼2013-02-08 16:23:39
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hanye0517 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭