话务量计算编程实现（MATLAB）

</pre><pre code_snippet_id="269865" snippet_file_name="blog_20140402_1_5289452" class="plain" name="code">Matlab

clc;clear;
%话务量计算程序：
n = input('Please Input the number of channels:\n'); %输入信道数n
b = input('Please Input the GOS:\n'); %输入呼损率b
k = 1;
num =1;
syms A;

%公式中的分母表达
for i= 1:n
k = k*i; %计算i的阶乘：i!,个人认为不比Matlab中的阶乘函数麻烦，所以没选择阶乘函数
num = num + (A^i)/k;
end
%公式中的分子表达
den = (A^n)/k;

%呼损率公式
B = den/num;
%解方程求呼损率Aa
Aa = solve(B -b); %方程的解可能有多个值

m =length(Aa)

for j=1:m
if subs(Aa(j))>0 && imag(subs(Aa(j)))==0 %判断Aa的值是否大于0，且不为虚数
Aa(j)
end
end