在matlab中搭建分数阶变量

function Myexamp1 clear,clc, tspan = [0 20]% t 的取值范围 y0=[01]% 初始值 y(t0) 和 y'(t0) options = odeset('RelTol'磨蠢,1e-4,'AbsTol',[1e-4 1e-4])[t,y] = ode45(@fun,tspan,y0 ,options)plot(t,y(:,1),'-',t,y(:,2),'-.') title('y'瞎迟陪'旦饥(t)= -3 cos(2t) +2sin(t)+ t-3.8') function dy = fun(t,y) dy = zeros(2,1)dy(1) = y(2)dy(2) = -3-cos(2*t) + 2*sin(t)+t-3.8

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

---