Numerical Integration at different point

Pritam Sahoo

2018 7 5

1 个回答

回答已采纳

4 次查看（30 天）

显示更早的评论

0 个投票

syms g;

I0_dBm=-12:22;

for jj=1:length(I0_dBm)

    var_ni(jj)=(6.09*10^-10*(10^(I0_dBm(jj)/10))*g)+(8.28*10^-14);
    fun(jj)=(1/4)*(2.02*10^7)*1.3844*(g^(2.1417))*exp(-(10* (10^(I0_dBm(jj)/10))*(10^-3)*g/(2*sqrt(var_ni(jj))))^2)*erfc((log(g/(0.0068))+0.1821)/sqrt(0.1));
    ber(jj)=sym(int(fun(jj),g, 0, Inf));

end

figure

semilogy(I0_dBm,ber,'g-d','LineWidth',2);

I tried this integration but the error is-

Error using semilogyConversion to double from sym is not possible.*

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

请先登录，再进行评论。

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

请先登录再关注

采纳的回答

Torsten 2018-7-5

在 MATLAB Online 中打开

0 个投票

I0_dBm=-12:22;
for jj=1:length(I0_dBm)
    var_ni=@(g)(6.09*10^-10*(10^(I0_dBm(jj)/10))*g)+(8.28*10^-14);
    fun=@(g)(1/4)*(2.02*10^7)*1.3844*(g.^(2.1417)).*exp(-(10* (10^(I0_dBm(jj)/10))*(10^-3)*g./(2*sqrt(var_ni(g)))).^2).*erfc((log(g/(0.0068))+0.1821)/sqrt(0.1));
    ber(jj)=integral(fun, 0, Inf);
end
figure
semilogy(I0_dBm,ber,'g-d','LineWidth',2);