| 查看: 617 | 回复: 3 | ||
[求助]
一段用MATLAB编写的代码，结果错误，请各位帮忙找找错误出在哪了！！！ 已有1人参与
|
|
%%%主程序 XY=xlsread('data.xlsx','ND'); I3=xlsread('data.xlsx','NE'); UB=xlsread('data.xlsx','UB'); ND=33;%点数 NE=40;%单元数 ND1=24;%边值点数 UK=UK1(ND,NE,I3,XY);%总系数矩阵 [UK,UB1]=UBK(UK,UB,ND1);%将KU=0方程组变为KU=B形式的方程组 U=UK\UB1'%计算得到结果 %%%%%%%% 总系数矩阵程序 %%%%%%%% 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 %遍历一个单元的三个节点 L=I3(I,P); %确定一个单元的节点 X(I)=XY(1,L); %该节点的x坐标 Y(I)=XY(2,L); %该节点的y左边 end KE=UKE1(X,Y);%一个单元的单元系数矩阵 %下面这段程序是为了把单元系数矩阵放入定带宽矩阵中去 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; %%%%%%% 单元系数矩阵 %%%%%%% function y=UKE1(X,Y) K1=zeros(3,3);%定义一个3*3数组用于存放单元系数矩阵的值 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形式转换为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; %%%%%%% 所得到的错误结果 %%%%%% 程序结果 解析解 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 |
回复此楼» 猜你喜欢
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有32人回复
-
物理学I论文润色/翻译怎么收费?
已经有67人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有78人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
2楼2016-12-25 22:10:38
3楼2016-12-26 20:48:25
FMStation
至尊木虫 (知名作家)
- 应助: 591 (博士)
- 贵宾: 0.03
- 金币: 18349.2
- 红花: 97
- 帖子: 8891
- 在线: 1490.1小时
- 虫号: 2400059
- 注册: 2013-04-04
- 专业: 计算机应用技术
4楼2016-12-28 20:38:37
