24小时热门版块排行榜

返回列表

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

brqhl_ing

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 262.9
散金: 11
红花: 1
帖子: 115
在线: 19.5小时
虫号: 854693
注册: 2009-09-23
性别: GG
专业: 物理电子学

[交流] 【求助】新手求教关于MATLAB 解超越方程和画图已有3人参与

y=sqrt(0.75/(1-x^2))*atan(-sqrt((4*x^2-1)/(1-x^2)))+m*pi/2
求出X随Y 的变化图像

回复此楼

» 猜你喜欢

1楼 2010-05-01 09:02:22

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

robert2020

银虫 (著名写手)

仿真EPI: 2
应助: 0 (幼儿园)
贵宾: 1.307
金币: 292.3
散金: 9143
红花: 19
沙发: 6
帖子: 1950
在线: 620.1小时
虫号: 861676
注册: 2009-10-02

brqhl_ing(金币+2):非常感谢你！昨天你给我的那个程序，我自己弄了一下已经很接近了1 2010-05-02 08:27

引用回帖:

Originally posted by brqhl_ing at 2010-05-01 20:49:45:

另外这是正确的方程形式y=(1-x.^2).^(-1).*atan(-((4*x.^2-1).*(1-x.^2).^(-1)).^(1/2))+(m*pi/2)*(1-x.^2).^(-1);

不好意思，在编程时没有计算x的取值范围，致使了虚数的出现，经计算x的有效范围为 0.5< x < 1 或 - 1< x < - 0.5，lz原意程序是在实数范围内吧？

当 0.5< x < 1 时，程序如下：

for m=1:6
x=0.5:0.01:1;
y=(1-x.^2).^(-1).*atan(-((4*x.^2-1).*(1-x.^2).^(-1)).^(1/2))+(m*pi/2)*(1-x.^2).^(-1);
plot(x,y)
hold on
end

当 - 1< x < - 0.5 时，程序如下：

for m=1:6
x=-1:0.01:-0.5;
y=(1-x.^2).^(-1).*atan(-((4*x.^2-1).*(1-x.^2).^(-1)).^(1/2))+(m*pi/2)*(1-x.^2).^(-1);
plot(x,y)
hold on
end

程序中m取值为1 2 3 4 5 6，把六个曲线画在一张图中，是楼主的原意吗？
还是要画六张图？