【happyz】数字信号处理MATLAB学习-代码整理01

stem(n,x)；n为横轴，x为纵轴的线性图形。

N=length(t)；计算时间向量t的个数并赋给变量N。

axis([x1,x2,y1,y2])；横坐标从x1—x2，纵坐标从y1—y2。

x=zeros(1,n)；产生n个0的序列。

y=ones(1,n)；产生n个1的序列。

%抽样
%function[x,n]=chouyang(n1,n2,n0)

n=[n1:n2];
x=[(n-n0)==0];%n=n0时，n-n0=0,==0,此时x=1
[x,n]=chouyang(0,6,4);stem(n,x)

%指数

n=[-10:10];x=0.8.^n;stem(n,x)
n=[-10:10];
alpha=-0.1+0.3*j;
x=exp(alpha*n);
realx=real(x);%实部
imagex=imag(x);%虚部
magx=abs(x);%幅度
phasex=angle(x);%相位
subplot(2,2,1); stem(n,realx)
subplot(2,2,2);  stem(n,imagex)
subplot(2,2,3);  stem(n,magx)
subplot(2,2,4);  stem(n,phasex)

%正弦余弦

n=[0:10];
x=3*cos(0.1*pi*n+pi/3);
%   sin
stem(n,x)

%随机

rand(1,N);   %长度N，随机
randn(1,N);   %长度N,高斯随机

%周期

x=[1,2,3];               %一个x(n)
xn=x'*ones(1,3);         %生成p列x(n),符号'表示转置,3x3=9(序列长度9
xn=xn(:)'  ;              % :将p列串接成长列序列，'转置
stem(xn)

%单位抽样

N=5;n=0:N-1;
x=zeros(1,N);
x(2)=1;
stem(n,x,'.');grid;

%单位阶跃 函数u.m

n1=-2;n2=8;k=0;
n=n1:n2;
nt=length(n);
nk=abs(k-n1)+1;
x=[zeros(1,nk-1),ones(1,nt-nk+1)];
stem(n,x,'filled');
title('单位阶跃序列');
xlabel('时间(n)');
ylabel('幅度x(n)');

%移位

n0=0:5;x=sin(n0);[y,n]=sigshift(x,5,n0);stem(n,y)

%find函数功能：寻找非零元素的索引号。

find((n>=min(n1))&(n<=max(n1)));%在符号关系运算条件的范围内寻找非零元素的索引号。

%序列相加

 n1=0:10;x1=sin(n1*0.25*pi);
 n2=5:15;x2=sin(n2*0.25*pi);
 [y,n]=sigadd(x1,n1,x2,n2);stem(n,y)

%相乘

n1=0:10;x1=sin(n1*0.25*pi);
n2=5:15;x2=sin(n2*0.25*pi);
 [y,n]=sigcheng(x1,n1,x2,n2);stem(n,y)

%z变换
%ztrans 求无限长序列的z变换。

X=ztrans(x);

%iztrans 求函数X(z)的z反变换。

x= iztrans (X);

%syms定义符号对象。

sysm t,i,x;%将变量t,i,x声明为符号变量。

%residuez有理多项式的部分分式展开。

[r,p,c]=residuez(b,a);%求解以系数向量b，a表示的系统函数的部分分式展开。

%zplane绘制零极点分布图。

zplane(z,p);%绘制由列向量z确定的零点、列向量p确定的极点构成的零极点分布图。
zplane(b,a);%绘制由行向量b,a构成的系统函数确定的零极点分布图。
[hz,hp,ht]=zplane(z,p);%获得三个句柄向量：hz为零点线句柄；hp为极点线句柄；ht为坐标轴、单位圆及文本对象的句柄。