实验四源码的.docVIP

实验四： 第1小题源码： % x1=input(请输入初值：x1=); % x2=input(请输入初值：x2=); % p=input(请输入最大迭代次数：p=); % w=input(请输入迭代误差w=); clc clear x1=3;x2=3;p=100;w=0.0001; s=0; a=1; for k=1:1:p s=s+1; j11=2*x1-5;j12=2*x2; j21=2*x1;j22=2*x2-3; J=-[j11 j12;j21 j22]; f11=x1^2+x2^2-5*x1+1; f21=x1^2+x2^2-3*x2-3; f=[f11;f21]; [n,m]=size(J); for i=1:n J(i,i)=1/J(i,i); for j=i+1:n J(i,j)=J(i,j)*J(i,i); end for k=i+1:n for j=i+1:n J(k,j)=J(k,j)-J(k,i)*J(i,j); end end end for i=1:n f(i)=f(i)*J(i,i); for j=i+1:n f(j)=f(j)-J(j,i)*f(i); end end for i=n-1:-1:1 for j=i+1:-1:2 f(j-1)=f(j-1)-J(j-1,i+1)*f(i+1); end end x1=x1+f(1,1) x2=x2+f(2,1) z1(a)=x1; z2(a)=x2; a=a+1; s %方程组迭代次数 if x1^2+x2^2-5*x1+1w x1^2+x2^2-3*x2-3w break; else continue; end end plot(z1) hold on plot(z2,r) 第2小题例19.1源码： % n=input(请输入节点数：n=); % nl=input(请输入支路数：nl=); % isb=input(请输入平衡母线节点号：isb=); % pr=input(请输入误差精度：pr=); % B1=input(请输入由支路参数形成的矩阵:B1=); % B2=input(请输入各节点参数形成的矩阵:B2=); % X=input(请输入由节点号及其对地阻抗形成的矩阵：X=); n=5; nl=5; isb=1; pr=0.00001; B1=[1 2 0.03i 0 1.05 0;2 3 0.08+0.3i 0.5i 1 0;2 4 0.1+0.35i 0 1 0;3 4 0.04+0.25i 0.5i 1 0;3 5 0.015i 0 1.05 1]; B2=[0 0 1.05 1.05 0 1;0 3.7+1.3i 1 0 0 2;0 2+1i 1 0 0 2;0 1.6+0.8i 1 0 0 2;5 0 1.05 1.05 0 3]; X=[1 0;2 0;3 0;4 0;5 0]; Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);O=zeros(1,n);S1=zeros(nl); for i=1:n if X(i,2)~=0; p=X(i,1); Y(p,p)=1./X(i,2); end end for i=1:nl if B1(i,6)==0; p=B1(i,1);q=B1(i,2); else p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4)./2; Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4)./2; end G=real(Y);B=imag(Y); for i=1:n e(i)=real(B2(i,3)); f(i)=imag(B2(i,3)); V(i)=B2(i,4); end for