| 6 | 1/1 | 返回列表 |
| 查看: 626 | 回復: 5 | ||
阿毛蟲銀蟲 (正式寫手)
|
[求助]
蟲友們,看一下這個程序,運行一下 已有1人參與
|
| 這個程序我就想將輸出的四個系數(shù)值,全部編輯到coefficient的矩陣中,然后一起輸出來,不想一次輸出一個來,這樣的話,后面緊接著畫圖,感覺不能直接畫在一幅圖上,還有就是為什么這個程序最后會輸出個ans=。。。。,求助謝謝 |

銀蟲 (正式寫手)
| 不好意思,附件沒傳上 |

| 改了一下,你可以看看 |
銀蟲 (正式寫手)

|
function [orthogonality, coefficient_all]=clamped_coefficient(c) format long clear all G=8.0e10; K=1.2; E=2.06e11; P=7.9e3; l=4; D=0.6; H=0.2; k=1/K; I=D*(H^3)/12; A=D*H; figure; coefficient_all = []; r=sqrt(I/A); s=sqrt(E*I/(k*A*G)); syms x ; syms c; frequency=[405.35308 1093.84137 2087.15447 3341.29751]; for ii=1:length(frequency) b=sqrt(P*A/(E*I))*frequency(ii); alfa=b*sqrt(-(r^2+s^2)+sqrt((r^2-s^2)^2+4./(b^2)))/sqrt(2); beta=b*sqrt(r^2+s^2+sqrt((r^2-s^2)^2+4./(b^2)))/sqrt(2); y=cosh(alfa*x)+(alfa*(beta^2-b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sinh(alfa*x)..., -cos(beta*x)-(beta*(alfa^2+b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sin(beta*x); y_int=int(y^2,0,l); fai=(alfa^2+b^2*s^2)/alfa*sinh(alfa*x)+(alfa*(beta^2-b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2))*(alfa^2+b^2*s^2)/alfa*cosh(alfa*x)..., -(b^2*s^2-beta^2)/beta*sin(beta*x)-(beta*(alfa^2+b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2))*(b^2*s^2-beta^2)/beta*cos(beta*x); fai_int=int(fai^2,0,l); orthogonality=c^2*y_int+c^2*r^2*fai_int-1/(P*A); coefficient_result=solve(orthogonality,c) ; coefficient_results=vpa(coefficient_result,6); coefficient = abs(coefficient_results(1)); coefficient_all = [coefficient_all; coefficient]; x_shuzhi=0:0.01:4; y_shuzhi=cosh(alfa.*x_shuzhi)+(alfa*(beta^2-b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sinh(alfa.*x_shuzhi)..., -cos(beta.*x_shuzhi)-(beta*(alfa^2+b^2*s^2)*(cos(beta*l)-cosh(alfa*l)))/((beta^2-b^2*s^2)*alfa*sinh(alfa*l)-beta*sin(beta*l)*(alfa^2+b^2*s^2)).*sin(beta.*x_shuzhi); hold on; h(ii) = plot(x_shuzhi,y_shuzhi); end |
銀蟲 (正式寫手)

| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 081700化工學碩調劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 0703化學調劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 0817 化學工程 299分求調劑 有科研經(jīng)歷 有二區(qū)文章 +18 | rare12345 2026-03-18 | 18/900 |
|
|
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +8 | zlingli 2026-03-13 | 8/400 |
|
|
[考研] 085601材料工程專碩求調劑 +10 | 慕寒mio 2026-03-16 | 10/500 |
|
|
[考研] 一志愿天津大學化學工藝專業(yè)(081702)315分求調劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|
[考研] 304求調劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
[考研] 324分 085600材料化工求調劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[教師之家] 焦慮 +9 | 水冰月月野兔 2026-03-13 | 13/650 |
|
|
[考研] 085410人工智能專碩317求調劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 0854可跨調劑,一作一項核心論文五項專利,省、國級證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考博] 環(huán)境領域全國重點實驗室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[考研] 考研化學學碩調劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 0856專碩279求調劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 本科南京大學一志愿川大藥學327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 266求調劑 +4 | 學員97LZgn 2026-03-13 | 4/200 |
|