demo

demo1，使用FFT函数进行频谱分析

t = 0:.001:.25;x = sin(2*pi*50*t) + sin(2*pi*120*t);

首先创建一些数据。 考虑以1000 Hz采样的数据。 从形成我们的数据的时间轴开始，从t = 0到t = 0.25，步长为1毫秒。 然后形成包含50Hz和120Hz的正弦波的信号x。

y = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')

添加具有标准偏差2的一些随机噪声以产生噪声信号y。通过绘制它看看这个嘈杂的信号y。

Y = fft(y,251);

对y进行傅里叶变换

Pyy = Y.*conj(Y)/251;f = 1000/251*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('Frequency (Hz)')

使用复共轭（CONJ）计算功率谱密度，各种频率下的能量的测量。 形成前127个点的频率轴，并使用它绘制结果。 （其余点是对称的。）

plot(f(1:50),Pyy(1:50))title('Power spectral density')xlabel('Frequency (Hz)')

hz放大五十倍

demo2，如何使用2维中的三次样条插值数据、

n = 8;x = rand(n,1);y = rand(n,1);plot(x,y,'.')axis([0 1 0 1])
建立坐标系

t = 1:n;ts = 1:1/10:n;xs = spline(t,x,ts);ys = spline(t,y,ts);hold onplot(xs,ys,'r');hold off

使用两个样条插补点，通过以更细的间距评估它们。 用红线绘制插值曲线。

输入一个行矩阵

可以分行输入行矩阵

MATLAB提示出现错误

用函数zeros生成全零阵

用函数zeros生成全零阵

矩阵加减运算

两个矩阵的乘法运算

矩阵的数乘运算

向量的点积

向量的混合积

左除和右除

矩阵的乘方

矩阵的转置

对奇异矩阵求逆时matlab给出的警告信息

用初等变换的方法求逆矩阵

以有理格式输出结果