求 matlab 中bode源代码 matlab中分数阶线性系统的bode图程序看不懂，求高手指...
matlab \u7528bode()\u547d\u4ee4\u7ed8\u5236\u7684bode\u56fe \u5982\u4f55\u63d0\u53d6bode\u56fe \u6570\u636e\u554a\uff1f\u5728matlab \u63d0\u53d6bode\u56fe\u6570\u636e\u5982\u4e0b\uff1a
1\u4e36\u6211\u4eec\u73b0\u5728matlab\u91cc\u9762\u7f16\u5199\u597d\u81ea\u5df1\u8981\u7ed8\u5236\u56fe\u5f62\u7684\u4ee3\u7801\u3002
2\u4e36\u5728\u4ee3\u7801\u540e\u9762\u52a0\u4e0a[mag,phase,w] = bode(sys);
3\u4e36\u8fd9\u65f6\u5019\u6211\u4eec\u5728workplace\u754c\u9762\u5c31\u53ef\u4ee5\u770b\u5230\u6211\u4eec\u9700\u8981\u7684\u6570\u636e\u4e86\u3002
\u8fd9\u6837\u5c31\u89e3\u51b3\u4e86\u5728matlab \u63d0\u53d6bode\u56fe\u6570\u636e\u7684\u95ee\u9898\u3002
function H=f_bode(a,eta,b,gamma,w,Gc)
%\u5206\u6570\u9636\u7ebf\u6027\u7cfb\u7edf\u7684bode\u56fe\u7a0b\u5e8f
%a,b,eta,gamma\u5747\u4e3a\u6574\u6570\uff0cw\u5b9e\u503c\u77e2\u91cf\uff0cGc\u6743\u91cd\uff0c\u8f93\u51faH\u662ffrd\u6a21\u578b\uff0c\u4f8b\u5982\uff1a
%a=3;b=1;eta=3;gamma=2;freq = logspace(1,2);w = 0.05*(freq);f_bode(a,eta,b,gamma,w)
%a=3;b=1;eta=3;gamma=2;freq = 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); %\u751f\u6210\u9891\u7387\u54cd\u5e94\u6570\u636e\u6a21\u578b
if nargin==6,
H=H*Gc;
end
if nargout==0,
bode(H);
end
%BODE Bode frequency response of LTI models.
%
% BODE(SYS) draws the Bode plot of the LTI model SYS (created with
% either TF, ZPK, SS, or FRD). The frequency range and number of
% points are chosen automatically.
%
% BODE(SYS,{WMIN,WMAX}) draws the Bode plot for frequencies
% between WMIN and WMAX (in radians/second).
%
% BODE(SYS,W) uses the user-supplied vector W of frequencies, in
% radian/second, at which the Bode response is to be evaluated.
% See LOGSPACE to generate logarithmically spaced frequency vectors.
%
% BODE(SYS1,SYS2,...,W) graphs the Bode response of multiple LTI
% models SYS1,SYS2,... on a single plot. The frequency vector W
% is optional. You can specify a color, line style, and marker
% for each model, as in
% bode(sys1,'r',sys2,'y--',sys3,'gx').
%
% [MAG,PHASE] = BODE(SYS,W) and [MAG,PHASE,W] = BODE(SYS) return the
% response magnitudes and phases in degrees (along with the frequency
% vector W if unspecified). No plot is drawn on the screen.
% If SYS has NY outputs and NU inputs, MAG and PHASE are arrays of
% size [NY NU LENGTH(W)] where MAG(:,:,k) and PHASE(:,:,k) determine
% the response at the frequency W(k). To get the magnitudes in dB,
% type MAGDB = 20*log10(MAG).
%
% For discrete-time models with sample time Ts, BODE uses the
% transformation Z = exp(j*W*Ts) to map the unit circle to the
% real frequency axis. The frequency response is only plotted
% for frequencies smaller than the Nyquist frequency pi/Ts, and
% the default value 1 (second) is assumed when Ts is unspecified.
%
% See also BODEMAG, NICHOLS, NYQUIST, SIGMA, FREQRESP, LTIVIEW, LTIMODELS.
% Old help
%warning(['This calling syntax for ' mfilename ' will not be supported in the future.'])
%BODE Bode frequency response for continuous-time linear systems.
% BODE(A,B,C,D,IU) produces a Bode plot from the single input IU to
% all the outputs of the continuous state-space system (A,B,C,D).
% IU is an index into the inputs of the system and specifies which
% input to use for the Bode response. The frequency range and
% number of points are chosen automatically.
%
% BODE(NUM,DEN) produces the Bode plot for the polynomial transfer
% function G(s) = NUM(s)/DEN(s) where NUM and DEN contain the
% polynomial coefficients in descending powers of s.
%
% BODE(A,B,C,D,IU,W) or BODE(NUM,DEN,W) uses the user-supplied
% frequency vector W which must contain the frequencies, in
% radians/sec, at which the Bode response is to be evaluated. See
% LOGSPACE to generate logarithmically spaced frequency vectors.
% When invoked with left hand arguments,
% [MAG,PHASE,W] = BODE(A,B,C,D,...)
% [MAG,PHASE,W] = BODE(NUM,DEN,...)
% returns the frequency vector W and matrices MAG and PHASE (in
% degrees) with as many columns as outputs and length(W) rows. No
% plot is drawn on the screen.
%
% See also LOGSPACE, SEMILOGX, MARGIN, NICHOLS, and NYQUIST.
% J.N. Little 10-11-85
% Revised A.C.W.Grace 8-15-89, 2-4-91, 6-21-92
% Revised Clay M. Thompson 7-9-90
% Revised A.Potvin 10-1-94
% Copyright 1986-2003 The MathWorks, Inc.
% $Revision: 1.1.8.2 $ $Date: 2005/12/22 17:45:13 $
ni = nargin;
no = nargout;
% Check for demo and quick exit
if ni==0,
eval('exresp(''bode'')')
return
end
error(nargchk(2,6,ni));
% Determine which syntax is being used
switch ni
case 2
if size(a,1)>1,
% SIMO syntax
a = num2cell(a,2);
den = b;
b = cell(size(a,1),1);
b(:) = {den};
end
sys = tf(a,b);
w = [];
case 3
% Transfer function form with time vector
if size(a,1)>1,
% SIMO syntax
a = num2cell(a,2);
den = b;
b = cell(size(a,1),1);
b(:) = {den};
end
sys = tf(a,b);
w = c;
case 4
% State space system without iu or time vector
sys = ss(a,b,c,d);
w = [];
otherwise
% State space system, with iu but w/o time vector
if min(size(iu))>1,
error('IU must be a vector.');
elseif isempty(iu),
iu = 1:size(d,2);
end
sys = ss(a,b(:,iu),c,d(:,iu));
if ni<6,
w = [];
end
end
if no==0,
bode(sys,w)
else
[magout,phase,w] = bode(sys,w);
[Ny,Nu,lw] = size(magout);
magout = reshape(magout,[Ny*Nu lw]).';
phase = reshape(phase,[Ny*Nu lw]).';
end
% end bode
扩展阅读：matlab nextpow2 ... matlab while ... matlab ode45 ... matlab figure ... matlab shading interp ... matlab adobe rgb ... matlab function ... matlab save txt ... matlab vpa det ...
