 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 612 | 回復(fù): 3 | ||
spvfly金蟲 (小有名氣)
|
[求助]
一段用MATLAB編寫的代碼，結(jié)果錯(cuò)誤，請(qǐng)各位幫忙找找錯(cuò)誤出在哪了！�。� 已有1人參與
|
|
%%%主程序 XY=xlsread('data.xlsx','ND'); I3=xlsread('data.xlsx','NE'); UB=xlsread('data.xlsx','UB'); ND=33;%點(diǎn)數(shù) NE=40;%單元數(shù) ND1=24;%邊值點(diǎn)數(shù) UK=UK1(ND,NE,I3,XY);%總系數(shù)矩陣 [UK,UB1]=UBK(UK,UB,ND1);%將KU=0方程組變?yōu)镵U=B形式的方程組 U=UK\UB1'%計(jì)算得到結(jié)果 %%%%%%%% 總系數(shù)矩陣程序 %%%%%%%% function y=UK1(ND,NE,I3,XY) X=[0,0,0]; Y=[0,0,0]; SK=zeros(ND,ND); KE=zeros(3,3); for P=1:NE %遍歷各個(gè)單元 for I=1:3 %遍歷一個(gè)單元的三個(gè)節(jié)點(diǎn) L=I3(I,P); %確定一個(gè)單元的節(jié)點(diǎn) X(I)=XY(1,L); %該節(jié)點(diǎn)的x坐標(biāo) Y(I)=XY(2,L); %該節(jié)點(diǎn)的y左邊 end KE=UKE1(X,Y);%一個(gè)單元的單元系數(shù)矩陣 %下面這段程序是為了把單元系數(shù)矩陣放入定帶寬矩陣中去 for J=1:3 for K=1:3 NJ=I3(J,P); NK=I3(K,P); SK(NJ,NK)=SK(NJ,NK)+KE(J,K); end end end y=SK; %%%%%%% 單元系數(shù)矩陣 %%%%%%% function y=UKE1(X,Y) K1=zeros(3,3);%定義一個(gè)3*3數(shù)組用于存放單元系數(shù)矩陣的值 A=[0,0,0]; B=[0,0,0]; C=[0,0,0]; A(1)=Y(2)-Y(3); A(2)=Y(3)-Y(1); A(3)=Y(1)-Y(2); B(1)=X(3)-X(2); B(2)=X(1)-X(3); B(3)=X(2)-X(1); S=2*(A(1)*B(2)-A(2)*B(1)); for I=1:3 for J=1:3 K1(I,J)=(A(I)*A(J)+B(I)*B(J))/S; end end y=K1; %%%%%%%%% 將KU=0形式轉(zhuǎn)換為KU=B形式 %%%%%%%%% function [X,Y]=UBK(UK,UB,ND1) UB1=zeros(1,33); for I=1:ND1 D=UB(1,I); V=UB(2,I); UK(D,D)=UK(D,D)*1e10; UB1(D)=V*UK(D,D); end X=UK; Y=UB1; %%%%%%% 所得到的錯(cuò)誤結(jié)果 %%%%%% 程序結(jié)果 解析解 0.354616826 0.302042805 0.455989065 0.355568018 0.562915374 0.41747372 0.675049586 0.482845483 0.767370579 0.526986487 0.761868629 0.517705616 0.621009799 0.455065182 0.46388813 0.373490832 0.342221541 0.301003461 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有6人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有212人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請(qǐng)問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學(xué)工程學(xué)碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有65人回復(fù)
-
華東師范大學(xué)芯片設(shè)計(jì)徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
丁丁糖27601
新蟲 (初入文壇)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 21.8
- 帖子: 30
- 在線: 4.9小時(shí)
- 蟲號(hào): 3984580
- 注冊(cè): 2015-07-22
- 性別: GG
- 專業(yè): 空間物理
2樓2016-12-25 22:10:38
spvfly
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1131
- 帖子: 148
- 在線: 54.3小時(shí)
- 蟲號(hào): 812244
- 注冊(cè): 2009-07-20
- 性別: GG
- 專業(yè): 應(yīng)用地球物理學(xué)
3樓2016-12-26 20:48:25
FMStation
至尊木蟲 (知名作家)
- 應(yīng)助: 591 (博士)
- 貴賓: 0.03
- 金幣: 18344.7
- 紅花: 97
- 帖子: 8891
- 在線: 1489.9小時(shí)
- 蟲號(hào): 2400059
- 注冊(cè): 2013-04-04
- 專業(yè): 計(jì)算機(jī)應(yīng)用技術(shù)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
https://www.mathworks.com/help/matlab/debugging-code.html https://www.mathworks.com/help/m ... s-and-features.html |
4樓2016-12-28 20:38:37
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 spvfly 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)調(diào)劑 ，六級(jí)已過(guò)，有科研經(jīng)歷 +13 | 曦熙兮 2026-03-15 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +21 | rare12345 2026-03-18 | 21/1050 |
|
|
|
[考研] 一志愿中南化學(xué)337求調(diào)劑 +4 | niko- 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 材料考研調(diào)劑 +3 | xwt。 2026-03-19 | 3/150 |
|
|
|
[考研] 332求調(diào)劑 +3 | ydfyh 2026-03-17 | 3/150 |
|
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +5 | hanamiko 2026-03-18 | 5/250 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|