Пример кода приведен ниже:
Код: Выделить всё
clc; clear; close all;
SNR=10;
sko_=sqrt((1/(2*10^(SNR/10))));
sko=sym(sko_);
A=sym(1);
M=sym(4);
m=sym(4);
syms x;
L=0;
R=2*pi;
f=exp(((abs(A)^2)/(2*sko^2))*cos(x-2*pi*m/M))/(2*pi*besseli(0,(abs(A)^2)/sko^2));
I1=int(f,'x',L,R);
disp(strcat('рез-т интегрирования:'));
disp(vpa(I1,10));