版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

spvfly

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1131
帖子: 148
在線: 54.3小時(shí)
蟲號(hào): 812244
注冊: 2009-07-20
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

[求助] 一段用MATLAB編寫的代碼，結(jié)果錯(cuò)誤，請各位幫忙找找錯(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論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有299人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(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ù)

1樓 2016-12-25 21:47:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

丁丁糖27601

新蟲 (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 21.8
帖子: 30
在線: 4.9小時(shí)
蟲號(hào): 3984580
注冊: 2015-07-22
性別: GG
專業(yè): 空間物理

有限元？

發(fā)自小木蟲Android客戶端

回復(fù)此樓

2樓2016-12-25 22:10:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

spvfly

金蟲 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 1131
帖子: 148
在線: 54.3小時(shí)
蟲號(hào): 812244
注冊: 2009-07-20
性別: GG
專業(yè): 應(yīng)用地球物理學(xué)

引用回帖:

2樓: Originally posted by 丁丁糖27601 at 2016-12-25 22:10:38
有限元？

是滴

發(fā)自小木蟲Android客戶端

回復(fù)此樓

3樓2016-12-26 20:48:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

FMStation

至尊木蟲 (知名作家)

應(yīng)助: 591 (博士)
貴賓: 0.03
金幣: 18344.7
紅花: 97
帖子: 8891
在線: 1489.9小時(shí)
蟲號(hào): 2400059
注冊: 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

贊一下

回復(fù)此樓

4樓2016-12-28 20:38:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 spvfly 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 279分求調(diào)劑一志愿211 +9	chaojifeixia 2026-03-19	10/500	2026-03-20 12:29 by lature00
[考研] 317求調(diào)劑 +4	申子申申 2026-03-19	8/400	2026-03-20 11:20 by 申子申申
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +5	Ncdx123456 2026-03-19	5/250	2026-03-20 11:15 by wangy0907
[考研] 304求調(diào)劑 +5	曼殊2266 2026-03-18	5/250	2026-03-20 09:00 by ZHANG0tao
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考研] 081700化工學(xué)碩調(diào)劑 +3	【1】 2026-03-16	3/150	2026-03-19 23:40 by edmund7
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +4	生物工程調(diào)劑 2026-03-16	12/600	2026-03-19 11:10 by 生物工程調(diào)劑
[考研] 328求調(diào)劑，英語六級551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 268求調(diào)劑 +6	簡單點(diǎn)0 2026-03-17	6/300	2026-03-18 09:04 by 無際的草原
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 11408 一志愿西電，277分求調(diào)劑 +3	zhouzhen654 2026-03-16	3/150	2026-03-17 07:03 by laoshidan
[考研] 機(jī)械專碩325，尋找調(diào)劑院校 +3	y9999 2026-03-15	5/250	2026-03-16 19:58 by y9999
[基金申請] 今年的國基金是打分制嗎？ 50+3	zhanghaozhu 2026-03-14	3/150	2026-03-16 17:07 by 北京萊茵潤色
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6	太想進(jìn)步了0608 2026-03-16	6/300	2026-03-16 16:13 by kykm678
[考研] 0856專碩279求調(diào)劑 +5	加油加油！? 2026-03-15	5/250	2026-03-15 11:58 by 2020015
[考研] 22408總分284求調(diào)劑 +3	InAspic 2026-03-13	3/150	2026-03-15 11:10 by zhq0425
[考研] 一志愿哈工大材料324分求調(diào)劑 +5	閆旭東 2026-03-14	5/250	2026-03-14 14:53 by 木瓜膏
[考研] 297求調(diào)劑 +4	學(xué)海漂泊 2026-03-13	4/200	2026-03-14 11:51 by 熱情沙漠

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

spvfly

[求助] 一段用MATLAB編寫的代碼，結(jié)果錯(cuò)誤，請各位幫忙找找錯(cuò)誤出在哪了！�。� 已有1人參與

» 猜你喜歡

丁丁糖27601

spvfly

FMStation

【答案】應(yīng)助回帖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

[求助] 一段用MATLAB編寫的代碼，結(jié)果錯(cuò)誤，請各位幫忙找找錯(cuò)誤出在哪了！�。� 已有1人參與