原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
CAAC实验报告-第次【DOC精选】
计算机辅助电路分析第1次上机报告
姓名：：：上机的心得体会及改进意见Matlab程序求解节点1、2、3、4的电压；
解：（1）混合方程矩阵形式（要求写出考虑不同支路情况时的步骤）
①不含无伴受控电压源和无伴受控电流源时的节点方程：
②接入无伴受控电压源时的节点方程：
③接入无伴受控电流源时的节点方程：
（2）matlab程序：
G1=0.2;G2=0.1;G3=0.5;w=10000;
C1=10e-6;C2=22e-6;L=0.1e-3;
u=0.5;g=1.5;
Is=10*exp(j*pi/9)
Us=100*exp(j*5*pi/18)
Y=[G1+j*w*C1 -j*w*C1 0 0 0;
-j*w*C1 j*w*(C1+C2)+1/(j*w*L) -j*w*C2 0 1;
-g -j*w*C2+g G2+j*w*C2 -G2 0;
0 0 -G2 G2+G3 0;
0 1 u -u 0]
I=[Us*G1 0 0 Is 0].
U=Y\I
Uabs=abs(U)
Uang=angle(U)*180/pi
（3）结果（写出时域表达式）
U =1.0e+002 *
0.5536 + 0.5926i
0.2068 + 0.7710i
-0.3083 - 1.7819i
0.1052 - 0.2400i
-0.1915 + 0.1281i
幅值：
Uabs =
81.0982
79.8204
180.8371
26.2038
23.0413
幅角：
Uang =
46.9488
74.9868
-99.8162
-66.3225
146.2121
各节点的电压时域表达式：
U1(t)=81.0982sin(10000t+46.9488°);
U2(t)= 79.8204sin(10000t+74.9868°);
U3(t)= 180.8371sin(10000t-99.8162°);
U4(t)= 26.2038sin(10000t-66.3225°);
Ik=23.0413sin(10000t+146.2121°)。
(4) 绘制出各节点电压及Ik的波形图
二、分别采用区间二分法和N-R法计算方程在[0，1]内的根，精确到小数点后4位，N-R法的迭代初值取x=0.2，要求用MATLAB求解。
要求：程序应正确完整，应为文本形式而非图片体现在本答案中，可以直接复制到Matlab命令窗口中运行。
解：
（1）区间二分法
求解函数：
function[x,yc]=bisection(fun_2,a,b,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
待求函数：
function y=fun_2(x)
y=exp(x)+2*x-4;
命令框函数：
[x,y]=bisection(fun_2,0,1,1e-3)
计算结果：
x =
0.8409
y =
4.3568e-004
（2）N-R法
评语及成绩评定
4
文档评论（0）