无法计算出均值，方差，均方值，求大神们指点迷津，感激不尽！！

matlabchen

%滤波器
%巴特沃斯模拟滤波器生成数字滤波器
fp=1000;fs=2000;
rp=1;rs=35;
Fs1=40000;
wp=2*pi*fp/Fs1;
ws=2*pi*fs/Fs1;
wap=tan(wp/2);
was=tan(ws/2);
Fs1=Fs1/Fs1;
[N,Wn]=buttord(wap,was,rp,rs,'s');
[z,p,k]=buttap(N);
[bp,ap]=zp2tf(z,p,k);%得到传输函数
[bs,as]=lp2lp(bp,ap,wap);%低通到低通，频谱变换
[bz,az]=bilinear(bs,as,Fs1/2); %将模拟滤波器传递函数转换为数字滤波器传递函数
[H,w]=freqz(bz,az,256,Fs1*10000);%特性分析
figure(11)
plot(w,abs(H));%画出数字滤波器的波形图
title('低通滤波器的频谱');
xlabel('f/hz');
grid on;
%%
%信号通过低通滤波器
a=filter(bz,az,x);
disp('输入a均值:');
a_junzhi=mean(a);%均值
disp('输入a方差:');
a_fangcha=var(a);%方差
disp('输入a均方值:');
a_junfangzhi=std(a);%均方值

happy

信号x没有定义

matlabchen

happy 发表于 2012-11-22 20:18

                               
登录/注册后可看大图

信号x没有定义

这是程序的一段，x前面定义过了
sine=sin(2*pi*100*t);%正弦信号
x=awgn(sine,3);%在正弦信号加入高斯白噪声

happy

matlabchen 发表于 2012-11-23 08:42

                               
登录/注册后可看大图

这是程序的一段，x前面定义过了
sine=sin(2*pi*100*t);%正弦信号
x=awgn(sine,3);%在正弦信号加入高斯白 ...

不是程序无法计算，是你的输出写的有问题，可修改为：

1. a=filter(bz,az,x);
   
2. a_junzhi=mean(a);%均值
   
3. disp(['输入a均值:',num2str(a_junzhi)]);
   
4. a_fangcha=var(a);%方差
   
5. disp(['输入a方差:',num2str(a_fangcha)]);
   
6. a_junfangzhi=std(a);%均方值
   
7. disp(['输入a均方值:',num2str(a_junfangzhi)]);

复制代码

matlabchen

happy 发表于 2012-11-23 09:08

                               
登录/注册后可看大图

不是程序无法计算，是你的输出写的有问题，可修改为：

非常感谢，可以了，能说说为什么吗？

happy

matlabchen 发表于 2012-11-23 09:25

                               
登录/注册后可看大图

非常感谢，可以了，能说说为什么吗？

这个问什么还真不知道怎么回答
简单来说说计算机太听话了，你没让他显示，他是坚决不会让你显示的

如果以后经常用到，建议一本入门的书籍系统的学习一下
如果只是现在临时用到，建议doc看一下帮助吧，比如disp命令应该怎么用