頭驗目的：一個二維靜電場，電位函數(shù)為＜p(x,y),邊界條件如圖所示，將正方形場域分成25個網(wǎng)格，有16個內(nèi)部網(wǎng)格點，假定16個網(wǎng)格點的初始值都為0，試用超松弛法確定16個內(nèi)網(wǎng)格的點位值。余7一4_丄_二P二P:二余7一4_丄_二P二P:二100V實驗程序用c語言寫:#include"stdio.h"#include"stdio.h"#include"math.h"Main()//頭文件//引用math頭文件//主函數(shù){inti,u;int{inti,u;intm,n,w;floata[6][6];floatp,q;w=0;for(i=0;i<6;i++){for(u=0;u<6;u++){a[i][u]=0;}}for(i=0;i<6;i++)//定義循環(huán)體變量i，u〃定義循環(huán)體變量m,n,定義計數(shù)變量w//定義浮點型二維數(shù)組a[6][6];〃定義浮點型變量p,q；//計數(shù)器賦初值//嵌套循環(huán)體給二維數(shù)組賦初值{a[0][i]=100;a[5][i]=50;a[i][5]=100;}do{p=a[4][4];for(m=1;m<5;m++){for{a[0][i]=100;a[5][i]=50;a[i][5]=100;}do{p=a[4][4];for(m=1;m<5;m++){for(n=1;n<5;n++)//定義頂部邊界電勢為100//定義底部邊界電勢為50//定義右邊邊界電勢為100〃設定dowhile循環(huán)〃用a[4][4]作為比較器//此處嵌套循環(huán)體對象是中間的16個點{a[m][n]=(a[m-1][n]+a[m+1][n]+a[m][n-1]+a[m][n+1])/4;}}q=abs(a[4][4]-p);w=w+1;}while(q>0.000001);for(i=0;i<6;i++)//有限差分法〃每循環(huán)一次，用新的a[4][4]與循環(huán)之前的a[4][4]進行對比，并把絕對值賦給q變量//每循環(huán)一次，計數(shù)器加一//當誤差小于0.000001時跳出循環(huán)，誤差控制在百萬分之一{for(u=0;u<6;u++){printf("%f",a[i][u]);} //打印輸出整個矩陣printf("\n");printf("%d",w);//打印輸出計數(shù)器，共進行了w次循環(huán)程序運行結果:C程序運行結果:

OOTOSOSosOSOS9OOTZiZEZLEOSE'CG99T6'ZfrEBTSTE0￡OOT8SZ9T8TB9S'OS09S€K0VOOTgOTE'ASSTElzfrAZSES'SSSASO'SE0EOOTBTBTE699TrSBOESOTZ￡Z￡ZYS0EOOTOOTOOTOOTOOTOOTI9■V￡EIojqnopgxgnpaujeuun田un:s>c|dpi|SAc|\]生 ￡￡呂日:4》暫罰」o/Vk辰\|即叫皿[7ppv MOHP珂n沖ol|SAno|\sja5n\：3:J3*P|CJjUBJJnQ越丨m保嶂丨心3■龜覽□:d|BHMopu麗do琲總召罔|目」它召Ein耳總口53iqde」芋ms氐衛(wèi)卩了引日:郢糜刼審日VI丄VIAI結果分析：1?通過C語言程序設計，使用循環(huán)體實現(xiàn)了有限差分法通過設定q變量，每次與循環(huán)之前對比，使用dowhile語句實現(xiàn)9Njj屮N-ljJvW保證了迭代誤差不超過允許范圍?？梢造`活i,j-i,j設置q的允許值來控制誤差通過設定w變量，實現(xiàn)計數(shù)功能，本次仿真，誤差控制在在1/1000000內(nèi)，共循環(huán)35次。4?將C程序得出的值輸入MATLAB，繪圖仿真出電磁場。總結和體會：通過這一