一维抛物有限元编程，检查半天也不知道哪错了，大家可以帮忙看看吗，谢谢了

justlikethisis · 发表于 2020-9-11 14:12

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

function L2=lunwen(N,M)
lx=0;rx=1;lt=0;rt=1;
h=(rx-lx)/N;
th=(rt-lt)/M;
x=0:h:1;
t=0:th:1;
A=zeros(N+1,N+1);
B=zeros(N+1,N+1);
D=zeros(N+1,N+1);
E=zeros(2*N+2,2*N+2);
F=zeros(N+1,1);
u3=zeros(N+1,1);
P=zeros(2*N+2,1);
X=zeros(2*N+2,1);
X1=zeros(2*N+2,1);
X2=zeros(2*N+2,M+1);
for i=0:N
u3(i+1,1)=sin(pi*i*h);
end
X(N+2:2*N+2,1)=u3(1:N+1,1);
for i=1:N
e=[x(i),x(i+1)];
A([i,i+1],[i,i+1])=A([i,i+1],[i,i+1])+ganddu1(e);
B([i,i+1],[i,i+1])=B([i,i+1],[i,i+1])+gangdu22(e);
D([i,i+1],[i,i+1])=D([i,i+1],[i,i+1])+gangdu33(e);
end
% E(1:N+1,1:N+1)=A;
% E(1:N+1,N+2:2*N+2)=-B;
% E(N+2:2*N+2,1:N+1)=-th*B;
% E(N+2:2*N+2,N+2:2*N+2)=B-th*B+th*D;
E=[A,-B;-th.*B,B-th.*B+th.*D];
for j=1:M
F=zeros(N+1,1); %%记得归零
y=exp(-t(j+1));
for i=1:N
e=[x(i),x(i+1)];
F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
end
X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
end
X1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
X2=[X1;X];
错误使用 vertcat
要串联的数组的维度不一致。

出错 lunwen (line 39)
X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
检查半天也不知道哪错了

qq666888qqw · 发表于 2020-9-11 15:11

你的程序没传完吧

somethingabc · 发表于 2020-9-11 16:35

程序好像没传完

帐号		自动登录	找回密码
密码			注册

一维抛物有限元编程，检查半天也不知道哪错了，大家可以帮忙看看吗，谢谢了

EDA365欢迎您登录！

浏览过的版块

推荐内容 /1