对一个复杂的代数方程的MATLAB图形绘制

[rocky14]

有这么个方程
b=2*atan((y.^2-x^2*1.483^2)./(x^2*1.49496^2-y.^2))./sqrt(x^2*1.49496^2-y.^2);
x=2*pi/1.550918;
b=0:0.002:8;
要绘制以b为变量，y为结果的图形。但是我无法化简，但matlab提示说因变量只能单独在一边，我没法了。
请大家帮我解决一下。谢谢了。

[laosam280]

应该是这样的吧？
b=2*atan((y.^2-x^2*1.483^2)./(x^2*1.49496^2-y.^2))./sqrt(x^2*1.49496^2-y.^2);
x=2*pi/1.550918;
y=0:0.002:8;
那你干脆用循环算了，将b的结果保存在一个矩阵里。
b=ones(length(y),1);
for i=1:400
y=0:0.02:8
x=2*pi/1.550918;
b(i)=2*atan((y(i)^2-x^2*1.483^2)./(x^2*1.49496^2-y(i)^2))./sqrt(x^2*1.49496^2-y(i)^2);
end
plot(b,y,'.-')

呵呵，大概就是这样吧

[kaklx]

不好，你的程序有错误！

[yindawei123]

你先找人把那个式子化简了，不化简有MATLAB很难弄出来的，至少我不会；
只要把那个式子化简了，就好做了

[rocky14]

就是想借用MATLAB的强大功能化简哈 ，哎 我的论文该怎么办呢

[laosam280]

这是一个分式方程，主要需要考虑分子的取值。分子含有周期函数，因此求解出来的结果也是周期性很强的。我画出了整体的解图像和单个周期的图像。在附件里。