y=abs(g);%求其幅频响应
subplot(2,1,1),
plot(w,x*180/pi)%画相频响应曲线
xlabel('w');
ylabel('角度');
grid;
subplot(2,1,2),
plot(w,y)%画幅频响应曲线
xlabel('w');
ylabel('abs(H)');
grid
compass(x,y)
软件截图如图9所示：
4
4.1
（1） 相频特性曲线分析：由图可知， 相频特性曲线是一条下降的曲线，说明相角随着频率的增加而减小但是其值总是大于0.当频率相当大时相角越来越趋近于0。
幅频特性曲线分析：由图可得， 的幅频曲线是一条上升的曲线，说明电压传输系数H的幅值随ω的增大而增大，当ω为0时H为0，当ω趋近于无穷大时H的幅值越来越趋近于0.6.
而MATLAB实现功能主要依靠编程，它是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。
y=abs(g);%求其幅频响应
subplot(2,1,1),
plot(w,x*180/pi)%画相频响应曲线
xlabel('w');
ylabel('角度');
grid;
subplot(2,1,2),
plot(w,y)%画幅频响应曲线
xlabel('w');
ylabel('abs(H)');
grid
[3].张学敏.MATLAB基础及应用.北京：中国电力出版社，2009
[4].张德丰.MATLAB建模与仿真实例精讲.北京：机械工业出版社.2010
[5].汪建.电路原理.北京：清华大学出版社.2007
本科生课程设计成绩评定表
姓名
性别
专业、班级
课程设计题目：
课程设计答辩或质疑记录：
1:根据自己选择的题目，说明建立数学模型的步骤
clear
U=40.*2^0.5.*(0.2j-0.2)./(3.2j+2.65);
theta=angle(U)*180/pi;
h=abs(U)
图5
（4）求 ：
= = ，又因为 =40 所以得到:
=40 *(0.2j-0.2)/(5.85J-0.55)
利用matlab编程如下：
clear
U=40.*2^0.5.*(0.2j-0.2)./(5.85j-0.55);
2：在仿真过程中，你遇到了哪些问题？你是如何解决的？
3：通过本次课程设计，如何理解MATLAB/SIMULINK在本专业领域的应用？
成绩评定依据：
序号
评定项目
3.2MATLAB
3.2.1
（1）程序设计框图：
（2）获取 函数的频率响应的编程：
Clear
fz=[0.2 0.2 0];% 表达式分子的系数向量
fm=[0.35 3.23];% 表达式分母的系数向量
w=0:0.1:50;%设定频率的变化范围
g=freqs(fz,fm,w);%求其频率响应
x=angle(g);%求其相频响应
（5） 的频率响应曲线：
图2
3.2.2
（1）根据给出的各个元件参数先求出 ：
=1.5+
利用如下编程得到其化简后的值：
clear
z=1.5+(0.2.*i-0.2)./(1.9+2.*i);
theta=angle(z)*180/pi;
h=abs(z)
软件截图如图3：
图3
（2）求 ：
= =40 /（1.5026+0.1025j）
1.2
Simulink是MATLAB最重要的组件之一，他的实现功能是图像模拟，它提供一个动态系统建模、仿真和综合分析的集成环境。在该环境中，无需大量书写程序，而只需要通过简单直观的鼠标操作，就可构造出复杂的系统。Simulink具有适应面广、结构和流程清晰及仿真精细、贴近实际、效率高、灵活等优点，并基于以上优点Simulink已被广泛应用于控制理论和数字信号处理的复杂仿真和设计。同时有大量的第三方软件和硬件可应用于或被要求应用于Simulink。
关键词：电路原理，传递函数，MATLAB，向量图
基于MATLAB的电路频率响应分析
1
1.1
本次基础强化训练，我将利用MATLAB软件进行对频率响应电路进行原理分析、建模以及必要的推导和可行性分析。通过这一过程，掌握MATLAB软件的基本操作，体会MATLAB在实际中的应用价值，并且在以后的学习工作中利用MATLAB为自己服务，解决自己遇到的问题。
（2） 相频特性曲线分析：由图可知， 相频特性曲线是一条下降的曲线，说明 的相角随着频率的增加而减小，当频率趋近于无穷时，相角趋近于-90度。
幅频特性曲线分析:由图可知， 的幅频特性曲线是一条先上升后下降的曲线，当频率接近3.2的时候H的幅值最大。
4.2
利用matlab编程化简复数得到各个结果与人工笔算的结果一直，验证了编程的正确性，并且由上面的过程可见，用matlab来求解复数运算是相当简单的。而后将求解到的五个结果利用编程的方法调用compass函数来实现获取向量图，与正确的结果一致，由此可见利用matlab编程画向量图可以节省大量的时间和精力，是一种快捷方便的方法。
3.2.3
将3.2.3中的五个向量的实部与虚部分别构成两个向量x,y,然后调用compass函数画出其向量图。
获取向量图的编程如下：
clear;
clc;
x=[37.4727 0.3605 2.0973 -0.8684 38.34];
y=[-2.5562 3.8340 1.7368 1.0487 -3.6050];
clear
I=(0.3605+3.834j)./(0.1j);
theta=angle(I)*180/pi;
h=abs(I)
软件截图如图8所示：
图8
综合上述所得： =37.4727-2.5562j
=38.3400-3.6050j
=0.3605+3.8340j
=-0.8684+1.0487j
=2.0973+1.7368j
（3） 的频率响应曲线如图1所示：
图1
（4）获取 函数的频率响应的编程：
Clear
fz=[0.2 0.2 0];% 表达式分子的系数向量
fm=[量
w=0:0.1:50;%设定频率的变化范围
g=freqs(fz,fm,w);%求其频率响应
x=angle(g);%求其相频响应
2.3
第（1）小题：我们约定电感的支路为支路1，对应的阻抗记为 ，电容所在的支路记为支路2，其阻抗记为 ，总阻抗记为 。于是：
=jω*0.1，
=2+ ，
=1.5+ 。
然后列出传递函数： =
=
=
= = =
=
=
第（2）小题：已知电源的参数和各个元件的具体的数值，所以只要将各个元件的等效阻抗算出来，根据分压原理算得各个的电压值就可以把所求的各个量给解出来了。
利用matlab编程如下：
clear
I=40.*2^0.5/(1.5026+0.1025j);
theta=angle(I)*180/pi;
h=abs(I)
软件截图如图4所示：
图4
（3）求 ：
= ，又因为 =40 所以得到:
=40 *(0.2j-0.2)/(3.2j+2.65)
利用matlab编程如下：
总体上在这次基础强化训练中，我学到了关于Matlab的一些基础知识，但是这还远远不够，知识了解了一些基本的功能和操作，要想解决一些高级的问题我的知识还远远不够，我要在日后进一步学习，更好的掌握Matlab。
6
[1].邱关源.电路.北京：高等教育出版社，2006
[2].同济大学数学系.高等数学.北京：高等教育出版社，2007
3
3.1
（1）对于第一小题，我采用编程的方法来求得频率响应的曲线，具体的程序设计思路如下：首先输入H的关于频率的表达式，然后分别列出分子分母的系数行列式，调用angle函数来计算电路的相频特性，然后再调用abs函数来计算电路的幅频特性，调用freqs函数计算频率特性，最后用plot函数来描绘曲线。
（2）对于第二小题，由第一小题求出了 和 函数，根据其关系式，将ω=1带入可以得到 和 的向量值，接着可以求出所要求的各个值，所得到的表达式用matlab化简，最后得出结果。
2
2.1
2.2
（1）频率响应：电路中的感抗、容抗跟随激励源的频率变化，导致电路的工作状态亦随着频率变化，称此为频率特性，又称为频率响应。
（2）网络函数的定义：在线性正弦稳态网络中，当只有一个独立的激励作用时，网络中的某一响应（电压或电流）与网络输入之比，称为该响应的网络函数。
（3）本题中的电路中，只有一个激励源作用，并且电路属于线性稳态电路，电路中的电感 ,电容的等效阻抗为 ,这两个原件的工作状态会随着频率的变化而变化，进而改变电路的工作状态，通过列出对应的网络函数就可以容易的得出对应的幅频响应和相频响应。
5
通过Matlab多天的学习研究，我对Matlab有了初步的认识，我掌握了Matlab的基本操作，并且学会了用Matlab解决一些电路和数学上的问题。