%?---?減法濾波器matlab輔助設計文件?---
%?問題：濾波后相位發生了變化

clc;clear;close?all;

%%?----?根據傳遞函數繪出幅頻特性?---?%%
%?設置減法濾波器的傳遞函數系數
a1=1;?b1=[1?0?0?0?0?0?-1];?%傳遞函數?H（z）=1-?z^（-6）
f=0:1:600;?%頻率范圍
h1=abs（freqz（b1a1f1200））;%返回在向量f范圍內，采樣頻率為1200Hz下的幅頻特性

%?由傳遞函數系數確定傳遞函數的幅頻特性
H1=h1/max（h1）;?%?單位化，其實h1即是幅頻特性

%?繪出幅頻特性
figure;
plot（fH1）;
xlabel（‘f/Hz‘）;?ylabel（‘H1‘）;

%%?----?濾波仿真效果?----?%%
%?模擬輸入參數
N=24;
t1=（0:0.02/N:0.04）;?%兩個基頻周期
m=size（t1）;?%?m=[149]?1行49列

%?基波電壓
Va=100*sin（2*pi*50*t1）;

%疊加直流分量和48次諧波分量
Va1=35+100*sin（2*pi*50*t1）+30*sin（2*pi*200*t1）+10*sin（2*pi*400*t1）;

%采用減法濾波器濾掉Va1的直流分量和48次諧波分量
Y=zeros（16）;
for?jj=7:m（2）?%7:49?濾波器為6階，故從7開始?
????Y（jj）=（Va1（jj）-Va1（jj-6））/1.414;?%?為什么要除以sqrt（2）
????%Y（jj）=（Va1（jj）-Va1（jj-6））;
end

%?輸出波形
figure;
plot（t1Va‘-ro‘?t1Va1‘-bs‘?t1Y‘-g*‘）;
xlabel（‘t/s‘）;?ylabel（‘v/V‘）
grid?on;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????117??2014-08-18?21:13??book1.txt

?????文件??????35804??2014-08-14?16:31??example5_1_three_fault.mdl

?????文件??????42416??2014-08-18?21:23??example5_2_Generator_fault.mdl

?????文件??????63102??2014-08-14?19:45??example5_3_1_Ground_fault.mdl

?????文件??????65099??2014-08-14?19:50??example5_3_2_Ground_fault.mdl

?????文件???????1062??2014-08-15?11:03??example7_1_minus_filter.m

?????文件????????468??2014-10-18?13:13??example7_2_frontier_amplitude.m

?????文件???????1386??2014-08-18?08:53??example7_2_frontier_method.m

?????文件???????1247??2014-08-15?15:32??example7_2_twopoints_method.m

?????文件??????80175??2014-08-18?15:30??example7_3_2_direction_relay.mdl

?????文件??????74922??2014-08-18?15:30??example7_3_direction_relay.mdl

?????文件???????1047??2014-08-18?10:35??example7_3_frontier_method.m

?????文件??????50933??2014-08-18?17:00??example7_4_1_transformer.mdl

?????文件??????55165??2014-08-18?17:07??example7_4_2_transformer.mdl

?????文件??????46048??2014-08-18?21:23??example7_4_transformer.mdl

?????文件??????51568??2014-08-18?20:43??example7_5_travelling_wave.mdl

?????文件??????45653??2014-08-18?22:01??example9_1_wind_turbine.mdl

?????文件??????45495??2014-08-18?22:08??example9_2_wind_turbine_fault.mdl

?????文件?????528496??2014-08-19?09:25??example9_3_two_wind_turbine.mdl

?????文件???23332292??2015-05-16?21:19??MATLAB_Simulink電力系統建模與仿真?于群.pdf

?????文件?????358798??2014-05-19?21:29??MATPOWER手冊（中文版）.pdf

?????文件????????968??2014-08-18?21:17??xingbo.m

?????文件?????594208??2014-08-18?19:45??xingbo.mat

?????文件???????3683??2010-03-10?13:08??matpower4.0b2\@opf_model\add_constraints.m

?????文件???????4625??2010-03-10?13:08??matpower4.0b2\@opf_model\add_costs.m

?????文件???????2007??2010-03-10?13:08??matpower4.0b2\@opf_model\add_vars.m

?????文件???????3296??2010-03-10?13:08??matpower4.0b2\@opf_model\build_cost_params.m

?????文件???????2991??2010-03-10?13:08??matpower4.0b2\@opf_model\compute_cost.m

?????文件???????2890??2010-03-10?13:08??matpower4.0b2\@opf_model\display.m

?????文件????????583??2010-03-10?13:08??matpower4.0b2\@opf_model\get.m

............此處省略230個文件信息