急求用matlab画幅频特性和相频特性图。 请问怎么用matlab画这个函数的幅频特性曲线和相频特性曲线...
\u5982\u4f55\u7528matlab\u7ed8\u5236\u5e45\u9891\u548c\u76f8\u9891\u7279\u6027\u56fe1\u3001\u9996\u5148\u53cc\u51fbmatlab\u8f6f\u4ef6\u56fe\u6807\uff0c\u6253\u5f00matlab\u8f6f\u4ef6\uff0c\u53ef\u4ee5\u770b\u5230matlab\u8f6f\u4ef6\u7684\u754c\u9762\u3002
2\u3001\u4f7f\u7528\u76f4\u63a5\u8f93\u5165\u6cd5\u521b\u5efa\u4e00\u4e2a\u6570\u7ec4u \uff1au=[11,7,16,15,23]; \u8be5\u6570\u7ec4\u7528\u6765\u8868\u793a\u4e00\u7ec4\u6570\u636e\u3002
3\u3001\u4f7f\u7528\u76f4\u63a5\u8f93\u5165\u6cd5\u521b\u5efa\u4e00\u4e2a\u6570\u7ec4e \uff1ae=[2,1.3,1.5,3,1]; \u8be5\u6570\u7ec4\u7528\u6765\u8868\u793a\u6570\u7ec4u\u7684\u8bef\u5dee\u3002
4\u3001\u4f7f\u7528\u51fd\u6570errorbar(u,e);\u7ed8\u5236\u8bef\u5dee\u6761\u56fe\u3002
5\u3001\u4f7f\u7528\u51fd\u6570title()\u7ed9\u8bef\u5dee\u6761\u56fe\u6dfb\u52a0\u6807\u9898\u3002\u4f7f\u7528\u51fd\u6570xlabel('u')\u548cylabel('e'); \u7ed9\u5750\u6807\u8f74\u6dfb\u52a0\u6807\u7b7e\u3002
6\u3001\u6700\u540e\u67e5\u770b\u7ed8\u5236\u7684\u8bef\u5dee\u6761\u56fe\uff0c\u5c31\u5b8c\u6210\u4e86\u3002
1\u3001\u9996\u5148\u5b9a\u4e49bernstein\u57fa\u51fd\u6570\uff0c\u7528\u4e8e\u8ba1\u7b97\u5728\u7ed9\u5b9at\u70b9\u5217\u7684\u5bf9\u5e94\u7684bernstein\u57fa\u51fd\u6570\u7684\u503c\u3002
2\u3001\u7ed9\u51fa\u63a7\u5236\u9876\u70b9\u7684\u7684xy\u5750\u6807\u3002
3\u3001\u5b9a\u4e49\u53c2\u6570t\u7684\u70b9\u5217\uff0c\u5b9a\u4e49xx\u548cyy\u4e3a0\uff0c\u5206\u522b\u5b58\u50a8\u8ba1\u7b97\u5f97\u5230\u7684bezier\u66f2\u7ebf\u5bf9\u5e94\u7684x\u5750\u6807\u548cy\u5750\u6807\u7684\u503c\u3002
4\u3001\u8ba1\u7b97bezier\u66f2\u7ebf\u5728\u5bf9\u5e94\u70b9\u5217\u7684\u503c\u3002\u6839\u636ebezier\u66f2\u7ebf\u7684\u516c\u5f0f\u3002
5\u3001\u753b\u51fabezier\u66f2\u7ebf\u548c\u5bf9\u5e94\u7684\u63a7\u5236\u591a\u8fb9\u5f62\u3001\u63a7\u5236\u9876\u70b9\u3002
6\u3001\u7ed8\u5236\u51fa\u7684\u51fd\u6570\u56fe\u50cf\u5982\u56fe\u6240\u793a\u3002
clf;
Fs =128; %采样频率
T = 1/Fs; % 采样周期
N = 128; % 采样点数
t = (0:N-1)*T; % 时间，单位：S
x=zeros(N);
for n=0:N-1
b=fix((n)/(N/4));
m=n+1;
A=1/(N/4);
if b==0
x(m)=A*n;
elseif b==1||b==2
x(m)=A*(N/2-n);
elseif b==3
x(m)=A*(n-N);
end;
end;
n=0:N-1;
subplot(3,1,1)
plot(t,x);
xlabel('时间/S');
ylabel('振幅');title('时域波形');grid on;
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y)*2/N; %求取Fourier变换的振幅;*2/N转变为真实幅值
f=n*Fs/N;
subplot(3,1,2)
plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('幅值谱');grid on;
p=mod(angle(y)*180/pi,360);
subplot(3,1,3)
plot(f(1:N/2),p(1:N/2)); %绘出Nyquist频率之前随频率变化的相位
xlabel('频率/Hz');
ylabel('振幅');title('相位谱');grid on;
求取系统对数频率特性图（波德图）和系统奈奎斯特图（幅相曲线图或极坐标图）可以用函数bode()，格式如下：
bode(num,den);
bode(num,den,w);
[mag,phase,w]=bode(num,den,w);
%例如要画G（s）=10/(3*s^2+s+1)的伯德图，用以下程序：
num=[10]; %分子的常数
den=[3 2 1]; %分母的系数；
bode（num，den）
绛旓細1. 棣栧厛锛屽畾涔塀ernstein鍩哄嚱鏁帮紝鐢ㄤ簬璁＄畻缁欏畾t鍊兼椂瀵瑰簲鐨凚ernstein鍩哄嚱鏁扮殑鍊笺2. 缁欏嚭鎺у埗椤剁偣鐨剎鍜寉鍧愭爣銆3. 瀹氫箟涓绯诲垪鍙傛暟t鐨勫笺傚垵濮嬪寲xx鍜寉y鏁扮粍锛岀敤浜庡瓨鍌ㄨ绠楀緱鍒扮殑Bezier鏇茬嚎瀵瑰簲鐨剎鍜寉鍧愭爣鍊笺4. 璁＄畻Bezier鏇茬嚎鍦ㄧ粰瀹歵鍊肩偣鍒椾笂鐨勫搴旂偣銆5. 鏍规嵁Bezier鏇茬嚎鐨勫叕寮忥紝璁＄畻鏇茬嚎涓婃瘡涓偣鐨勫潗鏍...
绛旓細1銆佸湪matlab涓粯鍒跺箙棰戝拰鐩搁鐗规у浘鐨勫懡浠ゆ槸area锛屽懡浠ょ畝鍗曪紝鐢ㄦ硶澶氭牱銆2銆佸惎鍔╩atlab锛屽湪鍛戒护琛岀獥鍙ｈ緭鍏ヤ互涓嬪懡浠ゃ3銆佸緱鍒板搴旂殑鍥惧舰銆4銆佷笅闈㈠皢鍩哄间笅绉讳袱涓崟浣嶃5銆佸彲浠ョ湅鍒伴潰绉浘鐨勭旱鍧愭爣鎬讳綋涓嬬Щ2涓崟浣嶃6銆佷笅闈㈢殑鍛戒护灏嗛潰绉浘鐨勭嚎鍨嬬敱瀹炵嚎鏀逛负铏氱嚎銆7銆佸彲浠ョ湅鍒伴潰绉浘鐨勮竟缂樼嚎鍙樻垚浜嗚櫄绾裤8...
绛旓細5銆佺敾鍑篵ezier鏇茬嚎鍜屽搴旂殑鎺у埗澶氳竟褰佹帶鍒堕《鐐广6銆缁樺埗鍑虹殑鍑芥暟鍥惧儚濡傚浘鎵绀恒
绛旓細姹傚彇绯荤粺瀵规暟棰戠巼鐗规鍥撅紙娉㈠痉鍥撅級鍜岀郴缁熷濂庢柉鐗瑰浘锛骞鐩告洸绾垮浘鎴栨瀬鍧愭爣鍥撅級鍙互鐢ㄥ嚱鏁癰ode()锛屾牸寮忓涓嬶細bode(num,den);bode(num,den,w);[mag,phase,w]=bode(num,den,w);渚嬪瑕佺敾G锛坰锛=10/(3*s^2+s+1)鐨勪集寰峰浘锛岀敤浠ヤ笅绋嬪簭锛歯um=[10]; %鍒嗗瓙鐨勫父鏁 den=[3 2 1]; %鍒嗘瘝鐨...
绛旓細鎴戠敤鐨勬槸matlab2011b锛屼緵鍙傝冿紝鎴戝氨褰撲綘杩欎釜骞呴鏄凡缁忚绠椾簡20*Log锛堬級鍚庣殑浜嗭紝鐩搁宸茬粡杞寲涓哄姬搴﹀埗锛歠igure(1)plot(w,l)figure(2)plot(w,phi)l,phi涓骞呴鍜鐩搁锛寃涓洪鐜囷紝浠ヤ笂瀹屼簨浼氱敓鎴愪袱寮犲浘锛屽彲浠ラ塭dit-銆媋xis銆傘傘傜劧鍚庢妸涓嬮潰x閫夐」鍗￠噷鐨刲inear鏀规垚log灏辫浜 ...
绛旓細骞呴鐗规у拰鐩搁鐗规т竴鑸弿杩扮殑鏄煇涓鏁板瓧淇″彿绯荤粺鐨勬ц川锛屼緥濡傛垜浠婊ゆ尝鍣ㄧ綉缁滅殑骞呴鐗规с傚彲浠ユ牴鎹郴缁熷嚱鏁板湪S骞抽潰鐨勯浂鏋佺偣鍒嗗竷缁樺埗骞呴鐗规у拰鐩搁鐗规ф洸绾裤傛ゼ涓荤殑闂搴旇鏈夐棶棰橈紝绂绘暎搴忓垪鍙互缁樺埗鐨勬槸骞呭害鍜岀浉浣嶇殑娉㈠舰銆MATLAB涓垎鍒敤 abs(x)鍜(180/pi)*angle(x)鏉ユ眰锛岀劧鍚庣敤stem鏉嗗浘缁樺埗銆傜粰浣...
绛旓細鐢ㄨ繖涓嚱鏁癧Pxx,f] = periodogram(x,window,f,fs) x鏄椂鍩熺殑鏁版嵁 杩欓噷window鏄獥鍑芥暟锛屾埅鍙栦綘瑕侀儴鍒嗭紝閮借鍙互杩欐牱鍐 L = length(x) window = rectwin(L); f鏄鍩熺殑鑼冨洿鏍规嵁浣犵殑闇瑕侀 f = 0:01:100; fs = 1/T; T鏄椂鍩熶俊鍙锋椂闂撮棿闅matlab鐢诲箙棰戠壒鎬鏇茬嚎 ...
绛旓細榛樿鎯呭喌涓婲=512銆傜敱浜庤皟鐢╢reqz()鍑芥暟鍙兘姹傚嚭绂绘暎绯荤粺棰戠巼鍝嶅簲鐨勬暟鍊硷紝涓嶈兘鐩存帴缁樺埗鏇茬嚎鍥撅紝鍥犳锛屽彲浠ュ厛鐢╢reqz()鍑芥暟姹傚嚭绯荤粺棰戠巼鍝嶅簲鐨勫硷紝鐒跺悗鍐鍒╃敤MATLAB鐨abs()鍜宎ngle()鍑芥暟浠ュ強plot()鍛戒护锛屽嵆鍙粯鍒跺嚭绯荤粺鍦ㄦ垨鑼冨洿鍐呯殑骞呴鐗规у拰鐩搁鐗规ф洸绾裤
绛旓細璋冪敤bode鍑芥暟灏卞彲浠ュ緱鍒 渚嬪锛>> s=tf('s');>> G=(s+8)/(s*(s^2+0.2*s+4)*(s+1)*(s+3));>> bode(G)灏卞彲浠ュ緱鍒拌繖涓や釜鍥
绛旓細鐢ㄨ繖涓嚱鏁癧Pxx,f]= periodogram(x,window,f,fs)x鏄椂鍩熺殑鏁版嵁 杩欓噷window鏄獥鍑芥暟锛屾埅鍙栦綘瑕侀儴鍒嗭紝閮借鍙互杩欐牱鍐 L = length(x)window = rectwin(L);f鏄鍩熺殑鑼冨洿鏍规嵁浣犵殑闇瑕侀 f = 0:01:100;fs = 1/T;T鏄椂鍩熶俊鍙锋椂闂撮棿闅matlab鐢诲箙棰戠壒鎬鏇茬嚎 ...
