C語言課程設(shè)計電場線

C語言課程設(shè)計PAGE1電場線一目的設(shè)計程序，模擬兩個等量異號電荷電場線，并在計算機(jī)上顯示出來。二硬件環(huán)境CPU(PentiumDual-CoreE5400)，內(nèi)存（2GB），硬盤（500GB），顯示器(ASUS)等。三軟件環(huán)境WindowsXP，MyTC5.4.1計算相應(yīng)位置場強(qiáng)Ex，Ey，E確定兩正負(fù)電荷的位置計算相應(yīng)位置場強(qiáng)Ex，Ey，E確定兩正負(fù)電荷的位置開始以兩電荷為圓心，分別畫圓電場線條數(shù)i=0，1，2，…29由A=Pi/15，計算電場線的初始位置x1=200+10*cosA，y1=250-10*sinA計算下一位置x2=x1+Ex/E，y2=y1+Ey/E畫線（生成電場線）x1=x2，y1=y2結(jié)束(x-x2)*(x-x2)+(y-y2)*(y-y2)>r*r五調(diào)試情況1、有一些輸入錯誤，經(jīng)提示找到語法錯誤并改正。2、正電荷水平向左發(fā)出的電場線無法回到負(fù)電荷是圖像不太對稱，故用if(A!=PI)將其除去。六程序運(yùn)行結(jié)果七心得體會通過這次課程設(shè)計，我了解學(xué)習(xí)了C語言繪圖函數(shù)及其功能，復(fù)習(xí)了C語言編程的基本知識，運(yùn)用這些函數(shù)，按照老師要求基本完成了這次課程設(shè)計的任務(wù)，雖然過程中出現(xiàn)了許多錯誤，經(jīng)過自己查閱相關(guān)資料和在同學(xué)的相互幫助下，終于完成了，有一種成就感。我更認(rèn)識到合作的重要性，做事要有耐心，有想法要勇于嘗試，才能得到正確結(jié)果，同時要注重理論聯(lián)系實(shí)際學(xué)以致用。八附錄（源程序代碼）#include<stdio.h>#include<math.h>#include<graphics.h>#definePI3.141main(){intgd=DETECT,gm;/*進(jìn)行圖形初始化*/doublex1,x2,y1,y2,r1,r2;doubleEx1,Ex2,Ey1,Ey2,Ex,Ey,E,A=0.0;inti,j;initgraph(&gd,&gm,"");setbkcolor(0);setlinestyle(1,0,1);line(300,0,300,500);/*畫y坐標(biāo)軸(虛線)*/setlinestyle(0,0,1);setfillstyle(1,RED);circle(200,250,10);floodfill(200,250,15);line(195,250,205,250);line(200,245,200,255);/*畫正電荷*/setfillstyle(1,2);circle(400,250,10);floodfill(400,250,15);line(395,250,405,250);/*畫負(fù)電荷*/for(i=0;i<30;i++)/*定義電場線條數(shù)*/{if(A!=PI)/*除去與x軸夾角為PI的電場線*/x1=200+10*cos(A);y1=250-10*sin(A);for(j=0;j<30000;j++){r1=sqrt(pow(x1-200,2)+pow(250-y1,2));Ex1=(x1-200)/pow(r1,3);Ey1=(y1-250)/pow(r1,3);/*計算正電荷電場線在(x1,y1)處場強(qiáng)的x軸和y軸分量*/r2=sqrt(pow(400-x1,2)+pow(250-y1,2));Ex2=(400-x1)/pow(r2,3);Ey2=(250-y1)/pow(r2,3);/*計算負(fù)電荷電場線在(x1,y1)處場強(qiáng)的x軸和y軸分量*/Ex=Ex1+Ex2;Ey=Ey1+Ey2;E=sqrt(pow(Ex,2)+pow(Ey,2));/*計算總場強(qiáng)*/x2=x1+Ex/E;y2=y1+Ey/E;line(x1,y1,x2,y2);x1=x2;y1=y2;if(pow(x1-400,2)+pow(y1-250,2)<=100){break;}/*生成電場線并設(shè)置電場線的終止于負(fù)電荷*/}A=A+PI/15;}while(!kbhit());closegraph();}波的合成一目的設(shè)計程序，模擬兩列波的疊加，并在屏幕上顯示出來。二硬件環(huán)境CPU(PentiumDual-CoreE5400)，內(nèi)存（2GB），硬盤（500GB），顯示器(ASUS)等。三軟件環(huán)境操作系統(tǒng):WindowsXP 軟件：MyTC5.4.1開始輸入a1,a1,w1,w2,v1,v2的值開始輸入a1,a1,w1,w2,v1,v2的值定義波形個數(shù)t=1，2,…..,N定義波的傳播距離x=0,1,2,……,999計算y1=a1*cos(w1*(t-x/v1))+80計算y2=a2*cos(w2*(t+x/v2))+200計算y=y1+y2+100畫線從（x,y10）到（x+1,y1）畫線從（x,y20）到(x+1,y2)畫線從（x,y0）到(x+1,y)結(jié)束五程序調(diào)試情況1、有一些輸入錯誤，經(jīng)提示找到語法錯誤并改正。2、設(shè)置顏色函數(shù)位置有誤，通過查閱越資料，我們將該函數(shù)放在所畫線前，該問題得到了解決。六程序運(yùn)行結(jié)果七心得體會通過這次課程設(shè)計，我了解學(xué)習(xí)了C語言繪圖函數(shù)及其功能，復(fù)習(xí)了C語言編程的基本知識，運(yùn)用這些函數(shù)，按照老師要求基本完成了這次課程設(shè)計的任務(wù)，雖然過程中出現(xiàn)了許多錯誤，經(jīng)過自己查閱相關(guān)資料和在同學(xué)的相互幫助下，終于完成了，有一種成就感。我更認(rèn)識到合作的重要性，做事要有耐心，有想法要勇于嘗試，才能得到正確結(jié)果，同時要注重理論聯(lián)系實(shí)際學(xué)以致用。八附錄（源程序代碼）#include<stdio.h>#include<math.h>#include<graphics.h>#defineN100main(){intx,x0,y,y0,y1=0,y2=0,y10,y20,t=0;doublea1,a2,w1,w2,v1,v2;intgd=DETECT,gm;printf("Pleaseinputa1a2w1w2v1v2:");scanf("%lf%lf%lf%lf%lf%lf",&a1,&a2,&w1,&w2,&v1,&v2);/*輸入數(shù)據(jù)*/initgraph(&gd,&gm,"");setbkcolor(BLACK);setcolor(WHITE);/*設(shè)置背景和線條顏色*/for(t=0;t<N;t++)/*設(shè)置動態(tài)波形*/{delay(30000);cleardevice();setlinestyle(1,0,1);setcolor(8);line(0,80,1000,80);line(0,200,1000,200);line(0,380,1000,380);setlinestyle(0,0,1);y10=a1*cos(w1*(t-1/v1))+80;y20=a2*cos(w2*(t-1/v1))+200;y0=y1+y2+100;for(x=0;x<1000;x++)/*定義波傳播距離*/{y1=a1*cos(w1*(t-x/v1))+80;y2=a2*cos(w2*(t-x/v2))+200;y=y1+y2+100;setcolor(3);line(x,y10,x+1,y1);setcolor(3);line(x,y20,x+1,y2);setcolor(5);line(x,y0,x+1,y);/*生成波形圖*/y10=y1;y20=y2;y0=y;}}while(!kbhit());closegraph();}分子熱運(yùn)動一目的1、學(xué)會數(shù)值求解運(yùn)動方程的verlet方法。2、掌握分子熱運(yùn)動的基本思想，并使用verlet方法編寫分子動力學(xué)基本程序，模擬微小團(tuán)簇的熱運(yùn)動，并輸出。二硬件環(huán)境CPU(PentiumDual-CoreE5400)，內(nèi)存（2GB），硬盤（500GB），顯示器(ASUS)等。三軟件環(huán)境操作系統(tǒng):WindowsXP 軟件：MyTC5.4.1開始定義七個粒子初始位置令初始速度均為零開始定義七個粒子初始位置令初始速度均為零計算初始位置（a=0）時的位置、速度、作用力、總勢能調(diào)用子函數(shù)Fxyz()、Fv()、Ff()、U()計算下一位置（a++）的相應(yīng)參數(shù)輸入n的值結(jié)束a<n？是是否否否否五原理在一個孤立系統(tǒng)中，N個粒子給予初始位置，初始速度。有牛頓第二定律可以解出和然后由verlet公式（h為時間間隔，令為0.00466.）可以由n時刻求出n＋1時刻的位置和速度。兩粒子間的勢能只與這兩粒子的距離有關(guān)：。由公式可求出其中是參量，是兩粒子間距離。分別表示第n時刻第i個粒子在x，y，z方向受到的合力。通過約化單位制＝1，＝1，m＝1。把以上三個公式代入varlet公式即可求出n時刻每個粒子的位置和速度。六程序調(diào)試情況1、有一些輸入錯誤，經(jīng)提示找到語法錯誤并改正。2、起初使用遞歸方法經(jīng)常益處，最后用直接賦值，避免遞歸的方法解決了該問題。七程序運(yùn)行結(jié)果八心得體會通過這次課程設(shè)計，我了解學(xué)習(xí)了C語言繪圖函數(shù)及其功能，復(fù)習(xí)了C語言編程的基本知識，運(yùn)用這些函數(shù)，按照老師要求基本完成了這次課程設(shè)計的任務(wù)，雖然過程中出現(xiàn)了許多錯誤，經(jīng)過自己查閱相關(guān)資料和在同學(xué)的相互幫助下，終于完成了，有一種成就感。我更認(rèn)識到合作的重要性，做事要有耐心，有想法要勇于嘗試，才能得到正確結(jié)果，同時要注重理論聯(lián)系實(shí)際學(xué)以致用。九、附錄（源程序代碼）#include<stdio.h>#include<math.h>#include<graphics.h>#defineN7#defineh0.005466doubleVx[N],Vy[N],Vz[N],Fx[N],Fy[N],Fz[N],rij,u,n;doublex[N]={-1.343659,0.393647,0.756033,0.741549,-0.034382,-0.354385,0.040196};doubley[N]={0.445998,-0.328215,-0.465412,0.556876,1.312412,-0.539929,0.318270};doublez[N]={1.419669,1.712890,0.663500,1.155653,0.807389,0.927735,0.349577};/*定義并設(shè)定七個分子的初始位置*/inti,j;doublea;voidFxyz();/*定義位置坐標(biāo)函數(shù)*/voidFf();/*定義作用力函數(shù)*/voidFv();/*定義速度函數(shù)*/voidU();/*定義系統(tǒng)總勢能函數(shù)*/main()/*主函數(shù)*/{intgd=DETECT,gm;printf("Pleaseinputn(n>10000):");scanf("%lf",&n);initgraph(&gd,&gm,"");for(a=0;a<n;a++){Fxyz();/*位置、速度、作用力子函數(shù)*/Fv();Ff();}U();/*系統(tǒng)總勢能子函數(shù)*/for(i=0;i<N;i++){printf("x=%9lf,y=%9lf;z=%9lf\n",x[i],y[i],z[i]);}printf("u=%lf\n",u);while(!kbhit());closegraph();}voidFxyz(){for(i=0;i<N;i++){setcolor(1);circle(300+20*x[i],300+20*y[i],2);/*畫七個分子并設(shè)定其大小*/x[i]+=h*Vx[i]+h*h*Fx[i]/2;/*計算每個分子下一位置的坐標(biāo)*/y[i]+=h*Vy[i]+h*h*Fy[i]/2;z[i]+=h*Vz[i]+h*h*Fz[i]/2;setcolor(15);/*分子顏色設(shè)置為白色*/circle(300+20*x[i],300+20*y[i],2);floodfill(300+20*x[i],300+20*y[i],15);}}voidFf()/*計算每個分子受到其他六個分子的作用力*/{for(i=0;i<N;i++){Fx[i]=0.0,Fy[i]=0.0,Fz[i]=0.0;for(j=i+1;j<N;j++){rij=sqrt(pow(x[i]-x[j],2)+pow(y[i]-y[j],2)+pow(z[i]-z[j],2));Fx[i]+=48*(x[i]-x[j])*(pow(1/rij,14)