Unable to convert the following expression into double array

Question

chaaru datta am 24 Jun. 2022

0
Verknüpfen

Direkter Link zu dieser Frage

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

Kommentiert: chaaru datta am 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 Kommentare
1 älteren Kommentar anzeigen1 älteren Kommentar ausblenden

chaaru datta am 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 am 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.

Melden Sie sich an, um zu kommentieren.

Melden Sie sich an, um diese Frage zu beantworten.

Answer 1

Torsten am 24 Jun. 2022

0
Verknüpfen

Direkter Link zu dieser Antwort

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

Bearbeitet: Torsten am 24 Jun. 2022

In MATLAB Online öffnen

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