%?分析RLC串聯諧振和并聯諧振特性
%?2018.12.11

clear
clc
close?all

R?=?10;
fn?=?6.78e6;?%?A4WP工作頻率6.78MHz
wn?=?2*pi*fn;
theta?=?0.1;??%?阻尼比越小能量衰減越慢，但起振速度也會越慢；對比分析0.01~0.1之間的阻尼比情況。

%%?RLC串聯諧振
L_s?=?R/（2*theta*wn）;
C_s?=?1/（wn^2*L_s）;
num_s?=?[1/L_s?0];
den_s?=?[1?R/L_s?1/（L_s*C_s）];
g_s?=?tf（num_sden_s）;
figurestep（g_s）

%%?RLC并聯諧振
C_p?=?1/（2*theta*wn*R）;
L_p?=?1/（wn^2*C_p）;
num_p?=?[1/（R*L_p*C_p）];
den_p?=?[1?1/（R*C_p）?1/（L_p*C_p）];
g_p?=?tf（num_pden_p）;
figurestep（g_p）

%%?諧振頻率激勵
%?模型離散化
M?=?200;
fs?=?M*fn;?%?采樣率設為諧振頻率的M倍
[n