c語言差分法求解波動方程程式,c語言差分法求解波動方程程式

2021-03-04 05:35:27 字數 1303 閱讀 6504

1樓:挪威的森林徐浩

第二個邊界條件ux(1,t)=exp(-t) 是不是寫錯了?應該是u(1,t)=exp(-t)吧。。

二維波動方程怎麼用差分法求解?就是把位移場或者波動場用波動方程表示出來,急求啊。 32

2樓:阿三五三四

若那個t=0的點為0.05則方程為 0.1cos(4π/5t+2π/5x-1π/30)單位自己腦補lhrd

二維波動方程的有限差分程式(詳細的matlab或者fortran程式) 15

3樓:鬥破了啊

^dx=8000;

dy=4000;

ix=400;

iy=80;

nx=dx/ix+1;

ny=dy/iy+1;

vel=1500*ones(ny,nx);

vel((fix(ny/3)):(fix(2*ny/3)))=1500;

vel((fix(2*ny/3)):ny)=1700;

fre=80;

it=0.01;

u_0=zeros(ny,nx);

u_1=zeros(ny,nx);

for i=1:nx;

for j=1:ny;

u_1(1,i)=-(it)^2*sin(2*pi*fre*it)*exp(-2*pi*fre*(it));

endend

u_1(1,1) = u_1(1,2);

u_1(1,nx) = u_1(1,(nx-1));

for k=2:n-1

for j=1:ny

for i=2:nx-1

if j=1

u_2(j,i) = ((it*vel(j,i))^2)/(ix^2)*(u_1(j,(i-1))+u_1(j,(i+1))-2*u_1(j,i)) + ((it*vel(j,i))^2)/(iy^2)*(u_1((j-1),i)+u_1((j+1),i)-2*u_1(j,i)) + 2*u_1(j,i) - u_0(j,i);

endu_2(1,i)= ((it*vel(1,i))^2)/(ix^2)*(u_1(1,(i+1))+u_1(1,(i-1))-2*u_1(j,j)+ ((it*vel(1,i))^2)/(iy^2)*2*(u_1(2,i)-u_1(1,i))+ 2*u_1(1,i)-u_0(1,i)- (it)^2*sin(2*pi*fre*k*it)*exp(-2*pi*fre*(it*k));

endelseif j=ny

u_2(i,j-1)=u_1(i,j)end

求解答c語言程式題,求解答C語言程式題

內容來自使用者 kcd9272 三 程式題 1 程式設計序解百雞問題 100元錢買100隻雞,公雞每只3元,母雞每只5元,小雞3只1元,問公雞 母雞 小雞各多少只?include int main 2 求出100 999的所有的水仙花數 水仙花數既xyz x x x y y y z z z incl...

C語言求解,C語言求解

我去,你也太大意了吧!1 定義的bml變數是小寫的 float w,h,bmi 在後面各處使用該變數的時候卻是大寫的 bmi fbmi w,h 2 scanf函式名你又打多了乙個字母t scantf f,f w,h 3 定義的函式fbmi,呼叫函式的時候同樣大小寫不匹配 bmi fbmi w,h 4...

C語言程式程式設計,如何編寫C語言程式

把你的思路想好,想要達到什麼樣的功能,思路清晰後就可以一個一個功能的去寫,出錯了不用急,想一下自己錯在什麼地方,好好的看一下 雖然這樣做會花很多 的時間,但是這樣做會提高自己的閱讀分 能力,啊?書上不是有麼 自己學著寫啊 10.include main printf n 如何編寫c語言程式?首先,你...