日志
并联谐振频率随R的变化曲线Matlab仿真
| |
天气: 舒适
心情: 高兴
参考书籍：RF circuit design：Theory and Applications 这个是书上的例题，程序自己写了，主要难点在于解复数方程，下面是源码，对于其他单谐振点的网络来说，稍微修改一下就可以使用，原图如下：
程序源码如下：
%并联谐振频率随电阻R变化关系图
%该程序只考虑了单个谐振点的情况，多个谐振点等遇到了再去考虑修改程序
clear all;
clc;
%给定元件值
L=10E-9;
C=1E-12;
R=1:200;
tic;
%解复数方程
syms w real %定义符号变量，real为了提取复数实部或虚部
for i=1:length(R)
Z(i)=1/(1/(j*w*L+R(i))+j*w*C);%并联阻抗表达式
a=(eval(solve(imag(Z(i)),'w')))'; %solve函数解出多个符号根；
%eval函数将其转变为数值根（方便参与运算）
%得到列向量，转置为行向量方便后续索引
index=find(a>0);%找出向量a中大于零的下标
%判断是否存在大于零的项，如果存在，则该值为角频率w；否则，令其为0
if index
root(i)=a(index);
else
root(i)=0;
end
end
toc;
omiga=root;
display(omiga)
f=omiga/(2*pi)/10e6 %MHz
plot(R,f,'r');
title('并联谐振频率随R变化的曲线图'),xlabel('R/Ohm'),ylabel('f/MHz');
运行结果：
/1 
eetop公众号
创芯大讲堂
创芯人才网