plot symbolic function real part

Question

john_arle 2019-3-31

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/453558-plot-symbolic-function-real-part

评论： Star Strider 2019-3-31

采纳的回答： Star Strider

在 MATLAB Online 中打开

Hello there,

I would like to plot the real part of this symbolic function I have just solved.

Parameters 'a' and 'e' are varying between 5 and 9, and 0.05 and 0.15 respectively.

Unfortunately, it does not work.

Here you go with my code ( solution is a ratio of sum of Bessel functions of order zero).

clc
clear
set(groot,'DefaultAxesFontSize',14) 
set(groot,'DefaultLineLineWidth',1.5)
syms a e
A = [besselj(0,-1i*a*e), bessely(0,-1i*a*e); ...
    besselj(0,-1i*a), bessely(0,-1i*a)];
b = [-1; -1];
X = linsolve(A, b);
disp(X)
a = linspace(5,9);
e = linspace(0.05, 0.15);
plot(a, e, X)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2019-3-31

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/453558-plot-symbolic-function-real-part#answer_368321

在 MATLAB Online 中打开

You are plotting a function of two variables, so you need to use a plot that accommodates them, and that your function is symbolic.

Try this:

syms a e
A = [besselj(0,-1i*a*e), bessely(0,-1i*a*e); ...
    besselj(0,-1i*a), bessely(0,-1i*a)];
b = [-1; -1];
X = linsolve(A, b);
disp(X)
a = linspace(5,9);
e = linspace(0.05, 0.15);
fsurf(real(X), [5 9, 0.05 0.15])
xlabel('a')
ylabel('e')