 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 357 | 回復(fù): 0 | |||
hawaiicn新蟲 (初入文壇)
|
[交流]
請(qǐng)大家?guī)臀铱聪挛业腗ATLAB卷積分問題出在哪里？謝謝！
|
|
用MATLAB自帶卷積分函數(shù)計(jì)算彈性半空間在一三角脈沖荷載下自由表面的豎向位移，與文獻(xiàn)正確結(jié)果相差10的5次方倍，但自己編寫卷積分代碼，反而得到文獻(xiàn)結(jié)果。請(qǐng)大家?guī)臀铱聪聠栴}出在哪？謝謝！ 另外，為何MATLAB自帶卷積分函數(shù)結(jié)果矩陣的維數(shù)會(huì)是被卷積兩矩陣維數(shù)之和減1？ 1. 問題描述，如下圖1所示： 2. 格林函數(shù)，如下圖2所示： 3. 參考文獻(xiàn)正確解，如下圖3所示： 4. 本人直接用自編卷積分MATLAB代碼算的出解（如下圖4）及相應(yīng)代碼： clear; clc; %loading history dt=0.01; ti=0.0008; te=12; t=ti:dt:te; pt=(10*t./1.5).*(t>=0 & t<=1.5)+(20-10*t./1.5).*(t>1.5 & t<=3)+0.*(t>3); m=length(pt); %load distribution in space rp=0.1; ri=0; rc=0.1; dr=0.001; r=ri:dr:rc; pr=1*(r<=rp)+0*(r>rp); n=length(pr); %load function with respect to t and r p=pr.'*pt; %green's function G=1; cs=1; for i=1:1:m for j=1:1:n u(j,i)=heaviside(cs*t(i)-r(j))/(pi*G*sqrt(t(i)^2-(r(j)/cs)^2)); end end %convolution and response of displacement for i=1:1:m for j=1:1:n v(j,i)=0; for k=1:1:i for g=1:1:j v(j,i)=v(j,i)+p(g,k)*u(j-g+1,i-k+1)*dr*dt; end; end; end; end; %plot the response history plot(t,v(n,  );xlabel('t/s'); ylabel('v/m'); grid on; title('response of point A or C'); 5. 直接用MATLAB自帶卷積分算的出解（如下圖5）及相應(yīng)代碼： clear; clc; %loading history dt=0.01; ti=0.0008; te=12; t=ti:dt:te; pt=(10*t./1.5).*(t>=0 & t<=1.5)+(20-10*t./1.5).*(t>1.5 & t<=3)+0.*(t>3); m=length(pt); %load distribution in space rp=0.1; ri=0; rc=0.1; dr=0.001; r=ri:dr:rc; pr=1*(r<=rp)+0*(r>rp); n=length(pr); %load function with respect to t and r p=pr.'*pt; %green function G=1; cs=1; for i=1:1:m for j=1:1:n u(j,i)=heaviside(cs*t(i)-r(j))/(pi*G*sqrt(t(i)^2-(r(j)/cs)^2)); end end %convolution and response of displacement v=conv2(u,p); %plot the response history rr=r(n); [tpu,rpu]=meshgrid(ti:dt:2*te-dt,ri:dr:2*rc); [X,Y,Z]=meshgrid(linspace(min(tpu( ),max(tpu()),linspace(min(rpu(),max(rpu()),linspace(min(v(),max(v()));V=Y; h=contourslice(X,Y,Z,V,tpu,rpu,v,[0 0]+rr); set(h,'edgecolor','k'); contourslice(X,Y,Z,V,tpu,rpu,v,[0 0]+rr); xlabel('t/s'); ylabel('r/m'); zlabel('v/m'); axis([0 12 0 0.1 0 12e4]); view(0,0); grid on; title('displacement response history of point A or C'); |
回復(fù)此樓» 猜你喜歡
0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可
已經(jīng)有11人回復(fù)
-
南京大學(xué)化學(xué)調(diào)劑
已經(jīng)有15人回復(fù)
-
292求調(diào)劑
已經(jīng)有4人回復(fù)
-
材料專碩調(diào)劑
已經(jīng)有8人回復(fù)
-
085601材料工程找調(diào)劑
已經(jīng)有17人回復(fù)
-
面上5B能上會(huì)嗎？
已經(jīng)有6人回復(fù)
-
化學(xué)308分調(diào)劑
已經(jīng)有7人回復(fù)
-
化學(xué)工程085602 305分求調(diào)劑
已經(jīng)有25人回復(fù)
-
求化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
085602 化學(xué)工程專碩 340分求調(diào)劑
已經(jīng)有3人回復(fù)
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 hawaiicn 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085601材料工程找調(diào)劑 +16 | oatmealR 2026-03-29 | 17/850 |
|
|---|---|---|---|---|
|
|
[論文投稿]
arxiv背書人郵件求助
 20+3
|
redcom 2026-03-25 | 7/350 |
|
|
[考研] 一志愿北京工業(yè)大學(xué)，324分求調(diào)劑 +6 | 零八# 2026-03-28 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +15 | Micky11223 2026-03-25 | 20/1000 |
|
|
|
[考研] 299求調(diào)劑 +10 | 15188958825 2026-03-25 | 10/500 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|
|
|
[考研] 0856，材料與化工321分求調(diào)劑 +12 | 大饞小子 2026-03-27 | 13/650 |
|
|
|
[考研] 085602 307分 求調(diào)劑 +7 | 不知道叫什么！ 2026-03-26 | 7/350 |
|
|
|
[材料工程] 一志愿C9材料與化工專業(yè)總分300求調(diào)劑 +8 | 曼111 2026-03-24 | 9/450 |
|
|
|
[考研] 330一志愿中國(guó)海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
|
[考研] 0703化學(xué)338求調(diào)劑！ +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +6 | 林之夕 2026-03-24 | 6/300 |
|
|
|
[考研] 336材料求調(diào)劑 +7 | 陳瀅瑩 2026-03-26 | 9/450 |
|
|
|
[考研] 一志愿河工大 081700 276求調(diào)劑 +4 | 地球繞著太陽(yáng)轉(zhuǎn) 2026-03-23 | 4/200 |
|
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|