I am unable to plot this mat lab program

1 次查看（过去 30 天）

Shashibhushan Sharma 2017-12-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/372730-i-am-unable-to-plot-this-mat-lab-program

编辑： Walter Roberson 2017-12-14

My Program included the triple integral:

clearvars; close all; clc;
Bi=[0.001 0.1:0.1:0.9 0.999];N0=0.01;n=0.7;Rth=1;
Omg_sr=1;Omg_rd=1;p1=5;p=10^(0.1*p1);N=3;
sop_th=zeros(1,length(Bi));
for a=1:length(Bi)
    b=Bi(a);
    k1=(1-b)*p;k2=N0/(k1*Omg_sr);k3=n*b*(1-b)*p;k4=n*b*N0;k5=(1-b)*N0;    
    fun=@(z,y,x) ((((k4*x+k5)./(k3*x)).*(exp(-(((y./z).*((k4*x+k5)./(k3*x)))+x/Omg_sr)))).*(N*((1-exp(-((z*k2)./(z+1)))).^(N-1)).*(exp(-k2*z)).*((k2./(z+1))+(1./(z+1).^2))));
    zmin=0;zmax=Inf;
    ymin=0;ymax=@(z)(2^Rth*(1+z)-1);
    xmin=0;xmax=100;
    sop_th(a)=(integral3(fun, zmin,zmax,ymin,ymax,xmin,xmax))^N;
end
semilogy(Bi,sop_th,'r>-');
grid on;hold on;

If any mistake in my program please help me.

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

Shashibhushan Sharma 2017-12-13

In first factor of ((k4*x+k5)./(k3*x)), if I remove the term './(k3*x)', then it is integrated. Otherwise, it creates the problem. In this, variables are x,y and z.

Shashibhushan Sharma 2017-12-14