磁耦合谐振式充电系统RLC谐振仿真

% 分析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