随机实验报告1CJ.doc

程序 1N=1024;fs=1000;n=0:N-1;signal=rand(1,N);%randn(1,N),exprnd(2,1,N),raylrnd(2,1,N),chi2rnd(2,1,N),signal_mean=mean(signal);signal_var=var(signal);signal_corr=xcorr(signal,signal,unbiased);signal_density=unifpdf(signal,0,1);signal_power=fft(signal_corr);%[s,w]=periodogram(signal);[k1,n1]=ksdensity(signal);[k2,n2]=ksdensity(signal,function,cdf);figure(1)subplot(2,3,1);hist(signal);title(图);subplot(2,3,2);plot(signal);title(曲线);subplot(2,3,3);f=n*fs/N; %频率序列plot(abs(signal_power));title(曲线);subplot(2,3,4);plot(angle(signal_power));title(曲线);subplot(2,3,5);plot(1:2047,signal_corr);title(图);subplot(2,3,6);plot(n1,k1);title(密度);figure(2)plot(n2,k2);结果程序2N=1024;signal_1=rand(1,N);signal_2=rand(1,N);signal_3=rand(1,N);signal_4=rand(1,N);signal_5=rand(1,N);signal=signal_1+signal_2+signal_3+signal_4+signal_5;[k1,n1]=ksdensity(signal);figure(1)subplot(1,2,1);hist(signal);title(叠加均匀分布随机数直方图);subplot(1,2,2);plot(n1,k1);title(叠加均匀分布的概率密度);结果2程序3clear all;clc; t=0:320; x=sin(2*pi*t*25); x1=wgn(1,321,0); z=x+x1; y=trapz(t,z); %y=int(z,x,0,t); subplot(3,2,1),plot(z); title(随机信号序列) meany=mean(z); figure(1)subplot(3,2,3),plot(t,meany,.); title(随机信号均值) vary=var(y); %方差 subplot(3,2,4),plot(t,vary,.); title(随机信号方差) cory=xcorr(z,unbiased);%自相关函数 subplot(3,2,2),plot(cory); title(随机信号自相关函数) covv=cov(y); subplot(3,2,5),plot(t,covv,.);title(随机信号协方差)t=[0:0.0005:0.045];X1=sin(2*pi*25*t);%正弦 figure(2)subplot(3,4,1);plot(t,X1);grid title(正弦函数序列); X2=randn(1,length(t)); %%X2=normrnd(2,0.04); %subplot(3,4,2); plot(t,X2); title(高斯噪声序列); X=X1+X2; %混合随机信号X(t)subplot(3,4,3); plot(t,X);gridtitle(混合随机信号); meany1=mean(X1); %subplot(3,4,6),plot(t,meany1);title(原信号均值); vary1=var(X1); %subplot(3,4,7),plot(t,vary1);title(原信号方差); cory1=xcorr(X1,unbiased); %subplot(3,4,8),plot(cory1);title(原信号自相关函数); meany=mean(X); %subplot(3,4,10),plot(t,meany);title(混合信号均值); vary=var(X); %subplot(3,4,11),plot(t,vary); title(混合信号方差)