CAAC实验报告-第次参考答案【DOC精选】.doc

计算机辅助电路分析第1次上机报告 姓名：：：上机的心得体会及改进意见Matlab程序求解节点1、2、3、4的电压； 解：（1）混合方程矩阵形式 （2）matlab程序： G1=0.2,G2=0.1,G3=0.5,C1=10e-6,C2=22e-6,L=1e-4,miu=0.5,g=1.5; Is=10*exp(20/180*pi*j); Us=100*exp(50/180*pi*j); w=10000; Y=[G1+j*w*C1 -j*w*C1 0 0 0;-j*w*C1 j*w*C1+1/(j*w*L)+j*w*C2 -j*w*C2 0 1;-g -j*w*C2+g j*w*C2+G2 -G2 0;0 0 -G2 G2+G3 0;0 1 miu -miu 0]; I=[Us*G1;0;0;Is;0]; U=inv(Y)*I Uabs=abs(U) Uang=angle(U)*180/pi （3）结果（写出时域表达式） Uabs = 81.0982 79.8204 180.8371 26.2038 23.0413 Uang = 46.9488 74.9868 -99.8162 -66.3225 146.2121 时域U1=81.0982sin(10000t+46.9488); U2=79.8204sin(10000t+74.9868); U3=180.8371sin(10000t-99.8162); U4=26.2038sin(10000t-66.3225); Ik=23.0413sin(10000t+146.2121) 另解： 本题评点： 1）所谓混合节点方程，体现在待求变量向量中不只是原来的节点电压，还可能包括增加的支路电流，而在等式右侧的独立源向量中，可能同时包含有独立电流源和独立电压源，在导纳矩阵中的元素不仅有导纳，还有没有量纲的常量，或有关的转移导纳，电流比，电压比等。 2）有伴电压源的处理，解法一中节点1与参考节点间采用了有伴独立电压源，在列写节点1的节点电压方程时要注意转换为等效独立电流源，体现在等式右侧包含独立源。有伴电压源在改进节点分析中往往处理为无伴电压源，这时会增加一个独立节点，相应会增加一个支路电流，所以解法二中矩阵方程的规模比解法一中增大了2阶。解法二中，有伴电压源分离为无伴电压源支路和电阻支路，所以不需再对此有伴电压源进行诺顿模型转换。 3）注意区分受控源和独立源，只有独立源（包括独立电压源和独立电流源）才会影响等式右侧的电源向量。 4）注意参考方向对于相关系数的影响，如解法一中g前面的符号。 5）注意回顾电路原理相关知识中，电容和电感的阻抗和导纳的相量模型 6）注意本题中为正弦交流分析(jw)，不同于暂态分析（时域t或复频域s） 7）注意矩阵方程的matlab求解，建议使用inv函数，因为“/”“\”容易混淆，即使用错了也可能得到结果，不一定是求逆运算。 8）abs，angle等函数的基本用法 9）运算中将角度（degree,deg）转换为弧度（radian,rad），在后面的simulink中，由于软件指定了单位为deg时，要采用角度，或许是软件内部进行了转换。 10）一些特殊字符的输入，在matlab中不能输入时，用相近的字符代替，如表示角频率的ω可用w代替，也可拼音或英文表示，如omega，还有μ，可表达为miu，但不要用u表示，因为电压用u。 二、分别采用区间二分法和N-R法计算方程在[0，1]内的根，精确到小数点后4位，N-R法的迭代初值取x=0.2，要求用MATLAB求解。 解：（1）区间二分法 因f(0)=-1，f(1)=0.8415，故在[0,1]间存在实根。计算结果如表 MATLAB程序： function [x,yc]=bisection(fun_2,a,b,tol) % fun_2：f(x) % a b：求解区间 % x ：f(x)=0的根 % yc：f(c) % tol：精度 ya=feval(fun_2,a); yb=feval(fun_2,b); if ya*yb0,return,end while abs(ya-yb) tol x=(a+b)/2; yc=feval(fun_2,x); if ya*yc0 a=x;ya=yc; else b=x;yb=yc; end end 采用bisection求解例3-2-1。 先定义函数： function y=f(x) y=exp(x)+2*x-4; 执行MATLAB命令： [x,y]=bisection(f,0,1,1e-4) 由此，解得：x=0.8408，y= -2.5582e-0