【求助】新手求教关于MATLAB 解超越方程和画图 - 仿真模拟 - MATLAB

Truthfulness,Compassion,Tolerance,Beauty

11楼2010-05-01 15:25:18

brqhl_ing

银虫 (小有名气)

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

楼上的 x的取值是0到10 m取0,1,2,3,4
但运行提示 Imaginary parts of complex X and/or Y arguments ignored.
画出来的图好像不对。我的m是分立取值的不是循环

12楼2010-05-01 20:43:59

brqhl_ing

银虫 (小有名气)

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

引用回帖:

Originally posted by brqhl_ing at 2010-05-01 20:43:59:
楼上的 x的取值是0到10 m取0,1,2,3,4
但运行提示 Imaginary parts of complex X and/or Y arguments ignored.
画出来的图好像不对。我的m是分立取值的不是循环

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

13楼2010-05-01 20:49:45

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，把六个曲线画在一张图中，是楼主的原意吗？
还是要画六张图？

Truthfulness,Compassion,Tolerance,Beauty

14楼2010-05-02 00:12:19

robert2020

银虫 (著名写手)

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

★ ★
nono2009(金币+2):Nice. 2010-05-02 08:37

引用回帖:

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);

如果是六张图，程序可修改如下：

当 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);
figure(m)
plot(x,y)
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);
figure(m)
plot(x,y)
end

赞一下(1人)

Truthfulness,Compassion,Tolerance,Beauty

15楼2010-05-02 00:15:05

brqhl_ing

银虫 (小有名气)

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

引用回帖:

Originally posted by robert2020 at 2010-05-02 00:15:05:

如果是六张图，程序可修改如下：

当 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);
figure(m ...

m从0 1 2 3 4 取五个值利用for循环可以汇在一个图上吧?x的取值可以是方程有虚数根。我做的的光波导里的东西。那个K是可以为虚数的。不胜感激！！！

16楼2010-05-02 08:29:51