数字信号处理利用matlab中的Fft分析信号

clc
clear
%2.(1)
X1=1:5
X2=1:0.500:2.5
X3=5:-1:1
%2.(2)
A1=zeros(3,3)
A2=zeros(3,4)
A3=ones(3,3)
A4=ones(3,4)
A5=ceil(rand(3,3)*10)
A6=ceil(rand(3,4)*20)
%2.(3)
A=[1 2 3
    4 5 6
    7 8 9];
B=[1.0 1.1 1.2
    2.0 2.1 2.2
    3.0 3.1 3.2];
B3=A+B
B4=A+B
B5=A-B
B6=A+3
B7=A-4
B8=A*B
B9=A.*B
B10=inv(A)
B11=A/B
B12=A./B
%2.(4)
x1=1:0.5:3;
C1=sin(x1)
C2=cos(x1)
C3=exp(x1)
C4=log(x1);
C5=sqrt(x1);
%3
x1=1:0.5:15;
stem(x1,sin(x1),'.')
stem(x1,cos(x1),'.')
stem(x1,exp(x1),'.')
stem(x1,log(x1),'.')
stem(x1,sqrt(x1),'.')
%4
n=0:10;
z1=-1/12+j*pi/6;
z2=1/12+j*pi/6;
z3=1/12;
z4=2+j*pi/6;
z5=j*pi/6;
stem(exp(real(z1)*n))
stem(exp(imag(z1)*n))
stem(exp(real(z2)*n))
stem(exp(imag(z2)*n))
stem(exp(real(z3)*n))
stem(exp(imag(z3)*n))
stem(exp(real(z4)*n))
stem(exp(imag(z4)*n))
stem(exp(real(z5)*n))
stem(exp(imag(z5)*n))
%5
n=0:1:30;
stem(n,1.5*sin(2*pi*0.1*n))
stem(n,1.5*sin(2*pi*0.9*n))
f1=100;
fs=1000;
N=100;
n=0:N-1;
t=n/fs;
plot(sin(2*pi*f1*t))