【无标题】

系列文章目录

提示：这里可以添加系列文章的所有文章的目录，目录需要自己手动添加
例如：第一章 Python 机器学习入门之pandas的使用

---

提示：写完文章后，目录可以自动生成，如何生成可参考右边的帮助文档

---

符号变量

用sym和syms定义
syms定义单个变量，sym定义多个变量

求导函数diff ( )

diff（自变量 应变量 求导阶数） 多元函数是输入一个因变量表示求偏导

示例：pandas 是基于NumPy 的一种工具，该工具是为了解决数据分析任务而创建的。

代码用法

> zeros（）和ones（） 表示生产全零或全一的数组 ones(行数，列数)
zeros (a ,b ) a为零矩阵行数，b为列数
>get ( ) 和 set ( ) 查变量属性和更改

**lsim (b,a x,t ) 解微分方程 b,a 是传递函数的分子系数和分母系数。x是输入信号，t是输入信号下标 。 lsim函数表示b，a确定的系统输入x，t确定的函数 的输出并绘制输出图形 ，如果不想输出图形，输入y=lsim (b,a x,t ) ；
**

step（b,a c） 和impulse(b,a,c) 函数 求阶跃和输出响应 b,a 表示传递函数的分子系数和分母系数 ，c表示坐标最大值

abs ( ) 取绝对值
ezplot（ ） 画符号变量**
ezplot (a ,b ) a表示符号变量，b表示自变量的范围
subplot ( ) 画分块图
subplot(2 2 1)
subplot(2 2 2)
subplot(2,2,[3 4])
@（输入参数列表）运算表达式
这样产生的函数句柄变量不指向特定的函数, 而是一个函数表达式
axis（） 限定x y坐标的范围
axis （a b c d）a b 是x的最小值和最大值，c d为y

square ()
square函数用于产生周期性矩形脉冲信号
一般square调用形式为y=square(t,duty)
duty指占空比，在这里调用时，可以理解为，在一个周期里幅值为1的脉冲的占比，默认50%

而square函数的默认周期为2*pi

故想要去自己想要的周期则需要先让 2*pi.*t 然后除以你的周期，就可以得到自己想要的周期。

据此可以做出任何自己想要的门函数等。

sinc(t) 画出sa（t）函数

x = sawtooth(t,xmax)
t是时间阵列，也就是时间轴；xmax这个参数的含义是这个锯齿波的峰值位置位于哪里，没有这个参数的话，其实默认为1，此时，峰值位于最右侧；如果设置为0，则峰值在左侧；可想而知，如果为0.5，则峰值位于中间。
rectpuls 是 MATLAB 中的一个函数，用于创建一个矩形脉冲信号。该函数的语法如下：
y = rectpuls(t, width)
其中，t 是时间轴，width 是脉冲的宽度。y 是输出的矩形脉冲信号。
tripuls 是 MATLAB 中的一个函数，用于创建一个三角信号。该函数的语法如下：
y = tripuls(t, width)
其中，t 是时间轴，width 是脉冲的宽度。y 是输出的矩形脉冲信号。

**impulse(b,a)**    	@绘出由向量a和b所定义的连续系统的冲激响应的时域波形。
impulse(b,a ,t0)        	@ 绘出由向量a和b所定义的连续系统在0 ~ t0时间范围内冲激响应的时域波形。
impulse(b,a,t1:p:t2)      	@绘出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的时域波形。
y=impulse(b,a,t1:p:t2);   	@求出由向量a和b所定义的连续系统在t1 ~ t2时间范围内，并且以时间间隔p均匀取样的冲激响应的数值解，但不绘出其相应波形。

.step函数（求阶跃响应）
step(num,den)
num表示分子各阶系数，den表示分母各阶系数

作用：已知闭环传递函数，求其阶跃响应
step(num,den)会输出图像；ct2=step(num,den)输出数值。

clear
clc   
num=[1 2 3];分子系数
den=[4 5 6] ;分母系数
step(num,den)%输出图像

**impz（b,a,n）b为分子系数，a为分母系数 n为下表**  求已知系统函数的冲击相应图 离散系统

**y = filter(b,a,x)**

其中，b和a是滤波器的系数向量，x是输入信号的向量。函数根据差分方程y[n] = b[1]*x[n] + b[2]*x[n-1] + ... + b[nb+1]*x[n-nb] - a[2]*y[n-1] - ... - a[na+1]*y[n-na]来计算输出信号y。

impulse 连续冲击
step 连续阶跃
lism 连续零状态
impz 离散冲击
filter 离散滤波器

roots（） 作用为求多项式的根
roots（[a,b,c]) a,b,c为方程由高到低次幂的系数。

freqs(B, A, w);
freqs 是 MATre 中频率值的个数。

freqz(B, A, w);
freqs 是 MATLAB 中的一个函数，用于计算给定连续时间复杂信号的频率响应。它的调用形式为：freqs(B, A, w)
其中，B 和 A 分别是系统的分子多项式和分母多项式系数向量，w 是一组频率值。函数返回一个包含频率响应的复数向量 H，其长度等于 w 中频率值的个数。

real() //求实部
imag() //求虚部
angle() // 求相位
abs() //求幅度

rand(n,m) //生产随机序列 nm
randn(n,m) //生产高斯随机序列 nm

ztrans(x); 求 式x的z变换
iztrans(X); 反z变换
laplace函数和ilaplace函数（拉氏逆变换）

Laplace
iLaplace
.tf函数（表示传递函数）
y=tf(num,den)
num表示分子各阶系数，den表示分母各阶系数
residue
[r,p,e]=residuez(b,a) b为分子系数，a为分母系数
r为幂级数系数，p为幂级数底数，e为常系数 由系统函数型转为并联型

zplane(b,a) b为分子系数，a为分母系数 求已知系统函数的零极点图
z zero零点
p polar 极点
lane 图

[b,a]=zp2tf(b,a ,k) 有零点序列和极点序列求系统函数分子和分母系数

conj函数：用于计算复数的共轭值

*dftmtx() 求序列x（n）的DFT　　
　dft（length（））x（n）‘
反dft xn=conj (dftmtx(4))/4*xk

tic ;toc;计算时间过了多少

nextpow2(N) ;求处大于等于N的最小二的系数。
pow2（N）； 求2的N次方。

巴特沃斯滤波器

wp=30;ws=35;Fs=100;rp=0.5;rs=40;
[n,Wn]=buttord(wp/(Fs/2),ws/(Fs/2),rp,rs,'z');//数字滤波器
[num,den]=butter(n,Wn);
[H,w]=freqz(num,den);plot(w*Fs/(2*pi),abs(H));

一型切比雪夫滤波器（带通滤波器）

Wp=[100 200]/500;Ws=[100-50 200+50]/500;
Rp=3;Rs=40;
[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs); %璋冪敤 Cheb1ord 鍑芥暟
[b,a]=cheby1(N,Rp,Wn);
freqz(b,a,512,1000);
title('Chebyshev Type I Bandpass Fitler')
axis([0,500,-80,0]);

二型切比雪夫滤波器（低通/ 高通）

fs=1000;fp=150;fr=100;
wp=fp/(fs/2);wr=fr/(fs/2);
Ap=1;Ar=40;
[N,wn]=cheb2ord(wp,wr,Ap,Ar);
[b,a]=cheby2(N,Ar,wn,'low'/'high');
[H,w]=freqz(b,a);
plot(w/pi*(fs/2),abs(H));

【n,wn】buttord(wp,ws,rp,rs，‘s’) 输入归一化的带通截至频率和带阻截至频率和带通衰减与带阻衰减，输出滤波器的阶数n和截止频率wn ‘s’：表示模拟滤波器（省略s表示数字滤波器）

【num,den】=butter（n，wm）输入滤波器的阶数n和截止频率wn，输出系统函数分子系数num，分子系数den；

【z,p,k】=buttap（n） 输入滤波器的系数 n，输出滤波器的零z,极点p，增益 k。
【b，a】=zp2tf(z,p,k) ,输入零极点和增益，输出系统函数的分子系数和分母系数。

boxcar：矩形窗 w = boxcar(n)
产生一个长度为n的矩形窗函数。

triang：三角窗。 w = triang(n)
triang(n)函数可产生一个长度为n的矩形窗函数。

bartlett：巴特利特窗。 w = bartlett(n)
bartlett(n)可得到n点的bartlett窗。bartlett窗函数系数为：

hamming：汉明窗。 w = hamming(n)
hamming(n)可产生n点的hamming窗。

hanning：汉宁窗。
w = hanning(n)
hanning(n)可产生n点的hanning窗。hanning窗的系数为：

blackman：布莱克曼窗。
w = blackman(n)
blackman(n)可产生n点的blackman窗。与等长度的hamming和hanning窗相比，blackman窗的主瓣稍宽，旁瓣稍低。

chebwin：切比雪夫窗。w = chebwin(n,r)
产生n点的窗函数，其傅里叶变换后的旁瓣波纹低于主瓣r dB。注意：当n为偶数时，窗函数的长度为(n+1)。

kaiser：凯瑟窗。w = kaiser(n,beta)

产生n点的kaiser窗函数，其中beta为影响窗函数旁瓣的函数，其最小的旁瓣抑制和之间的关系为：

printsys () 打印数据 打印传递函数

NUM = [1 0 2];
DEN = [2 3 1];
 
printsys(NUM,DEN,'s')
% printsys(NUM,DEN)    % 省略最后一个参数也是可以的

实现序列的周期延拓（mod）

DFT变换

DFT
xn% 信号
xk=dftmtx(length(xn))*xn;
Xn=conj(dftmtx(length(xn)))/length(xn)*xk

FFT K点序列的计算

求信号的功率 信号和信号的共轭相乘/ 信号总点数

1、传递函数转换为零极点增益（tf2zp）
由零极点增益转换为传递函数（zp2tf）

[z,p,k]=tf2zp(b,a)
[b,a]=zp2tf(z,p,k)
%a,b的长度要相等，不等的话要补零

2、由零极点增益转换为二次分式（zp2sos）
由二次分式转换为零极点增益（sos2zp）

[sos,g]=zp2sos(z,p,k),g为整个系统的增益，即H(z)=g*H1(z)*H2(z)…*HN(z)
sos2zp(z,p,k) =[sos,g],g为整个系统的增益,默认为1。

3、二次分式相乘转换为传递函数(sos2tf)
由传递函数转换为二次分式相乘 (tf2sos)

[b,a]=sos2tf(sos)
[sos,g]=tf2sos(b,a)

4.并联转直接型 residuez （也可以直接型转并联）

lsim函数：lsim函数是针对线性时不变模型，给定任意输入，得到任意输出。lsim函数表示任意输入函数的响应，连续系统对任意输入函数的响应可以利用lsim函数求取。

lsim(b,a,x,t)  @绘出由向量a和b所定义的连续系统在输入为x和t所定义的信号
y=lsim(b,a,x,t)  @不绘制出系统的**零状态响应**曲线，而只是求出与向量t定义的时间范围相一致的系统零状态响应的数值解。
lsim(A,B,C,D,e,t,X0)	@绘出由系数矩阵A,B,C,D所定义的连续时间系统在输入为e和t所定义的信号X0为初始状态
[Y,X]= lsim(A,B,C,D,e,t,X0) 	@不绘出全响应波形，而只是求出与向量t定义的时间范围相一致的系统输出向量Y的全响应以及状态变量X的数值解。

gensig
[u,t] = gensig(type,tau)
[u,t] = gensig(type,tau,Tf)
[u,t] = gensig(type,tau,Tf,Ts)
        u为信号序列，t为时间序列

        type为类型，包括：sin(正弦波)，square(方波)，pulse(周期脉冲)

        ta为type类型的周期，tf为持续时间，ts为采样时间。

H = heaviside(x)；并且heaviside是一个不连续的函数。它返回 0 表示 x < 0，返回 1/2 表示 x = 0，返回 1 表示 x > 0。