用matlab进行随机数处理作图与频谱分析

您好， 是这样的丛胡梁：

一、倒谱渗运分析原理

倒谱分析 matlab中的函数为

1) 复倒做念频谱 y=cceps(x) 2) 实倒频谱 y=rceps(x)

二、编程 ：

1.fft变换 t=0:0.01:2y1=sin(2*pi*30*t) subplot(2,1,1) plot(t,y1)

y2=abs(fft(y1,512)) f=100*(0:256)/512 subplot(2,1,2) plot(f,y2(1:257)) XlabelYlabel

2.功率谱变换

Fs=2000 NFFT=1024 t=0:1/Fs:1x=sin(2*pi*800*t)window=boxcar(length(x))periodogram(x,window,NFFT,Fs) 3.倒谱变换 Fs=200t=0:1/Fs:2

x=sin(2*pi*50*t).*sin(2*pi*40*t)z=cceps(x)subplot(2,1,1)plot(t,x)subplot(2,1,2)plot(t,z)

题主虚举想用matlab给出模型绘制一个图像。应该这样来实现：

1、首先需要有一系列数据 x、y

2、自定义模局返型函数表达式

func=@(a,x)a(1)/(1+(a(1)/x0-1)*exp(-a(2)*t))

3、利用nlinfit函数拟合系数Xm和r

a= nlinfit(x,y,func,a0)  a(1)=Xm，a(2)=r

4、求拟合值

yi=func(a,xi)

5、绘制模型图像

plot(x,y,xi,yi)

如差腊碧x0=10，r=20，xm=e^20，则其图形为

T=1954:2005

t=1:length(T)

x0=10r=20xm=exp(1)^20

x=xm./(1+(xm/x0-1)*exp(-r*t))

t=1954:2005

plot(t,x,'k-')

xlabel('t'),ylabel('x（t）')

---