Scribd
Upload
39 views

Assignment 3 Anurag Satpathy 116EE0343: For 3 Iterations

The document presents the results of an iterative algorithm to solve a system of equations over multiple iterations. The algorithm is run first without using a relaxation parameter, reaching convergence in 12 iterations. It is then run a second time using a relaxation parameter of 0.2, reaching convergence faster in 7 iterations. Each iteration displays the updated values of variables v2, v3 and v4 as well as the value of parameter q2.

Uploaded by

Anurag Satpathy
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
39 views

Assignment 3 Anurag Satpathy 116EE0343: For 3 Iterations

The document presents the results of an iterative algorithm to solve a system of equations over multiple iterations. The algorithm is run first without using a relaxation parameter, reaching convergence in 12 iterations. It is then run a second time using a relaxation parameter of 0.2, reaching convergence faster in 7 iterations. Each iteration displays the updated values of variables v2, v3 and v4 as well as the value of parameter q2.

Uploaded by

Anurag Satpathy
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Assignment 3

Anurag Satpathy 116EE0343


For 3 iterations
v1=1.04; p2=0.5; v2=1.04; v21=1; v3=1; v4=1; p3=-1; q3=0.5; p4=0.3; q4=-0.1; q2=0;
y=[3-9i -2+6i -1+3i 0; -2+6i 3.667-11i -0.667+2i -1+3i; -1+3i -0.667+2i 3.667-11i -2+6i; 0 -1+3i
-2+6i 3-9i];
for i=1:3
display('Iteration');
disp(i);
a=abs(v2);
q2 = -imag(conj(v2)*((y(2,1)*v1)+(y(2,2)*v2)+(y(2,3)*v3)+(y(2,4)*v4)))
if q2>=0.2 && q2<=1
v2=(((p2-q2*1i)/conj(v2))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2);
b=real(v2); c=imag(v2); t=atan(c/b);
v2=a*exp(1i*t)
end
if q2<0.2
q2=0.2
v2=(((p2-q2*1i)/conj(v21))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2)
v21=v2;
end
if q2>1
q2=1
v2=(((p2-q2*1i)/conj(v21))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2)
v21=v2;
end
v3=(((p3-q3*1i)/conj(v3))-(y(3,1)*v1)-(y(3,2)*v2)-(y(3,4)*v4))/y(3,3)
v4=(((p4-q4*1i)/conj(v4))-(y(4,1)*v1)-(y(4,2)*v2)-(y(4,3)*v3))/y(4,4)
end

Iteration 1
q2 = 0.2080
v2 = 1.0395 + 0.0335i
v3 = 1.0317 - 0.0894i
v4 = 1.0343 - 0.0151i

Iteration 2
q2 = -0.1563
q2 = 0.2000
v2 = 1.0669 + 0.0151i
v3 = 1.0470 - 0.0984i
v4 = 1.0541 - 0.0284i

Iteration 3
q2 = 0.1049
q2 = 0.2000
v2 = 1.0727 + 0.0080i
v3 = 1.0580 - 0.1056i
v4 = 1.0638 - 0.0361i
Without Alpha
v1=1.04; p2=0.5; v2=1.04; v21=1; v3=1; v4=1; p3=-1; q3=0.5; p4=0.3; q4=-0.1; q2=0; e1=[0,0,0,0];
y=[3-9i -2+6i -1+3i 0; -2+6i 3.667-11i -0.667+2i -1+3i; -1+3i -0.667+2i 3.667-11i -2+6i; 0 -1+3i
-2+6i 3-9i];
for i=1:50
v2oo=v2; v3oo=v3; v4oo=v4; q2oo=q2;
display('Iteration');
disp(i);
a=abs(v2);
q2 = -imag(conj(v2)*((y(2,1)*v1)+(y(2,2)*v2)+(y(2,3)*v3)+(y(2,4)*v4)))
if q2>=0.2 && q2<=1
v2=(((p2-q2*1i)/conj(v2))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2);
b=real(v2); c=imag(v2); t=atan(c/b);
v2=a*exp(1i*t)
end
if q2<0.2
q2=0.2
v2=(((p2-q2*1i)/conj(v21))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2)
v21=v2;
end
if q2>1
q2=1
v2=(((p2-q2*1i)/conj(v21))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2)
v21=v2;
end
v3=(((p3-q3*1i)/conj(v3))-(y(3,1)*v1)-(y(3,2)*v2)-(y(3,4)*v4))/y(3,3)
v4=(((p4-q4*1i)/conj(v4))-(y(4,1)*v1)-(y(4,2)*v2)-(y(4,3)*v3))/y(4,4)
e1(1)=v2oo-v2; e1(2)=v3oo-v3; e1(3)=v4oo-v4; e1(4)=q2oo-q2;
e=abs(max(e1));
if e<0.0001
break
end
end

Iteration 1
q2 = 0.2080
v2 = 1.0395 + 0.0335i
v3 = 1.0317 - 0.0894i
v4 = 1.0343 - 0.0151i

Iteration 2
q2 = -0.1563
q2 = 0.2000
v2 = 1.0669 + 0.0151i
v3 = 1.0470 - 0.0984i
v4 = 1.0541 - 0.0284i

Iteration 3
q2 = 0.1049
q2 = 0.2000
v2 = 1.0727 + 0.0080i
v3 = 1.0580 - 0.1056i
v4 = 1.0638 - 0.0361i

Iteration 4
q2 = 0.1293
q2 = 0.2000
v2 = 1.0775 + 0.0042i
v3 = 1.0636 - 0.1096i
v4 = 1.0693 - 0.0403i

Iteration 5
q2 = 0.1619
q2 = 0.2000
v2 = 1.0800 + 0.0021i
v3 = 1.0668 - 0.1117i
v4 = 1.0724 - 0.0427i

Iteration 6
q2 = 0.1787
q2 = 0.2000
v2 = 1.0814 + 0.0009i
v3 = 1.0686 - 0.1129i
v4 = 1.0741 - 0.0439i

Iteration 7
q2 = 0.1881
q2 = 0.2000
v2 = 1.0822 + 0.0003i
v3 = 1.0695 - 0.1136i
v4 = 1.0750 - 0.0446i

Iteration 8
q2 = 0.1934
q2 = 0.2000
v2 = 1.0826 - 0.0001i
v3 = 1.0701 - 0.1140i
v4 = 1.0756 - 0.0450i

Iteration 9
q2 = 0.1963
q2 = 0.2000
v2 = 1.0829 - 0.0003i
v3 = 1.0704 - 0.1142i
v4 = 1.0759 - 0.0453i

Iteration 10
q2 = 0.1980
q2 = 0.2000
v2 = 1.0830 - 0.0004i
v3 = 1.0706 - 0.1143i
v4 = 1.0760 - 0.0454i

Iteration 11
q2 = 0.1989
q2 = 0.2000
v2 = 1.0831 - 0.0004i
v3 = 1.0707 - 0.1143i
v4 = 1.0761 - 0.0454i

Iteration 12
q2 = 0.1994
q2 = 0.2000
v2 = 1.0831 - 0.0005i
v3 = 1.0707 - 0.1144i
v4 = 1.0762 - 0.0455i

With Alpha=0.2
v1=1.04; p2=0.5; v2=1.04; v21=1; v3=1; v4=1; p3=-1; q3=0.5; p4=0.3; q4=-0.1; q2=0; e1=[0,0,0,0];
alpha=0.2;
y=[3-9i -2+6i -1+3i 0; -2+6i 3.667-11i -0.667+2i -1+3i; -1+3i -0.667+2i 3.667-11i -2+6i; 0 -1+3i
-2+6i 3-9i];
for i=1:50
v2oo=v2; v3oo=v3; v4oo=v4; q2oo=q2;
display('Iteration');
disp(i);
a=abs(v2);
q2 = -imag(conj(v2)*((y(2,1)*v1)+(y(2,2)*v2)+(y(2,3)*v3)+(y(2,4)*v4)))
if q2>=0.2 && q2<=1
v2=(((p2-q2*1i)/conj(v2))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2);
b=real(v2); c=imag(v2); t=atan(c/b);
v2=a*exp(1i*t);
v2=v2+(alpha*(v2-v2oo));
b=real(v2); c=imag(v2); t=atan(c/b);
v2=a*exp(1i*t)
end
if q2<0.2
q2=0.2
v2=(((p2-q2*1i)/conj(v21))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2);
v2=v2+(alpha*(v2-v2oo))
v21=v2;
end
if q2>1
q2=1
v2=(((p2-q2*1i)/conj(v21))-(y(2,1)*v1)-(y(2,3)*v3)-(y(2,4)*v4))/y(2,2);
v2=v2+(alpha*(v2-v2oo))
v21=v2;
end
v3=(((p3-q3*1i)/conj(v3))-(y(3,1)*v1)-(y(3,2)*v2)-(y(3,4)*v4))/y(3,3);
v3=v3+(alpha*(v3-v3oo))
v4=(((p4-q4*1i)/conj(v4))-(y(4,1)*v1)-(y(4,2)*v2)-(y(4,3)*v3))/y(4,4);
v4=v4+(alpha*(v4-v4oo))
e1(1)=v2oo-v2; e1(2)=v3oo-v3; e1(3)=v4oo-v4; e1(4)=q2oo-q2;
e=abs(max(e1));
if e<0.0001
break
end
end

Iteration 1
q2 = 0.2080
v2 = 1.0392 + 0.0402i
v3 = 1.0380 - 0.1058i
v4 = 1.0461 - 0.0285i

Iteration 2
q2 = -0.2505
q2 = 0.2000
v2 = 1.0777 + 0.0021i
v3 = 1.0571 - 0.1079i
v4 = 1.0686 - 0.0416i

Iteration 3
q2 = 0.1883
q2 = 0.2000
v2 = 1.0789 + 0.0020i
v3 = 1.0681 - 0.1140i
v4 = 1.0738 - 0.0447i

Iteration 4
q2 = 0.1540
q2 = 0.2000
v2 = 1.0827 - 0.0004i
v3 = 1.0696 - 0.1142i
v4 = 1.0756 - 0.0454i

Iteration 5
q2 = 0.1989
q2 = 0.2000
v2 = 1.0828 - 0.0004i
v3 =1.0705 - 0.1145i
v4 = 1.0760 - 0.0455i

Iteration 6
q2 = 0.1966
q2 = 0.2000
v2 = 1.0831 - 0.0005i
v3 =1.0707 - 0.1144i
v4 = 1.0762 - 0.0455i

Iteration 7
q2 = 0.1999
q2 = 0.2000
v2 = 1.0832 - 0.0005i
v3 = 1.0708 - 0.1144i
v4 = 1.0762 - 0.0455i

Published with MATLAB® R2014a

You might also like