shixie2.m

用matlab来模拟振荡器谐振腔失谐情况

%无线送能模型的谐振分析
clear all;
L1 = 5e-6;      %发送端等效参数
C1 = 100e-12;    
Rs = 5;   
w1 = 1/sqrt(L1*C1);
Z11 = j*w1*L1-j/(w1*C1)+Rs;
Y11 = 1/Z11;
L2 = 5e-6;      %接收端等效参数
C2 = 100e-12;    
R2 = 100;            %按发光二级管正常发光时的等效电阻计算
w2 = 1/sqrt(L2*C2);
Z22 = j*w2*L2-j/(w2*C2)+R2;
Y22 = 1/Z22;

k = 0.1;            %线圈的耦合系数是0.1时
M = k*L1;
w0 = 1/sqrt(L1*C1);   %谐振点角频率
Zf2 = (w0*M)^2/Z11;
V0 = abs(j*w0*M/Z11*R2/(Z22+Zf2));

%失谐状态时，发送端由于受接收端影响，振荡频率回发生变化，假设变化是细小的，
%取振荡频率为：w0=sqrt(w1*w2)
q = 0.2:0.001:5;    %定义一个参量，表示接收端失谐量，q=w2/w1;
w = q*w1;
w0 = sqrt(w1.*w);
Z11 = j.*w0.*L1-j./(w0.*C1)+Rs;
Z22 = j.*w0.*L2-j./(w0.*C2)+R2;
V21 = j.*w0.*M./Z11;
Zf2 = (w0.*M).^2./Z11;
V = abs(V21.*R2./(Z22+Zf2))/V0;
plot(q,V,'r');
hold on;

%$$$$$$$$$$$$$$$$$$$$$$$$$$ k=0.2 $$$$$$$$$$$$$$$$$$$$$$$$$$
clear all;
L1 = 5e-6;      %发送端等效参数
C1 = 100e-12;    
Rs = 5;   
w1 = 1/sqrt(L1*C1);
Z11 = j*w1*L1-j/(w1*C1)+Rs;
Y11 = 1/Z11;
L2 = 5e-6;      %接收端等效参数
C2 = 100e-12;    
R2 = 100;            %按发光二级管正常发光时的等效电阻计算
w2 = 1/sqrt(L2*C2);
Z22 = j*w2*L2-j/(w2*C2)+R2;
Y22 = 1/Z22;

k = 0.2;            %线圈的耦合系数是0.2时
M = k*L1;
w0 = 1/sqrt(L1*C1);   %谐振点角频率
Zf2 = (w0*M)^2/Z11;
V0 = abs(j*w0*M/Z11*R2/(Z22+Zf2));

%失谐状态时，发送端由于受接收端影响，振荡频率回发生变化，假设变化是细小的，
%取振荡频率为：w0=sqrt(w1*w2)
q = 0.2:0.001:5;    %定义一个参量，表示接收端失谐量，q=w2/w1;
w = q*w1;
w0 = sqrt(w1.*w);
Z11 = j.*w0.*L1-j./(w0.*C1)+Rs;
Z22 = j.*w0.*L2-j./(w0.*C2)+R2;
V21 = j.*w0.*M./Z11;
Zf2 = (w0.*M).^2./Z11;
V = abs(V21.*R2./(Z22+Zf2))/V0;
plot(q,V,'g');

%$$$$$$$$$$$$$$$$$$$$$$$$$$ k=0.4 $$$$$$$$$$$$$$$$$$$$$$$$$$
%无线送能模型的谐振分析
clear all;
L1 = 5e-6;      %发送端等效参数
C1 = 100e-12;    
Rs = 5;   
w1 = 1/sqrt(L1*C1);
Z11 = j*w1*L1-j/(w1*C1)+Rs;
Y11 = 1/Z11;
L2 = 5e-6;      %接收端等效参数
C2 = 100e-12;    
R2 = 1000;            %按发光二级管正常发光时的等效电阻计算
w2 = 1/sqrt(L2*C2);
Z22 = j*w2*L2-j/(w2*C2)+R2;
Y22 = 1/Z22;

k = 0.4;            %线圈的耦合系数是0.1时
M = k*L1;
w0 = 1/sqrt(L1*C1);   %谐振点角频率
Zf2 = (w0*M)^2/Z11;
V0 = abs(j*w0*M/Z11*R2/(Z22+Zf2));

%失谐状态时，发送端由于受接收端影响，振荡频率回发生变化，假设变化是细小的，
%取振荡频率为：w0=sqrt(w1*w2)
q = 0.2:0.001:5;    %定义一个参量，表示接收端失谐量，q=w2/w1;
w = q*w1;
w0 = sqrt(w1.*w);
Z11 = j.*w0.*L1-j./(w0.*C1)+Rs;
Z22 = j.*w0.*L2-j./(w0.*C2)+R2;
V21 = j.*w0.*M./Z11;
Zf2 = (w0.*M).^2./Z11;
V = abs(V21.*R2./(Z22+Zf2))/V0;
plot(q,V,'k');
hold on;

%$$$$$$$$$$$$$$$$$$$$$$$$$$ k=0.5 $$$$$$$$$$$$$$$$$$$$$$$$$$
clear all;
L1 = 5e-6;      %发送端等效参数
C1 = 100e-12;    
Rs = 5;   
w1 = 1/sqrt(L1*C1);
Z11 = j*w1*L1-j/(w1*C1)+Rs;
Y11 = 1/Z11;
L2 = 5e-6;      %接收端等效参数
C2 = 100e-12;    
R2 = 100;            %按发光二级管正常发光时的等效电阻计算
w2 = 1/sqrt(L2*C2);
Z22 = j*w2*L2-j/(w2*C2)+R2;
Y22 = 1/Z22;

k = 0.5;            %线圈的耦合系数是0.1时
M = k*L1;
w0 = 1/sqrt(L1*C1);   %谐振点角频率
Zf2 = (w0*M)^2/Z11;
V0 = abs(j*w0*M/Z11*R2/(Z22+Zf2));

%失谐状态时，发送端由于受接收端影响，振荡频率回发生变化，假设变化是细小的，
%取振荡频率为：w0=sqrt(w1*w2)
q = 0.2:0.001:5;    %定义一个参量，表示接收端失谐量，q=w2/w1;
w = q*w1;
w0 = sqrt(w1.*w);
Z11 = j.*w0.*L1-j./(w0.*C1)+Rs;
Z22 = j.*w0.*L2-j./(w0.*C2)+R2;
V21 = j.*w0.*M./Z11;
Zf2 = (w0.*M).^2./Z11;
V = abs(V21.*R2./(Z22+Zf2))/V0;
plot(q,V,'b');