simulink中分数阶模块如何搭建

打开软件。

1、打开Simulink软件，创建一个新模型。

2、选择LibraryBrowser（库浏览器），打开Simulink模型库。

3、在Simulink模型库中搜索分数阶模型。

4、在库浏览器中找到FractionalOrderDynamics，然后选取需要的分数阶模块，拖动到模型窗口中。

5、连接其他模块或信号源以完善模型，连接方式和常规Simulink模型建立方式一样。

6、对分数阶模块进行参数设置，包括分数阶指数、采样周期等。

7、进行仿真，可通过Simulink模型中的Scope模块来展示模型的输出结果。

function Faf = frft(f, a)

% The fast Fractional Fourier Transform

% input: f = samples of the signal

%a = fractional power

% output: Faf = fast Fractional Fourier transform

error(nargchk(2, 2, nargin))

f = f(:)

N = length(f)

shft = rem((0:N-1)+fix(N/2),N)+1

sN = sqrt(N)

a = mod(a,4)

% do special cases

if (a==0), Faf = freturnend

if (a==2), Faf = flipud(f)returnend

if (a==1), Faf(shft,1) = fft(f(shft))/sNreturnend

if (a==3), Faf(shft,1) = ifft(f(shft))*sNreturnend

% reduce to interval 0.5 <a <1.5

if (a>2.0), a = a-2f = flipud(f)end

if (a>1.5), a = a-1f(shft,1) = fft(f(shft))/sNend

if (a<0.5), a = a+1f(shft,1) = ifft(f(shft))*sNend

% the general case for 0.5 <a <1.5

alpha = a*pi/2

tana2 = tan(alpha/2)

sina = sin(alpha)

f = [zeros(N-1,1) interp(f) zeros(N-1,1)]

% chirp premultiplication

chrp = exp(-i*pi/N*tana2/4*(-2*N+2:2*N-2)'.^2)

f = chrp.*f

% chirp convolution

c = pi/N/sina/4

Faf = fconv(exp(i*c*(-(4*N-4):4*N-4)'.^2),f)

Faf = Faf(4*N-3:8*N-7)*sqrt(c/pi)

% chirp post multiplication

Faf = chrp.*Faf

% normalizing constant

Faf = exp(-i*(1-a)*pi/4)*Faf(N:2:end-N+1)

function xint=interp(x)

% sinc interpolation

N = length(x)

y = zeros(2*N-1,1)

y(1:2:2*N-1) = x

xint = fconv(y(1:2*N-1), sinc([-(2*N-3):(2*N-3)]'/2))

xint = xint(2*N-2:end-2*N+3)

function z = fconv(x,y)

% convolution by fft

N = length([x(:)y(:)])-1

P = 2^nextpow2(N)

z = ifft( fft(x,P) .* fft(y,P))

z = z(1:N)

function H=f_bode(a,eta,b,gamma,w,Gc)

%分数阶线性系统的bode图程序

%a,b,eta,gamma均为整数，w实值矢量，Gc权重，输出H是frd模型，例如：

%a=3b=1eta=3gamma=2freq = logspace(1,2)w = 0.05*(freq)f_bode(a,eta,b,gamma,w)

%a=3b=1eta=3gamma=2freq = linspace(0.01,0.2)w = 0.05*(freq)f_bode(a,eta,b,gamma,w)

H=zeros(size(w))

for i=1:length(w)

P=b*((sqrt(-1)*w(i)).^gamma.')

Q=a*((sqrt(-1)*w(i)).^eta.')

H(i)=P/Q

end

H=frd(H,w)%生成频率响应数据模型

if nargin==6,

H=H*Gc

end

if nargout==0,

bode(H)

end

---