求用MATLAB设计程序

[N,Wc]=buttord()

%估算得到Butterworth低通滤波器的最小阶数N和3dB截止频率Wc

[a,b]=butter(N,Wc)%设计Butterworth低通滤波器

[h,f]=freqz()%求数字低通滤波器的频率响应

figure(2)% 打开窗口2

subplot(221)%图形显示分割窗口

plot(f,abs(h))%绘制Butterworth低通滤波器的幅频响应图

title(巴氏低通滤波器'')

grid%绘制带网格的图像

sf=filter(a,b,s)%叠加函数S经过低通滤波器以后的新函数

subplot(222)

plot(t,sf)%绘制叠加函数S经过低通滤波器以后的时域图形

xlabel('时间 (seconds)')

ylabel('时间按幅度')

SF=fft(sf,256)%对叠加函数S经过低通滤波器以后的新函数进行256点的基—2快速傅立叶变换

w= %新信号角频率

subplot(223)

plot())%绘制叠加函数S经过低通滤波器以后的频谱图

title('低通滤波后的频谱图')

sum=0

aver=0

for i=1:6

for j=1:6

sum=sum+a(i,j)

end

end

aver=sum/36

fprintf('所有元素之和为：%f\n',sum)

fprintf('所有元素之和的平均值为：%f\n',aver)