How can I plot the following function?

2023 12 10
0 个回答
更新时间：2023 12 10
2 次查看（30 天）
Ran in:

0 个投票

Ran in:
alpha=0.001:0.1:5;
s=1;sb=.1;s1=sb./s;j=0.1;EZ=1;zeta0=-1;beta=1;i=(-1).^(1./2);
sigma=0.1; ur=1;k=5;b1=0.1;b2=0.1;c=0.1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
k1=k.^2.*zeta0./(1+c);
T1=(1+c)./(2.*s.*(i.*alpha.^2.*j-4.*c));
T2=(4.*s.*c-i.*alpha.^2)./(1+c);
r11=(4.*s.*c+i.*alpha.^2.*(c-j.*s.*(1+c)))./(c.*(1+c));r22=alpha.^2.*s.*(4.*i.*c+alpha.^2.*j)./(c.*(1+c));
alpha11=((1./2).*(r11+(r11.^2-4.*r22).^(1./2))).^(1./2);alpha12=-((1./2).*(r11+(r11.^2-4.*r22).^(1./2))).^(1./2);
alpha22=((1./2).*(r11-(r11.^2-4.*r22).^(1./2))).^(1./2);alpha21=-((1./2).*(r11-(r11.^2-4.*r22).^(1./2))).^(1./2);
KK=-(k1.*(k.^2-s.*(4-i.*alpha.^2.*j.*c.^(-1))))./(k.^4-(r11).*k.^2+r22);
format long g
syms r
L3 = vpaintegral(r*besselk(1,alpha11*r),r,sigma,1);
w1 =2.*( -T1 .* alpha11 .* L3.*(alpha11 .^ 2 + T2));
plot(alpha,abs(w1),'-black','LineWidth',1.2);
ylabel('$w1(r)$','Interpreter','latex','FontSize',12,'FontName','TiemsNewRoman','FontWeight','Normal')
xlabel('$\alpha$','Interpreter','latex','FontSize',12,'FontName','TiemsNewRoman','FontWeight','Normal')
%%%%%%%%%%%

6 个评论
显示 4更早的评论 隐藏 4更早的评论

Torsten
Torsten 2023-12-10
I get a plot, but note that the elements of L3 and thus of w1 are complex. Thus you can plot abs(w1), real(w1) or imag(w1), but not w1 itself.
Torsten
Torsten 2023-12-10
So problem solved ?
Shreen El-Sapa
Shreen El-Sapa 2023-12-10
where? can you excute it? I couldnot get any thing.
Torsten
Torsten 2023-12-10
编辑：Torsten 2023-12-10
See above. I reduced the alpha resolution from 0.0001 to 0.1 to save computation time (and I think this resolution is sufficient).

请先登录，再进行评论。

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

回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Lighting, Transparency, and Shading 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by