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 on plot(xs,ys,'r'); hold off
使用两个样条插补点,通过以更细的间距评估它们。 用红线绘制插值曲线。
输入一个行矩阵
可以分行输入行矩阵
MATLAB提示出现错误
用函数zeros生成全零阵
用函数zeros生成全零阵
矩阵加减运算
两个矩阵的乘法运算
矩阵的数乘运算
向量的点积
向量的混合积
左除和右除
矩阵的乘方
矩阵的转置
对奇异矩阵求逆时matlab给出的警告信息
用初等变换的方法求逆矩阵
以有理格式输出结果