实验数据.docx
《实验数据.docx》由会员分享,可在线阅读,更多相关《实验数据.docx(26页珍藏版)》请在冰豆网上搜索。
实验数据
>>clear
G(1,1)=0;B(1,1)=-5.;G(1,2)=0;B(1,2)=2.5;
G(1,3)=0;B(1,3)=2.5;G(2,1)=G(1,2);B(2,1)=B(1,2);
G(2,2)=0;B(2,2)=-2.5;G(2,3)=0;B(2,3)=0;
G(3,1)=G(1,3);B(3,1)=B(1,3);G(3,2)=G(2,3);B(3,2)=B(2,3);
G(3,3)=0;B(3,3)=-2.5;
Y=G+j*B
delt
(1)=0;delt
(2)=0;u
(1)=1;u
(2)=1;
p
(2)=-1;q
(2)=-0.6;p
(1)=0.2;q
(1)=0.1;
k=1;
precision=1;
N1=2;
whileprecision>0.4
delt(3)=0;u(3)=1;
form=1:
N1
forn=1:
N1+1
pt(n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
qt(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
end
pp(m)=p(m)-sum(pt);
qq(m)=q(m)-sum(qt);
end
form=1:
N1
forn=1:
N1+1
h0(n)=u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
n0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
j0(n)=-u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
L0(n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
end
H(m,m)=sum(h0)+u(m)^2*B(m,m);N(m,m)=sum(n0)-u(m)^2*G(m,m);
J(m,m)=sum(j0)+u(m)^2*G(m,m);L(m,m)=sum(L0)+u(m)^2*B(m,m);
end
form=1:
N1
JJ(2*m-1,2*m-1)=H(m,m);
JJ(2*m-1,2*m)=N(m,m);
JJ(2*m,2*m-1)=J(m,m);
JJ(2*m,2*m)=L(m,m);
end
form=1:
N1
forn=1:
N1
ifm==n
else
H(m,n)=-u(m)*u(n)*(G(m,n)*sin(delt(m)-delt(n))-B(m,n)*cos(delt(m)-delt(n)));
J(m,n)=u(m)*u(n)*(G(m,n)*cos(delt(m)-delt(n))+B(m,n)*sin(delt(m)-delt(n)));
N(m,n)=-J(m,n);L(m,n)=H(m,n);
JJ(2*m-1,2*n-1)=H(m,n);
JJ(2*m-1,2*n)=N(m,n);
JJ(2*m,2*n-1)=J(m,n);
JJ(2*m,2*n)=L(m,n);
end
end
end
form=1:
N1
PP(2*m-1)=pp(m);
PP(2*m)=qq(m);
end
uu=-inv(JJ)*PP';precision=max(abs(uu));
forn=1:
N1
delt(n)=delt(n)+uu(2*n-1);
u(n)=u(n)+uu(2*n)*u(n);
end
k
precision
PP'
JJ
delta=delt',
ua=u'
k=k+1;
end
forn=1:
N1+1
U(n)=u(n)*(cos(delt(n))+j*sin(delt(n)));
end
form=1:
N1+1
I(m)=Y(3,m)*U(m);
end
S3=U(3)*sum(conj(I))
form=1:
N1+1
forn=1:
N1+1
S(m,n)=U(m)*(conj(U(m))-conj(U(n)))*conj(-Y(m,n));
end
end
S
Y=
0-5.0000i0+2.5000i0+2.5000i
0+2.5000i0-2.5000i0
0+2.5000i00-2.5000i
k=
1
precision=
0.7200
ans=
0.2000
0.1000
-1.0000
-0.6000
JJ=
-5.000002.50000
0-5.000002.5000
2.50000-2.50000
02.50000-2.5000
delta=
-0.3200
-0.7200
0
ua=
0.8000
0.5600
1.0000
k=
2
precision=
6.0265
ans=
0.3930
-0.1699
-0.5639
-0.3524
JJ=
-2.93010.19301.0316-0.4361
0.1930-3.46990.43611.0316
1.03160.4361-1.03160.4361
-0.43611.03160.4361-0.5364
delta=
0.3291
2.8735
0
ua=
2.5844
3.9348
1.0000
k=
3
precision=
0.5239
ans=
12.4074
-48.2028
-15.2958
-60.3298
JJ=
14.908312.2074-21.022314.2958
12.2074-81.6973-14.2958-21.0223
-21.0223-14.295821.0223-14.2958
14.2958-21.0223-14.2958-98.4372
delta=
0.0131
2.6176
0
ua=
1.4016
1.8735
1.0000
k=
4
precision=
0.5705
ans=
3.5129
-11.8594
-4.3586
-15.0155
JJ=
2.13683.3129-5.64053.3586
3.3129-21.7819-3.3586-5.6405
-5.6405-3.35865.6405-3.3586
3.3586-5.6405-3.3586-23.1904
delta=
-0.2228
2.5623
0
ua=
0.8072
0.8047
1.0000
k=
5
precision=
0.7532
ans=
1.2127
-2.7117
-1.5668
-3.7406
JJ=
-0.44641.0127-1.52180.5668
1.0127-6.0698-0.5668-1.5218
-1.5218-0.56681.5218-0.5668
0.5668-1.5218-0.5668-4.7594
delta=
-0.4786
2.9307
0
ua=
0.5368
0.1986
1.0000
k=
6
precision=
3.9133
ans=
0.7475
-0.4066
-0.9295
-0.9557
JJ=
-0.93420.5475-0.2571-0.0705
0.5475-1.94740.0705-0.2571
-0.25710.07050.25710.0705
-0.0705-0.25710.0705-0.4543
delta=
-0.5744
6.8440
0
ua=
0.5957
-0.1079
1.0000
k=
7
precision=
7.1054
ans=
0.8635
-0.4918
-0.8543
-0.6969
JJ=
-1.18240.6635-0.0678-0.1457
0.6635-2.36590.1457-0.0678
-0.06780.14570.06780.1457
-0.1457-0.06780.1457-0.1260
delta=
-0.4387
13.9494
0
ua=
0.7139
-0.3692
1.0000
k=
8
precision=
1.9082
ans=
0.3197
-0.6688
-0.3617
-0.7772
JJ=
-1.77930.11970.1636-0.6383
0.1197-3.31680.63830.1636
0.16360.6383-0.16360.6383
-0.63830.16360.6383-0.5181
delta=
-0.3666
15.8576
0
ua=
0.8631
-0.6750
1.0000
k=
9
precision=
3.8953
ans=
1.6922
-0.3437
-1.7189
-0.4723
JJ=
-3.28111.49221.26670.7189
1.4922-4.1685-0.71891.2667
1.2667-0.7189-1.2667-0.7189
0.71891.2667-0.7189-1.0113
delta=
-1.0340
17.0108
0
ua=
-0.6076
1.9543
1.0000
k=
10
precision=
1.4194
ans=
1.0341
-4.5803
-3.1393
-12.2058
JJ=
2.83460.8341-2.05782.1393
0.8341-6.5260-2.1393-2.0578
-2.0578-2.13932.0578-2.1393
2.1393-2.0578-2.1393-21.1538
delta=
0.2624
18.4302
0
ua=
-0.0963
0.9623
1.0000
k=
11
precision=
8.2038
ans=
0.4085
-0.3590
-1.1461
-3.0950
JJ=
0.41260.2085-0.18000.1461
0.2085-0.5054-0.1461-0.1800
-0.1800-0.14610.1800-0.1461
0.1461-0.1800-0.1461-4.8100
delta=
3.7983
26.6340
0
ua=
0.0347
0.2557
1.0000
k=
12
precision=
56.8973
ans=
0.2364
0.0104
-0.9834
-0.7782
JJ=
0.08360.0364-0.0148-0.0166
0.0364-0.09560.0166-0.0148
-0.01480.01660.01480.0166
-0.0166-0.01480.0166-0.3417
delta=
6.2276
83.5313
0
ua=
0.4146
0.2296
1.0000
k=
13
precision=
4.2239
ans=
0.4824
0.1972
-1.2248
-0.8100
JJ=
-0.95680.2824-0.07820.2248
0.2824-0.7624-0.2248-0.0782
-0.0782-0.22480.0782-0.2248
0.2248-0.0782-0.2248-0.3419
delta=
5.4982
85.6972
0
ua=
0.3246
-0.7403
1.0000
k=
14
precision=
2.7770
ans=
1.3722
0.0941
-1.5985
-2.0233
JJ=
-0.52111.1722-0.05300.5985
1.1722-0.5328-0.5985-0.0530
-0.0530-0.59850.0530-0.5985
0.5985-0.0530-0.5985-2.7936
delta=
3.0271
82.9202
0
ua=
-0.3502
-0.2819
1.0000
k=
15
precision=
3.0040
ans=
0.0591
0.3034
-0.7591
-0.8519
JJ=
-0.8164-0.1409-0.0532-0.2409
-0.1409-0.40970.2409-0.0532
-0.05320.24090.05320.2409
-0.2409-0.05320.2409-0.4505
delta=
2.5351
85.9242
0
ua=
-1.3029
-0.1852
1.0000
k=
16
precision=
2.3451
ans=
2.6542
-5.7938
-1.5977
-0.7681
JJ=
-2.59412.4542-0.08230.5977
2.4542-14.3817-0.5977-0.0823
-0.0823-0.59770.0823-0.5977
0.5977-0.0823-0.5977-0.2539
delta=
2.6796
85.8294
0
ua=
-0.8328
0.2492
1.0000
k=
17
precision=
2.1511
ans=
0.6120
-1.5568
-0.4840
-0.8082
JJ=
-1.81060.4120-0.0530-0.5160
0.4120-5.12430.5160-0.0530
-0.05300.51600.05300.5160
-0.5160-0.05300.5160-0.3634
delta=
2.7015
87.9806
0
ua=
-0.7546
0.4517
1.0000
k=
18
precision=
1.7663
ans=
1.3790
-0.2753
-1.3753
-0.3450
JJ=
-2.47201.17900.76520.3753
1.1790-3.2226-0.37530.7652
0.7652-0.3753-0.7652-0.3753
0.37530.7652-0.3753-0.2552
delta=
2.4779
87.0611
0
ua=
-0.3927
-0.3462
1.0000
k=
19
precision=
2.4657
ans=
0.8856
-0.2280
-1.0807
-1.2298
JJ=
-0.44320.6856-0.33020.0807
0.6856-1.0992-0.0807-0.3302
-0.3302-0.08070.3302-0.0807
0.0807-0.3302-0.0807-0.9294
delta=
2.0833
89.5268
0
ua=
-0.3213
0.1755
1.0000
k=
20
precision=
6.7584
ans=
0.9702
-0.1445
-1.0702
-0.7993
JJ=
-0.27160.7702-0.12230.0702
0.7702-0.7607-0.0702-0.1223
-0.1223-0.07020.1223-0.0702
0.0702-0.1223-0.0702-0.2763
delta=
2.5405
96.2852
0
ua=
-0.4498
-0.6442
1.0000
k=
21
precision=
1.0095
ans=
0.4867
0.6504
-0.6507
-1.0029
JJ=
-1.56220.28670.6347-0.3493
0.2867-0.46130.34930.6347
0.63470.3493-0.63470.3493
-0.34930.63470.3493-1.4404
delta=
3.0560
96.0991
0
ua=
-0.9040
-0.3727
1.0000
k=
22
precision=
0.9630
ans=
-0.3931
-1.4325
-0.2136
-0.6457
JJ=
-2.5532-0.59310.3016-0.7864
-0.5931-5.61820.78640.3016
0.30160.7864-0.30160.7864
-0.78640.30160.7864-0.3930
delta=
2.7823
97.0621
0
ua=
-0.8600
-0.6688
1.0000
k=
23
precision=
1.2868
ans=
1.0020
-0.1481
-1.0461
-0.2811
JJ=
-3.44980.80201.43710.0461
0.8020-3.9459-0.04611.4371
1.4371-0.0461-1.4371-0.0461
0.04611.4371-0.0461-0.7992
delta=
2.5570
96.1676
0
ua=
-0.3943
0.1918
1.0000
k=
24
precision=
3.9311
ans=
0.8564
-0.0072
-1.1125
-0.8439
JJ=
-0.67010.6564-0.15200.1125
0.6564-0.8845-0.1125-0.1520
-0.1520-0.11250.1520-0.1125
0.1125-0.1520-0.1125-0.3359
delta=
2.1559
100.0987
0
ua=
-0.3409
-0.5566
1.0000
k=
25
precision=
13.4887
ans=
0.6611
-0.4139
-0.7507
-1.7781
JJ=
-0.06710.4611-0.4035-0.2493
0.4611-1.09480.2493-0.4035
-0.40350.24930.40350.2493
-0.2493-0.40350.2493-1.9526
delta=
15.6446
111.7163
0
ua=
-3.4301
1.1257
1.0000
k=
26
precision=
0.4501
ans=
-8.6032
-47.7527
8.3458
-1.3506
JJ=
-10.9756-8.80322.4175-9.3458
-8.8032-106.68109.34582.4175
2.41759.3458-2.41759.3458
-9.34582.41759.3458-3.9188
delta=
15.5860
111.7421
0
ua=
-1.8863
0.6517
1.0000
k=
27
precision=
0.3809
ans=
-2.1262
-11.9926
1.8998
-0.6441
JJ=
-5.6984-2.32621.0176-2.8998
-2.3262-29.88352.89981.0176
1.01762.8998-1.01762.8998
-2.89981.01762.8998-1.1058
delta=
15.4910
111.9370
0
ua=
-1.1679
0.5392
1.0000
S3=
0.6285-0.3512i
S=
01.2745+2.4855i-0.6285+0.5586i
-1.2745-0.1974i00
0.6285-0.3512i00
>>