Can anyone please help how to simulate closed form expression in Matlab?

Question

Dhawal Beohar 2020-11-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/642290-can-anyone-please-help-how-to-simulate-closed-form-expression-in-matlab

评论： Dhawal Beohar 2020-11-9

采纳的回答： Walter Roberson

How to simulate closed form expression in Matlab?

5 个评论
显示 3更早的评论隐藏 3更早的评论

Dhawal Beohar 2020-11-9

在 MATLAB Online 中打开

Hi Walter! I am sharing my code below.

I want to calculate closed form for mu_o_fix. I have solved this mathematically and got equation (1) as final expression. Can please anyone help how to perform the simulation for closed form expression.

clc
clear all;
%Parameters defined
 
N_0=10^-11.4;                                                               
Xsr=1;                                                                      
Bmax=50;                                                                    
Bslot=3;                                                                
PmaxdB=[15];                                                                
Pmax=10.^(PmaxdB/10);                                                       
 
Qu=[-3:0.1:-1];
QoS_Component_u=10.^(Qu);
% Starting For Loop for QoS Component u  
  
for m=1:length(QoS_Component_u)
    u=QoS_Component_u(m);
% Starting For Loop for Expectation
    
    for n=1:100000                                                          % First Monte-Carlo Simulations  
 
% Equations for Rayleigh Fading h1
        hx1=(1/sqrt(2))*randn(1,1);
        hy1=(1/sqrt(2))*randn(1,1);
        h1=hx1.^2+hy1.^2; 
% Ending Rayleigh Fading h1                      
    end
    
% Ending Loop for Expectation 
        ee(m)=(-u*((Pmax*h1)+N_0)-h1);                                      %exponent expression
        mu_o_fix(m) = ((1/u)*log(ee(m))/(Pmax*u)+1);                   % equation (1)
        mu_o_fix_snr(m) = mu_o_fix(m)/N_0;
        mu_o_fix_dB(m)=10*log(mu_o_fix_snr(m))/log(10);
        mu_o_fix_dB_real(m) = real(mu_o_fix_dB(m));
end

Dhawal Beohar 2020-11-9

equation 1 is the closed expression i derived and now i want to simulate in Matlab. I don't know if this is the correct way i am doing.

Dhawal Beohar 2020-11-9

编辑：Dhawal Beohar 2020-11-9

i want to calculate closed form expression for the following equation:

; where

in Matlab.

请先登录，再进行评论。

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

Answer 1

Walter Roberson 2020-11-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/642290-can-anyone-please-help-how-to-simulate-closed-form-expression-in-matlab#answer_539345

在 MATLAB Online 中打开

Ten = sym(10);
N_0 = Ten^sym(-11.4);
PmaxdB = sym(15);
Pmax = Ten.^(PmaxdB/Ten);                                                       
syms gamma Qu
assume(Qu < 0)
QoS_Component_u = Ten.^(Qu);
u = QoS_Component_u;
Q = -1/u * log(int(exp(-u*(Pmax*gamma + N_0))*exp(-gamma), gamma, 0, inf));
Qs = simplify(Q);
disp(char(Qs))
(log(10*10^(1/2)*10^Qu + 1) + (10^(3/5)*10^Qu)/1000000000000)/10^Qu
disp(char(vpa(Qs,5)))
(log(31.623*10.0^Qu + 1.0) + 3.9811e-12*10.0^Qu)/10.0^(1.0*Qu)