如何用matlab编写一个解决线性代数系统的程序

第一步：打开芦激MATLAB：如下

2

第二部：计算如下的线性代数问题：

3

第三步：冲哗毁在MATLAB中输入如散备下式子：（注意输入正确的字符和式子）

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

---