你需要括号
2*n。
否则，Matlab将除以2，然后
乘
经过
n
。
百分比Chebshevy波兰人
Clc
清除
n = 3;％过滤器订单
波纹= 1;％ D b
epsilon = sqrt（10^（0.1*Ripple）-1）;
a =（1/n）*asinh（1/epsilon）;
为了k = 1：1：n
sigmak（k）= -sinh（a）*sin（（（（2*k -1）/（2*n））*pi）;
Omegak（k）= Cosh（a）*cos（（（（2*k-1）/（2*n））*pi）;
结尾
disp（sigmak）
disp（Omegak）
