| 24小時熱門版塊排行榜 |
| 5 | 1/1 | 返回列表 |
| 查看: 618 | 回復(fù): 3 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
spvfly金蟲 (小有名氣)
|
[求助]
一段用MATLAB編寫的代碼,結(jié)果錯誤,請各位幫忙找找錯誤出在哪了。。 已有1人參與
|
|
|
%%%主程序 XY=xlsread('data.xlsx','ND'); I3=xlsread('data.xlsx','NE'); UB=xlsread('data.xlsx','UB'); ND=33;%點數(shù) NE=40;%單元數(shù) ND1=24;%邊值點數(shù) UK=UK1(ND,NE,I3,XY);%總系數(shù)矩陣 [UK,UB1]=UBK(UK,UB,ND1);%將KU=0方程組變?yōu)镵U=B形式的方程組 U=UK\UB1'%計算得到結(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 %遍歷各個單元 for I=1:3 %遍歷一個單元的三個節(jié)點 L=I3(I,P); %確定一個單元的節(jié)點 X(I)=XY(1,L); %該節(jié)點的x坐標 Y(I)=XY(2,L); %該節(jié)點的y左邊 end KE=UKE1(X,Y);%一個單元的單元系數(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);%定義一個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; %%%%%%% 所得到的錯誤結(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 |
金蟲 (小有名氣)
新蟲 (初入文壇)
至尊木蟲 (知名作家)
|
https://www.mathworks.com/help/matlab/debugging-code.html https://www.mathworks.com/help/m ... s-and-features.html |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +7 | 十三加油 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 081700 調(diào)劑 267分 +7 | 迷人的哈哈 2026-03-23 | 7/350 |
|
|
[考研] 材料調(diào)劑 +3 | 匹克i 2026-03-23 | 3/150 |
|
|
[考研] 291求調(diào)劑 +8 | hhhhxn.. 2026-03-23 | 8/400 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 求老師收我 +3 | zzh16938784 2026-03-23 | 3/150 |
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
[考研] 石河子大學(xué)(211、雙一流)碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
[考研] 華東師范大學(xué)-071000生物學(xué)-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 一志愿華中科技大學(xué),080502,354分求調(diào)劑 +5 | 守候夕陽CF 2026-03-18 | 5/250 |
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
[考研] A區(qū)線材料學(xué)調(diào)劑 +5 | 周周無極 2026-03-20 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 261求B區(qū)調(diào)劑,科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|