使用符号方法的Matlab方程

时间:2018-02-12 22:01:29

标签: matlab

我试图在Matlab中使用符号方法来解决以下等式。

我的最终结果(见下面的变量f6)再次以sym变量的形式出现。 任何人都可以建议我如何以双打格式得到答案。我已经尝试过' vpa'和'简化'命令。我正在运行的代码也在下面给出。

enter image description here

f6 = (3884161996073403*symsum((5^jj*symsum((-1)^(jj - m)/factorial(jj - m), m, 0, jj))/factorial(jj), jj, 0, Inf))/576460752303423488

悲伤

clc;
  clear all;
  close all;
Pc = 45-30;
Pc_lin = 10^(Pc/10);
Po = 3-30;
Po_lin = 10^(Po/10);
lambda_m = 6e-6
alpha_l = 2;
alpha_o = 4;
varepsilon_2 = 0.5;
Xd = 5;
% T_mj = e^(-K)*K^j*(1/(factorial(j) * factorial(j-m)) )
beta = (2*pi/alpha_o)/sin(2*pi/alpha_o);
syms jj m
K = 5; 
    T = -5;
    T_lin = 10^(T/10);
    Cd = pi*(T_lin*varepsilon_2*Xd^alpha_l)^(2/alpha_o)*(lambda_m*(Pc_lin/Po_lin)^(2/alpha_o))*beta;
          A = Cd*T_lin^(2/alpha_o);
    B = exp(-1*Cd*T_lin^(2/alpha_o));         
    % T1
    syms jj m z ii
    f1 = (-1)^(jj-m)/(factorial(z-ii)*factorial(ii)) * gamma(2*ii/alpha_o)/gamma(2*ii/alpha_o-(jj-m)+1)
    f2 = symsum(f1, ii, 1, z)
    f3 = symsum(f2,z,1,jj-m);
    f4 = (-1)^(jj-m)*exp(-K)*K^jj/(factorial(jj)*factorial(jj-m))
    f5 = symsum(f4, m, 0, jj)
    f6 = symsum(f5, jj,0,inf)

0 个答案:

没有答案