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

justlikethisis

function L2=lunwen(N,M)
. t: ^2 v4 I% }0 R8 `  E8 hlx=0;rx=1;lt=0;rt=1;
; G% u9 g3 m- `' G$ zh=(rx-lx)/N;
th=(rt-lt)/M;
x=0:h:1;
t=0:th:1;
A=zeros(N+1,N+1);
5 C2 P2 G* u1 W; k+ W1 y4 cB=zeros(N+1,N+1);
D=zeros(N+1,N+1);
E=zeros(2*N+2,2*N+2);
1 F9 G4 j! `( I+ Z. O' J; mF=zeros(N+1,1);
/ ~* l% T$ D: ]+ @u3=zeros(N+1,1);
7 l9 J+ }# R( g# v: F4 ~% EP=zeros(2*N+2,1);
X=zeros(2*N+2,1);
7 \' c  Z* i3 @X1=zeros(2*N+2,1);
( p5 C5 M' |! m+ d8 }/ W. g6 j0 [) RX2=zeros(2*N+2,M+1);
for i=0:N
    u3(i+1,1)=sin(pi*i*h);
5 Q* w1 C4 T) i7 Tend
X(N+2:2*N+2,1)=u3(1:N+1,1);
/ F+ W) V, J' g: \4 ]7 bfor i=1:N
" w8 k% Z. j, F3 Q) ]   e=[x(i),x(i+1)];
    A([i,i+1],[i,i+1])=A([i,i+1],[i,i+1])+ganddu1(e);
* v; d7 E% J4 _1 z    B([i,i+1],[i,i+1])=B([i,i+1],[i,i+1])+gangdu22(e);
" ^0 ^( h& t# m8 B- {* L    D([i,i+1],[i,i+1])=D([i,i+1],[i,i+1])+gangdu33(e);
end
1 u) V7 l/ P- m( F  E# |( m7 @% E(1:N+1,1:N+1)=A;
% E(1:N+1,N+2:2*N+2)=-B;
  f+ t( \/ [" z+ F% Y/ W0 k% 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];
2 @2 m: f, l1 t" ~+ w7 gfor j=1:M
    F=zeros(N+1,1); %%记得归零
    y=exp(-t(j+1));
: a7 Q. W: G, A8 M6 V    for i=1:N
        e=[x(i),x(i+1)];
& D) {: j1 n3 ?. z, I* c' M/ ~        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
    end
3 Q% a3 j) Z$ C- X# s4 z    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)];
: K+ q5 p( y. @: F/ a. s5 t# \+ uX2=[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)];
& w; R% _' J/ K, Y/ r7 b检查半天也不知道哪错了
; j1 \) J7 B4 p, l1 R
! R0 j" }$ R1 N. w6 B  ^
& c! k- w; n' D) {

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完