Matlab 学习笔记

未分类

1

借助符号计算所求的极限与理论值一致；
用数值法近似计算的极限与理论不一致。

x=eps;
L1=(1-cos(2*x))/(x*sin(x)), 
L2=sin(x)/x,    
syms t
f1=(1-cos(2*t))/(t*sin(t));
f2=sin(t)/t;
Ls1=limit(f1,t,0)
Ls2=limit(f2,t,0)

2

• digits 显示当前环境下十进制符号数字的有效位数
  digits(n) 把十进制符号数字有效位数设定为n
  xs=vpa(x) 据表达式x得到digits指定精度下的符号数字xs
  xs=vpa(x,n) 据表达式x得到n位有效数字的符号数字xs

• sym(Num) 采用精准数值类数创建精准的符号数字（推荐格式!）
  sc=sym(Num) 采用精准数值类数创建精准的符号常数sc（推荐格式!）

3绘图符号

(x,y,'s')：平面绘线三元组，'s'缺省时为“蓝色细实线”。

表5.2-2 离散数据点形允许设置值

表5.2-3 连续线形允许设置值

表5.2-4 点线色彩允许设置值

• grid 是否画分格线的双向切换指令（使当前分格线状态翻转）
  grid on 画出分格线
  grid off 不画分格线（默认设置）

• box 坐标形式在封闭式和开启式之间切换指令
  box on 使当前坐标呈封闭形式（默认设置）
  box off 使当前坐标呈开启形式

4 鼠标获取数据点

[x,y]=ginput(n) 用鼠标从二维图形上获取n个点的数据坐标 (x,y)

5 脉冲 阶跃

• 单位阶跃函数 function dtstep

function Y = dtstep
Y=zeros(size(X))
Y(X>=0)=1
end
n=0:9
y=dtstep(n-2)-dtstep(n-7)
stem(n,y)

• 单位脉冲信号
  两个阶跃函数差

6离散时间信号

1时间变换

1 时间翻转Y[n]=X[-n]

2时间展缩 Y[n]=X[an]

3 时间平移Y[n]=X[n-n0]

2幅度变换

1幅度翻转 A为负数

3幅度缩放

2幅度平移 B

7 时域 频域

观察时域信号： 示波器 产生实时的电压——频率图
观察频域信号时候 频谱分析仪 增益电压与频率

1.

1. 傅里叶级数——————周期性信号
2. 傅里叶变换——————非周期性信号

8 傅里叶变换

1 N=128 的基FFT 算法计算矩形脉冲频谱

T2=0:0.1:1.9;
xs=rectpuls(T2-0.5);%creat discrete-time rectangular pulse
m=zeros(1,108);%to make N=128
xs=[xs,m];
XS=fft(xs,128);
figure(1),y=abs(XS);plot(y)
XKS=fftshift(XS);% center the frequency spectrum around 0
np=-64:63;
figure(2),stem(np,abs(XKS));

2 N=256 正弦信号的DFT

t=0:0.001:0.255;
y=cos(100*t);
Y=fft(y,256);
XKS=fftshift(Y);
np=-128:127;
figure(3),plot(np,abs(XKS));

DFT计算过程 有个 隐含加窗处理 Xn加一个矩形窗
会导致信号频谱扩散

3 汉明窗 汉宁窗

t=0:0.001:0.255;
y=cos(100*t);
xham=hamming(256)'.*y;
Y=fft(xham);
XKS=fftshift(Y);
np=0:length(Y)-1;
figure(3),stem(np,abs(XKS));

9 MATLAB 数据拟合

1. cftool