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語言程式?首先,你...