资源描述：
《昆明理工大学-电力工程-学院学生实验报告.docx》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、昆明理工大学电力工程学院学生实验报告实验课程名称：《电力系统计算机辅助分析》开课实验室：计算机中心4082014年12月9日年级、专业、班电自114学号姓名成绩实验项目名称电力系统计算机辅助分析实验指导教师谢静教师签名：年月日 实验一 MATLAB软件的基本操作及程序的基本结构实验内容：1.验证实验2.对于如下电路图，已知V，，，，，自己编写程序，画出和的波形（四个周期）程序：clearmUsm=200;aUsm=(pi/180)*(+90);Usm=mUsm*exp(j*aUsm);R1=22,L=0.05,R2=33,C=0.0001;w=100*pi;Z1=R1+(w*0.
2、05*j);Z2=(1/(1/R2+j*w*C));Ucm=Z2/(Z1+Z2)*Usm;mUcm=abs(Ucm);aUcm=angle(Ucm);T=4*2*pi/w;N=100;dt=T/N;t=0:dt:T;fork=1:101us(k)=mUsm*cos(w*t(k)+aUsm);uc(k)=mUcm*cos(w*t(k)+aUcm);endplot(t,uc,'r--');holdonplot(t,us);title('电压uc，电压us');xlabel('时间/s','fontsize',16);ylabel('电压');legend('uc','us');gri
3、don波形思考题：1．matlab软件是什么？有什么特点？答：MATLAB是美国MathWorks公司自20世纪80年代中期推出的数学软件，具有优秀的数值计算能力和卓越的数据可视化能力，可以提供与矩阵有关的强大的数据处理和图形显示功能。特点有：1.编程效率高，2.用户使用方便，3.扩充能力强，交互性好，4.移植性和开放性很好，5.语句简单，内涵丰富，6.高效方便的矩阵和数组运算，7.方便的绘图功能。3.如何获得命令或函数的帮助文档？答：MATLAB的各个函数，不管是内建函数、M文件函数、还是MEX文件函数等，一般它们都有M文件的使用帮助和函数功能说明，各个工具箱通常情况下也具有一
4、个与工具箱名相间的M文件用来说明工具箱的构成内容等，在MATLAB命令窗口中，可以通过指令来获取这些纯文本的帮助信息。通常能够起到帮助作用，获取帮助信息的指令有help、lookfor、which、doc、get、type等。实验二 电力系统计算中常用的数值算法及电力网络的数学模型1.用因子表法求解电路方程中的节点电压；程序：clearclcclose%U1=10%-U1/10+U2*(1/10+1/30+1/50)-U3/50=0%-U1/15-U2/50+U3*(1/15+1/50+1/30)=0U1=2;A=[100;-1/101/10+1/30+1/50-1/50;-1/
5、15-1/501/15+1/50+1/30]B=[U1;0;0][n,m]=size(A);fori=1:nA(i,i)=1/A(i,i);forj=i+1:nA(i,j)=A(i,j)*A(i,i);endfork=i+1:nforj=i+1:nA(k,j)=A(k,j)-A(k,i)*A(i,j);endendenddisp('矩阵A的因子表为：');disp(A)fori=1:nB(i)=B(i)*A(i,i);forj=i+1:nB(j)=B(j)-A(j,i)*B(i);endendfori=n-1:-1:1forj=i+1:-1:2B(j-1)=B(j-1)-A(j-
6、1,i+1)*B(i+1);endenddisp('在因子表的基础上求解线性方程组的解为：x=');disp(B)2、将图2.1中的电压源分别改为2V、4V、6V、8V时，用因子表进行求解，体会该算法带来的益处；3上机验证课本例17.6－17.7的计算结果①例17.6clearclcclosen=5;nl=5;B=[120.03i01.050;230.08+0.3i0.5i10;340.015i01.051;250.1+0.35i010;350.04+0.25i0.5i10];X=[10;20;30;40;50];Y=zeros(n);fori=1:nifX(i,2)~=0;p=
7、X(i,1);Y(p,p)=1./X(i,2);endendfori=1:nlifB(i,6)==0p=B(i,1);q=B(i,2);elsep=B(i,2);q=B(i,1);endY(p,q)=Y(p,q)-1./(B(i,3)*B(i,5));Y(q,p)=Y(p,q);Y(q,q)=Y(q,q)+1./(B(i,3)*B(i,5)^2)+B(i,4)./2;Y(p,p)=Y(p,p)+1./B(i,3)+B(i,4)./2;enddisp('导纳矩阵Y=:');disp
