t=0:0.001:6; x=sin(2*pi*50*t)+sin(2*pi*120*t)+randn(1,length(t)); % A = randn(m,n) 或A=randn([m n]):返回一个m*n的随机项矩阵。 % size(A):获取数组A的行数和列数 %length(A):获取数组长度(即行数或列数中的较大值) Y=fft(x,512); P=Y.*conj(Y)/512; %'.*'运算符用于矩阵间对应元素的相乘,或数与数之间,数与矩阵之间的相乘。 %注:a*b表示矩阵a与矩阵b进行矩阵相乘。 % a.*b表示矩阵a中的元素与矩阵b中的元素按照相同位置进行相乘,得到的结果作为新矩阵中相同位置的元素。 % conj(A):返回A的共轭值 f=1000*(0:255)/512; plot(f,P(1:256))12345678910111213
2.进行信号检测。分析信号频谱对应的数字频率和模拟频率之间的关系。模拟信号x(t)=2*sin(4*pi*t)+5*cos(8*pi*t),以t=0.01n(0