untitled11.m

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

for m=6:15;
for j=1:m-2
a(j)=(j*cos((pi/6)-(j*pi)/(6*(m-1))));
end
sum=a(j);
a(1:j);
j=0;
while j<m-3
    j=j+1;
    sum=sum+a(j);
end
sum;
a=(m-1)/2+sum;
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