| 网站首页 | 业界新闻 | 小组 | 威客 | 人才 | 下载频道 | 博客 | 代码贴 | 在线编程 | 编程论坛
欢迎加入我们,一同切磋技术
用户名:   
 
密 码:  
共有 601 人关注过本帖
标题:亲们 求围观 求火速围观 大神们都来指导指导 Matlab 吧·····
只看楼主 加入收藏
sugeladijq
Rank: 1
等 级:新手上路
帖 子:9
专家分:3
注 册:2011-4-1
结帖率:100%
收藏
 问题点数:0 回复次数:2 
亲们 求围观 求火速围观 大神们都来指导指导 Matlab 吧·····
%***************  生成输入信号  *****%
f = 2000                  %设定基频率为2000HZ
    Fs=44100;                 %设定采样频率Fs=44.1kHz
N=256;                    %取的样本点数N
n=0:N-1;                  %建立矩阵
t=n/Fs;                   %采样时间t
w = 2*pi*f;
T = 2*pi/w;                   %求出正弦函数的周期
x1 = sin(w*t);
x2 =awgn(x1,5,'measured');      %加入SNR为5dB的高斯白噪声得到输入信号x2
xn =x2-x1;                    %间接获得白噪声xn
%*****************   时域波形   *************%
figure(1);                  
plot(t,x1,'k-');          %正弦波信号时域波形
title('正弦信号时域波形');
xlabel('t'),ylabel('x1(t)');
hold on;
grid on;
figure(2);                 
plot(t,xn,'r-');        %高斯白噪声时域波形
title('高斯白噪声信号时域波形');
xlabel('t'),ylabel('xn(t)');
hold on;
grid on;
figure(3)               
plot(t,x2,'b-');       %输入信号时域波形
title('输入信号时域波形');
xlabel('t'),ylabel('x2(t)');
hold on;
grid on;
   
    %*****************  噪声的数字特性  ****************%
   
    %*****************    1.均值       *****************%
   
    E = mean(xn);      %求白噪声的均值
     figure(4);
    plot(t,E,'b-');          %白噪声均值的时域波形
    title('噪声均值波形');
    xlabel('t'),ylabel('E');
hold on;
grid on;
   
    %****************      2.方差       *****************%
    D = var(xn);
     figure(5);
    plot(t,D,'b-');          %白噪声方差的时域波形
    title('噪声方差波形');
    xlabel('t'),ylabel('D');
hold on;
grid on;  
   
    %****************     3.均方值      ****************%
   
    E2 = E+D;
    figure(6);
    plot(t,E2,'b-');          %白噪声均方值的时域波形
    title('噪声均方值波形');
    xlabel('t'),ylabel('E2');
hold on;
grid on;
     %***************   4.输入信号自相关函数    ****************%
     R=xcorr(xn);     %*******自相关函数*******%
     t1=(-length(xn)+1:length(xn)-1)/Fs;
     figure(7);
     plot(t1,R,'r-');          %白噪声自相关函数的时域波形
     title('噪声自相关波形');
     xlabel('t1'),ylabel('R');
hold on;
grid on;  
   
     
    %****************      5.概率密度      *****************%  这个多次尝试都做不出来
    figure()
   
   
   
   
    %****************       6.频谱       ****************%    这里编译之后显示有错误 可能要改一下
     figure(9)
     H_arr=fft(xn,N);      
     cm=abs(H_arr);
     fl=(0:length(H_arr)-1)'*44100/length(H_arr);
     plot(f,H_arr,'b-')
     title('线性信号的频谱')
     xlabel('f'),ylabel('S_x_i(f)');
     hold on;
     grid on
    %***************        7.功率谱密度      *****************%  还有这里 出图像之后是直线 不知道对不
    figure(10)
    X_arr=fft(R);      
    cm=abs(X_arr);
    fl=(0:length(X_arr)-1)'*44100/length(X_arr);
    plot(f,X_arr,'r-')
    title('线性信号的功率谱')
    xlabel('f'),ylabel('S_x_i(f)');
    hold on;
    grid on        
请大家帮帮忙吧 帮我把后边打注释的完善下  还有就是看看那些有疑问的地方 真心谢谢大家了
搜索更多相关主题的帖子: 样本 
2011-10-30 10:48
sugeladijq
Rank: 1
等 级:新手上路
帖 子:9
专家分:3
注 册:2011-4-1
收藏
得分:0 
沙发下 大家赶紧来吧

生活是美好的
2011-10-30 10:48
vipwdp
Rank: 2
等 级:论坛游民
帖 子:3
专家分:14
注 册:2011-12-12
收藏
得分:0 
路过。。学习下
2012-03-20 15:12
快速回复:亲们 求围观 求火速围观 大神们都来指导指导 Matlab 吧·····
数据加载中...
 
   



关于我们 | 广告合作 | 编程中国 | 清除Cookies | TOP | 手机版

编程中国 版权所有,并保留所有权利。
Powered by Discuz, Processed in 0.036751 second(s), 7 queries.
Copyright©2004-2024, BCCN.NET, All Rights Reserved