离散系统Matlab信号处理

一、离散时间信号

代码：

n=-2:7
x=[0 2 3 5 6 -1 -5 7 9 -2]

subplot(2,1,1);
stem(n,x);
xlabel('n');ylabel('x(n)');title("stem函数绘制离散信号")

subplot(2,1,2);
plot(n,x);
hold on
plot(n,zeros(1,length(n)));
xlabel('时间/s');ylabel('函数值');title("plot函数绘制离散信号");
axis tight;

二、采样定理

代码：

clear all;

f1=10;dt=0.01;
%原始信号100个点，dt表示每个点的间隔
n1=0:100;t1=n1*dt;

%10Hz的正弦信号
x1=sin(2*pi*f1*t1+0.5);

subplot(5,1,1);
plot(t1,x1);title("10Hz原始信号")

subplot(5,1,2)
%原始信号
plot(t1,x1);hold on;

%采样频率20HZ，ts表示采样的每个点的间隔
fs=20;ts=1/fs;
%每秒钟采样20个点
n2=0:20;t2=n2*ts;
%得到的采样点
x2=sin(2*pi*f1*t2+0.5);
plot(t2,x2,'ro');title("采样")
%采样得到的信 号
subplot(5,1,3);plot(t2,x2);
ylim([-1 1]);title('20Hz采样后的信号');xlabel("时间/s");

subplot(5,1,4)
%原始信号
plot(t1,x1);hold on;
%采样频率20HZ，ts表示采样的每个点的间隔
fs=40;ts=1/fs;
%每秒钟采样20个点
n3=0:40;t3=n3*ts;
%得到的采样点
x3=sin(2*pi*f1*t3+0.5);
plot(t3,x3,'ro');title("采样")
%采样得到的信 号
subplot(5,1,5);plot(t3,x3);
ylim([-1 1]);title('40Hz采样后的信号');xlabel("时间/s");

 三、单位脉冲信号

代码：

clear all
n=0:20;
x=zeros(1,length(n));x(1)=1;
subplot(1,2,1);stem(n,x);

grid on;axis([-1 21 0 1.1]);
xlabel('n');title('\delta(n)');

%移位序列
x=zeros(1,length(n));x(7)=1;
subplot(1,2,2);stem(n,x);
grid on;axis([-1 21 0 1.1]);
xlabel('n');title('\delta(n-6)');

 四、单位阶跃信号

代码：

clear all
n=0:20;
x=ones(1,length(n));
subplot(1,2,1);stem(n,x);

grid on;axis([-1 21 0 1.1]);
xlabel('n');title('u(n)');

%移位序列
x=[zeros(1,8),ones(1,length(n))];
n1=0:length(n)+7

subplot(1,2,2);stem(n1,x);
grid on;axis([-1 29 0 1.1]);
xlabel('n');title('\delta(n-8)');

 五、斜坡信号

 五、正余弦信号

 代码：

clear all
n=0:30;
x=3*sin(n*pi/5+pi/4);
subplot(2,1,1)
stem(n,x);
xlabel('n');ylabel('x(n)');title('正弦序列');
axis([0 31 -3.1 3.1])


x=2*cos(n*pi/3-pi/4);
subplot(2,1,2)
stem(n,x);
xlabel('n');ylabel('x(n)');title('余弦序列');
axis([0 31 -2.1 2.1])

六、指数序列

代码：

clear all
n=0:10;

x1=1.2.^n;
x2=(-1.2).^n;
x3=0.6.^n;
x4=(-0.6).^n;

subplot(2,2,1);stem(n,x1);
xlabel('n');ylabel("x1(n)");title("x1(n)=1.2^n");

subplot(2,2,2);stem(n,x2);
xlabel('n');ylabel("x2(n)");title("x2(n)=(-1.2)^n");

subplot(2,2,3);stem(n,x3);
xlabel('n');ylabel("x3(n)");title("x3(n)=0.6^n");

subplot(2,2,4);stem(n,x4);
xlabel('n');ylabel("x4(n)");title("x4(n)=(-0.6)^n");

 七、复指数序列

代码：

clear all
n=0:20;

x=5*exp((i*pi/3+1/4)*n)
subplot(2,2,1);
stem(n,real(x),'fill');
xlabel('n');title('实部');
grid on

subplot(2,2,2);
stem(n,imag(x),'fill');
xlabel('n');title('虚部');
grid on


subplot(2,2,3);
stem(n,abs(x),'fill');
xlabel('n');title('模');
grid on

subplot(2,2,4);
stem(n,angle(x),'fill');
xlabel('n');title('相角');
grid on

八、周期序列

九、序列的加法和乘法

代码：

clear all

n=0:7;
x1=[-3 5 6 7 8 0 1 9];
x2=[6 -2 -3 5 -1 7 2 -4];

y1=x1+x2;
y2=x1.*x2;

subplot(2,2,1);stem(n,x1);
xlabel('n');title('x1(n)');

subplot(2,2,2);stem(n,x2);
xlabel('n');title('x2(n)');
 

subplot(2,2,3);stem(n,y1);
xlabel('n');title('x1(n)+x2(n)');

subplot(2,2,4);stem(n,y2);
xlabel('n');title ('x1(n)*x2(n)')