有没高手能帮我修改一下 MATLAB 数字信号处理程序中的错误

n=20;
rp=0.5;
wn=[100 8000]/10000
[b,a]=cheby1(n,rp,wn);
[h,w]=freqz(b,a,512);
h1=20*log10(abs(h));
m=1,b1=abs(b);
while fix(b1)>0
b1=abs(b)/(2^m);
m=m+1;
end
c1=fix(b1*2^6);
c=sign(b)*c1*2^(m-6-1);
n=1;b2=abs(a);
while fix(b2)>0
b2=abs(b)/(2^n);
n=n+1;
end
d1=fix(b2*2^6);
d=sign(a)*d1*2^(m-6-1);
[hq,w]=freqz(c,d,512);
h2=20*log10(abs(hq));
subplot(121);
plot(w/pi,h1,'b',w/pi,h2,'r');
grid;
axis([0 1 -100 100]);
xlabel('归一化频率');ylabel('幅频响应特性');
legend('量化前','量化后');title('相位量化前后幅频特性分析');
subplot(122);
[z1,p1,k1]=tf2zp(b,a);
[z2,p2,k2]=tf2zp(c,d);
zplaneplot([z1,z2],[p1,p2],[k1,k2],{'s','+','d','*'});
legend('量化前的零点','量化前的极点','量化后的零点','量化后的极点');
title('相位量化前后零极点分析');
c=sig

我写了一篇英语日记,有高手能帮我修改一下吗? 帮我把这段英文翻译过来！如果 有高手能帮我修改一下那就更好了 有没有人能帮我一下 有朋友能帮我修改一下公司简介吗？ 有日语高手没?能帮我翻译一下:我爱叶巧~!~!么 高手过来帮我修改一下这首诗啊 请高手帮我检查 修改一下！ 高手帮我修改一下这篇文章~ 做诗高手帮帮忙~帮我修改一下 高手请进！！！帮我修改一下！