next up previous
Next: A. プログラムの解説 Up: 4 実行の仕方 Previous: 4.3 計算結果をファイルに保存する

4.4 固有値を調べる

save() した変数の内容は、 load() することでメモリー上に復元出来る。 例えば
» load('free_0.3_1280.mat');
とすると、N, poisson, d1280, e1280 が戻る。

free_0.3_1280.mat は 2GB 以上あるファイルなので、 読み込みには少し時間がかかる。

analyze_ev_03.m
load('free_0.3_10.mat');
load('free_0.3_20.mat');
load('free_0.3_40.mat');
load('free_0.3_80.mat');
load('free_0.3_160.mat');
load('free_0.3_320.mat');
load('free_0.3_640.mat');
load('free_0.3_1280.mat');
ev03=zeros(200,8);
ev03(:,8)=sort(diag(d1280));
ev03(:,7)=sort(diag(d640));
ev03(:,6)=sort(diag(d320));
ev03(:,5)=sort(diag(d160));
ev03(:,4)=sort(diag(d80));
ev03(:,3)=sort(diag(d40));
ev03(:,2)=sort(diag(d20));
ev03(1:100,1)=sort(diag(d10));
error03=ev03(:,2:8)-ev03(:,1:7);
save('eigen_only_03.mat', 'poisson', 'ev03', 'error03');

% 誤差減少の様子のグラフを描く
x=1:7;
semilogy(x,error03(4:40,:));

% 固有値の表 (N=1280の場合と誤差の目安N=1280,640の差)
ev_and_error03=zeros(200,2);
ev_and_error03(:,1)=ev03(:,8);
ev_and_error03(:,2)=error03(:,7);
ev_and_error03

» clear
» analyze_ev_03

ev_and_error03 =

        0.0013649376407771       0.00131804197613936
       0.00252915430712694        0.0024748440140326
       0.00325771474892962       0.00319239425393388
          181.394289625725       0.00344357463796996
          384.008733137516        0.0081658839260399
          589.042437123803        0.0113128114106757
          1211.09781000587        0.0212640410970835
          1211.09799702052        0.0214435521793348
          3732.34385261567         0.127309501172931
          3732.34425608862         0.127707324018047
          4055.89733718611        0.0879733601450425
          4797.64315733848         0.168206207419644
          5955.41240760091         0.186384941122014
          11121.7657408117         0.389936375522666
          11121.7658466479         0.390039700629131
           13714.172347162         0.811331795957813
          14992.3991893111         0.854916557067554
          17283.8066483796         0.970440969020274
          17283.8069222436         0.970708716016816
          23361.2214152101         0.999834179459867
          26083.3276248129          1.48965588299325
          28386.0371201005           1.4969686598015
           39310.360231204          3.45407351025642
          39310.3603609914          3.45420124251541
          41672.4812033997          3.71575093670981
            45771.96441852          3.87184302660171
          46402.3573217012          2.76687474390201
          46402.3573929108          2.76694505955675
           59010.201221058          4.88540001578804
          59010.2013349507          4.88551071860275
          79021.4588651732          5.67223714283318
          85182.5492425693          7.26532158453483
          86266.3258468876          10.8614341835782
          89110.8903863748          7.18668166708085
          89953.8014098738          10.9435912937624
          94142.9804709299          11.5854896803212
          94142.9806637669          11.5856804827199
          113184.534771899          13.5832112985227
          118845.562796414          13.7373521245318
          132483.357334267          11.9333889349946
          132483.357458445          11.9335115564463
          154502.637279212          17.9538902596105
          154502.637448147          17.9540589544922
          172603.205649117           28.455339094915
           172603.20568244          28.4553709959437
          175446.231584341          29.5457886184158
          183982.452198466          29.8637832001841
          202671.013491451           21.258094235498
          209398.755250829          33.3966757245944
          209398.755312708          33.3967338011716
          212597.495006249           25.367056076735
          218511.287627452          25.1077884892002
          257300.254592629          40.1929241697362
          265227.706948033           40.155794010323
           301615.69016376           65.450553792296
          304234.251293826          39.9083605690976
          304234.251409618          39.9084718467202
          308843.695026233          65.8211032238905
          315222.159033523          65.7414655012544
          315222.159175893          65.7416057654191
          340861.508010279          52.2170908502885
            340861.5080983          52.2171748431865
          348738.466003743           72.942456478253
          358994.366639234          73.6771637507482
          419704.205321572          84.1163488712627
          419704.205363347          84.1163893595804
          432104.023789524          61.9977013702155
          448111.341829276          71.4046083428548
          455965.418132773          71.4308832453098
          504634.116430834          135.696183689579
          504634.116466369          135.696218073484
           508160.78800335          138.175009226252
          514150.300607866          101.624177543446
          524777.019267532           101.13336351153
          525427.988618664          138.728924912633
          564402.506314476            147.3644847041
          564402.506347273          147.364514835528
          607789.318710032          103.214412132627
          607789.318856048          103.214549056836
          641333.033987161          163.424287841306
          654467.803352731          164.001635291846
          659188.365490078          128.792643907247
          659188.365648636          128.792798820301
          768422.310670925          193.086471016053
          768422.310790811          193.086585524608
           780531.36703622          259.381981076673
          792546.454849486          258.928816297092
          804256.037450477          260.146224832395
          804256.037537475          260.146310563898
          817611.079818579            154.6750980051
          840830.516183255          174.472524670418
          848652.098908719          191.766743696993
          854113.933689159          274.019409974455
          875676.327923699          258.257123250514
          932911.829440303          228.237137800781
          946511.095847984          227.158707594732
          962442.047892718          300.481696763309
          962442.047919214          300.481718725059
           1094545.4296956          239.937561749481
          1094545.42985362          239.937697523739
          1103469.15688749          336.423763490049
          1119916.46693493          336.414477792103
          1164122.03582189          302.081983955577
          1164122.03594479          302.082103944384
          1177361.39197426          447.988285851432
          1177361.39199494          447.988303680904
           1178941.6979411          471.518023404991
          1205687.86148877          473.418976456858
          1263972.77656606          486.819185405271
          1263972.77658351          486.819198070792
          1316580.48816828           392.58198495442
          1316580.48835604          392.582152920309
            1377083.519074          522.213879709598
          1394854.04056731          519.579673970584
          1417428.75785089          341.814259669511
          1449724.64227334          372.478431772441
          1467975.80916944          374.343155370792
          1562714.97466722          574.606707576429
          1562714.97473368          574.606766778277
          1571964.71023319          467.980706463102
          1588044.47145735          469.112491414649
           1681526.7354997          793.193972618319
          1702948.23858199          790.450701476773
          1712311.36762524          800.216244099196
          1712311.36772334          800.216339084553
          1787094.20459172          824.748921256047
          1793740.74541712          645.518612759421
          1812678.47124255          827.959851702675
          1813851.43055642          644.511388438288
          1829912.53382023          502.102205586154
          1829912.53391886          502.102297396399
           1920937.4777785          598.264840957941
          1920937.47791107          598.264967764262
          1943633.19966898          848.539444791386
          1943633.19978467          848.539557423908
          2116520.49852476          746.270226002205
           2116520.4985789          746.270276872441
          2139132.25741397          939.041253651027
          2161903.06769487          939.561404744629
          2298468.86795587           695.05418825848
          2341917.94824191          736.885544531513
           2362076.8942592          1285.83531035203
          2362076.89434457          1285.83539149584
          2362791.28140998          733.161924201995
          2364093.82400985          1304.59300984489
           2403233.6584324          1291.83177950885
          2430101.56330805          1056.20867442572
          2430101.56345172          1056.20881554438
          2490121.48398548          1325.23431929946
          2490121.48401201          1325.23434512643
          2500002.67967977          890.559409231413
          2522456.54589562          901.393844678532
          2642039.94794139          1396.41796060745
          2673147.42617723          1392.60824018018
          2779978.68258512          1168.14009857411
          2803861.47537892          1166.47869716678
          2880008.22459064          1085.32515348261
          2880008.22481691          1085.32537928689
          2900902.88705633          1365.83627556683
          2900902.88710822          1365.83631930966
          3002572.72507619          1088.17355911527
          3002572.72526047          1088.17374143843
          3199195.94487165          2034.78399233101
          3207914.54046035          1613.39477163414
          3230203.68868303          2030.78775129002
          3235270.86203712          1613.37468295498
          3236260.49722868          1805.42431836715
           3236260.4976304          1805.42471203301
          3252571.91919109          1585.44747315627
           3252571.9193132          1585.44759125123
          3344104.60860356          2086.99229426775
          3376658.85696765          2090.26340351207
          3542684.06584957          1291.12049713265
          3551430.88941027          2159.26846531546
          3551430.88961228          2159.26866538916
          3594929.08291571          1360.45273380075
          3619331.60375038          1354.95049784612
          3640072.77936095          1794.20450076321
          3640072.77940504          1794.20454354072
          3791936.05491453          1620.77019608021
           3808436.3115273          1752.08849012572
          3818141.36429521          2257.69092766475
          3863406.29007746          2120.54427834461
          4138771.69685334          2018.37305933097
          4164767.40287664           2031.7574648615
          4198638.16841503          2472.16444041394
          4198638.16874922          2472.16477098875
          4279248.60904028          3110.32035160065
          4282748.30440861          3030.52658050973
          4282748.30444707          3030.52661624085
          4339449.80930016          3094.37180056237
          4343133.85555124          1789.19096419308
          4343133.85595667          1789.19136499055
          4443028.19333915          3046.74781952798
          4443028.19338985          3046.74786863104
          4493959.81771926          2033.12655983679
           4493959.8178311          2033.12667025905
          4650158.54006466          3267.60115485732
          4660680.34175147          2664.27756253909
          4688987.09193959          3273.99175896775

>>


next up previous
Next: A. プログラムの解説 Up: 4 実行の仕方 Previous: 4.3 計算結果をファイルに保存する
桂田 祐史
2014-05-27