Unable to convert the following expression into double array

Question

chaaru datta on 24 Jun 2022

0
Link

Direct link to this question

https://ch.mathworks.com/matlabcentral/answers/1747050-unable-to-convert-the-following-expression-into-double-array

Commented: chaaru datta on 25 Jun 2022

Hello all, When I run my code, I keep getting this error: "Unable to convert expression containing remaining symbolic function calls into double array. Argument must be expression that evaluates to number."

clc;
clear all;
close all;
al = 0.5; 
n = 5; 
sig2_fsnrn = 1; 
sig2_ksln = 1; 
sig2_glnrn = 10; 
sig2_e_fsnrn = 0;
sig2_e_ksln= 0; 
sig2_e_glnrn =0; 
L = 1;
Nr = 1; 
gamma_th_dB = 3;
gamma_th = (10^(gamma_th_dB/10));
 vel_r =0; 
 c = 3*10^8; 
 fc = 915*10^6; 
 Rs = 9.6*10^3; 
 Z1 = ((2*pi*fc*vel_r)/(Rs*c)); 
 P_fsnrn = besselj(0,Z1); 
 P_ksln = besselj(0,Z1);
 P_glnrn = besselj(0,Z1); 
 
 v_phi_fsnrn = (1-P_fsnrn^(2*(n-1)))*sig2_fsnrn; 
 v_phi_ksln = (1-P_ksln^(2*(n-1)))*sig2_ksln; 
 v_phi_glnrn  = (1-P_glnrn ^(2*(n-1)))*sig2_glnrn ;
 
 rho2_fsnrn = (P_fsnrn)^(2*(n-1));
 rho2_ksln = (P_ksln)^(2*(n-1));
 rho2_glnrn = (P_glnrn)^(2*(n-1));
 
 
  OP_th = [];
  
  for j1 = -5:5:45 
     jj1 = 10^(j1/10); 
     jj1
        
     
      mu4_1 = (al^2)*v_phi_glnrn*(v_phi_ksln + (rho2_ksln*sig2_e_ksln));
      mu4_2 = (al^2)*rho2_glnrn*sig2_e_glnrn*(v_phi_ksln + (rho2_ksln*sig2_e_ksln));
      mu4_3 =  v_phi_fsnrn+(rho2_fsnrn*sig2_e_fsnrn+(1/jj1));
      mu4 = mu4_1+mu4_2+mu4_3;
      
      zeta = (al^2)*rho2_glnrn*rho2_ksln;
      mu1 = (al^2)*rho2_glnrn*(v_phi_ksln + (rho2_ksln*sig2_e_ksln));
      mu2 = (al^2)*rho2_ksln*(v_phi_glnrn + (rho2_glnrn*sig2_e_glnrn));
      mu3 = rho2_fsnrn;
     
       syms x1;
     part1 = exp((gamma_th*mu2*x1)/((zeta*x1-mu1)*sig2_glnrn));
     part2 = exp((gamma_th*mu4)/((zeta*x1-mu1)*sig2_glnrn));
     part3 = exp(-x1/sig2_ksln);
     part4_1 = (-1/sig2_fsnrn);
     part4_2 = (gamma_th*mu3)/((zeta*x1-mu1)*sig2_glnrn);
     part4 = -1/(part4_1+part4_2);
     
     pt_1 = part1*part2*part3*part4;
     I1 = int(pt_1,0,10);
     op_th1 = (1/(sig2_ksln*sig2_fsnrn))*I1;
     op_th = 1-op_th1;
     OP_th = [OP_th,op_th];            
  end
  SNRdB = -5:5:45; 
grid on;
semilogy(SNRdB,OP_th,'g-','LineWidth',1.1);

3 Comments
Show 1 older commentHide 1 older comment

chaaru datta on 24 Jun 2022

Thanks for your answer sir. It means MATLAB is unable to solve the integration. So does it mean that we cannot plot this expression in MATLAB.

KSSV on 24 Jun 2022

In a way yes.....try to find out the integration of the expression. Or read tips to see whether integration can be solved in matlab.

Sign in to comment.

Sign in to answer this question.

Answer 1

Torsten on 24 Jun 2022

0
Link

Direct link to this answer

https://ch.mathworks.com/matlabcentral/answers/1747050-unable-to-convert-the-following-expression-into-double-array#answer_992570

Edited: Torsten on 24 Jun 2022

Open in MATLAB Online

clc;

clear all;

close all;

al = 0.5;

n = 5;

sig2_fsnrn = 1;

sig2_ksln = 1;

sig2_glnrn = 10;

sig2_e_fsnrn = 0;

sig2_e_ksln= 0;

sig2_e_glnrn =0;

L = 1;

Nr = 1;

gamma_th_dB = 3;

gamma_th = (10^(gamma_th_dB/10));

vel_r =0;

c = 3*10^8;

fc = 915*10^6;

Rs = 9.6*10^3;

Z1 = ((2*pi*fc*vel_r)/(Rs*c));

P_fsnrn = besselj(0,Z1);

P_ksln = besselj(0,Z1);

P_glnrn = besselj(0,Z1);

v_phi_fsnrn = (1-P_fsnrn^(2*(n-1)))*sig2_fsnrn;

v_phi_ksln = (1-P_ksln^(2*(n-1)))*sig2_ksln;

v_phi_glnrn = (1-P_glnrn ^(2*(n-1)))*sig2_glnrn ;

rho2_fsnrn = (P_fsnrn)^(2*(n-1));

rho2_ksln = (P_ksln)^(2*(n-1));

rho2_glnrn = (P_glnrn)^(2*(n-1));

OP_th = [];

for j1 = -5:5:45

jj1 = 10^(j1/10);

%jj1

mu4_1 = (al^2)*v_phi_glnrn*(v_phi_ksln + (rho2_ksln*sig2_e_ksln));

mu4_2 = (al^2)*rho2_glnrn*sig2_e_glnrn*(v_phi_ksln + (rho2_ksln*sig2_e_ksln));

mu4_3 = v_phi_fsnrn+(rho2_fsnrn*sig2_e_fsnrn+(1/jj1));

mu4 = mu4_1+mu4_2+mu4_3;

zeta = (al^2)*rho2_glnrn*rho2_ksln;

mu1 = (al^2)*rho2_glnrn*(v_phi_ksln + (rho2_ksln*sig2_e_ksln));

mu2 = (al^2)*rho2_ksln*(v_phi_glnrn + (rho2_glnrn*sig2_e_glnrn));

mu3 = rho2_fsnrn;

%syms x1;

part1 = @(x1)exp((gamma_th*mu2*x1)./((zeta*x1-mu1)*sig2_glnrn));

part2 = @(x1)exp((gamma_th*mu4)./((zeta*x1-mu1)*sig2_glnrn));

part3 = @(x1)exp(-x1/sig2_ksln);

part4_1 = (-1/sig2_fsnrn);

part4_2 = @(x1)(gamma_th*mu3)./((zeta*x1-mu1)*sig2_glnrn);

part4 =@(x1) -1./(part4_1+part4_2(x1));

pt_1 = @(x1)part1(x1).*part2(x1).*part3(x1).*part4(x1);

I1 = integral(pt_1,0,10)

op_th1 = (1/(sig2_ksln*sig2_fsnrn))*I1;

op_th = 1-op_th1;

OP_th = [OP_th,op_th];

end