untitled33.m

来自「基本的迭代求和办法」· M 代码 · 共 23 行

for m=6:15;
for j=1:m-2
a(j)=(j*cos((pi/6)-(j*pi)/(6*(m-1))));
end
%s=a(j);
b=a(1:j);
b=b';
S=sum(b);
%j=0;
%while j<m-3
    %j=j+1;
   % s=s+a(j);
%end
%s;
a=(m-1)/2+S;
b=a^2;
c=27*128*(sin(pi/(12*(m-1)))^2);
d1=(pi^2)*[(4+sqrt(3))*((m-1)^2)];
d2=1+(11-6*sqrt(3))/(13*(m-1)^2);
d=d1*d2;
THD=sqrt(d/(c*b)-1);
end
THD