2 一維瞬態常係數導熱離散方程時間隱式格式

2021-10-10 06:14:11 字數 1148 閱讀 7442

#include double length;  //導熱體的長度

double tl,tr; //左側溫度,右側溫度

double den,c,k; //物質密度,比熱,導熱係數

double s; //熱源

double dt; //時間間隔

int n; //控制體的個數

int tmaxstep;

//-------------以上為和物理問題有關的變數--------------------

double dx;

double t0[103],t[103]; //從0-10共十乙個元素,即左右邊界兩個體積為零的控制體和中間9個控制體,內部的元素對應11個溫度的值

double ae0[103],aw0[103],ap0[103],ap1[103],b[103];

double ae1[103],aw1[103];

double a[103],b0[103],c0[103],d[103],x[103],p[103],q[103]; //前三為係數矩陣係數,d為等式右側矩陣係數,x為未知數

//三對角子函式

void tdma()

}//---------------------------------

int main()

{ int tstep; //時間步迴圈控制變數

int i; //控制體迴圈變數

length=3;

tl=3;tr=5;

den=100;c=1000;k=10;

s=10;

dt=100;

n=3; //當n值取得很大時,t值不變會引起非物理現象,即ap0為負數

tmaxstep=2000;

//----------------以上為計算過程中需要的變數,並給已知量進行賦值-----------------------------

dx=length/n;

/* for(i=2;i在此處開啟powershell視窗;

3)執行語句:

1dtr.exe >> t.txt

如果以上命令不行可以用.\1dtr.exe >> t.txt

4)完成,控制台內容儲存到日誌檔案t.txt中

*/

c語言一維陣列(2)

3.陣列引用 1 原則 通過下標引用陣列中每乙個單元資料 2 下標範圍 0 n 1 n為陣列長度 3 格式 下標必須為常量或常量表示式 4 例 int a 5 a 0 1 a 5 3 2 a a 1 3 陣列下標越界不是語法錯誤,編譯器無法檢查出來,但沒有意義 4.陣列的初始化 1 定義同時賦值 i...

《Python自動化運維 技術與最佳實踐》一導讀

為什麼要寫這本書 隨著資訊時代的迅速發展,尤其是網際網路日益融入大眾生活,作為這一切背後的it服務支撐,運維角色的作用越來越大,傳統的人工運維方式已經無法滿足業務的發展需求,需要從流程化 標準化 自動化去構建運維體系,其中流程化與標準化是自動化的前提條件,自動化的最終目的是提高工作效率 釋放人力資源...

陣列一維2維3維動態申請及用Vector的表示方法

1 變長一維陣列 這裡說的變長陣列是指在編譯時不能確定陣列長度,程式在執行時需要動態分配記憶體空間的陣列。實現變長陣列最簡單的是變長一維陣列,你可以這樣做 include using namespace std int main 注意int p new int len 這一句,你不能這樣做 int ...