How can I plot the following function?

5 次查看（过去 30 天）

Shreen El-Sapa 2023-12-10

评论： Shreen El-Sapa 2023-12-10

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

%%%%%%%%%%%

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

Shreen El-Sapa 2023-12-10

Thanks so much.

MATLAB Graphics 2-D and 3-D Plots Line Plots

在 Help Center 和 File Exchange 中查找有关 Line Plots 的更多信息

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

Translated by