鐑�浼犲�兼柟绋嬪悜鍚庡樊鍒嗘牸寮忕殑MATLAB绋嬪簭 - 鐧惧害鏂囧簱

鏂板缓

涓婁紶

棣栭〉

鍔╂墜

鏈�杩?/div>

璧勬枡搴?/div>

宸ュ叿

鍚戝悗宸�鍒嗘牸寮�

MATLAB

缂栫▼锛?/p>

 

c

lear;clc; 

format short e 

a=input('

璇疯緭鍏ョ郴鏁?/p>

a

鐨勫�?/p>

'); 

l=input('

璇疯緭鍏ラ暱搴?/p>

l

鐨勫�?/p>

'); 

M=input('

璇疯緭鍏ュ皢鍖洪棿

[0,1]

绛夊垎鐨勪釜鏁?/p>

M '); 

ot=input('

璇疯緭鍏ユ椂闂村�為�?/p>

ot

鐨勫�?/p>

'); 

n=input('

璇疯緭鍏ヨ繍琛屾�℃�?/p>

n

鐨勫�?/p>

'); 

ox=1/M; 

x0=zeros(M+1,1) 

for ii=1:M 

 

 

 

 

x0(ii+1)=ii*ox; 

end 

u=sin(pi*x0/l); 

r=a*ot/(ox)^2; 

for ii=1:n 

 

 

 

 

%

鏁版嵁鐨勮緭鍏?/p>

 

 

 

 

 

B=zeros(M-1,1); 

 

 

 

 

A=zeros(M-2,1); 

 

 

 

 

C=zeros(M-2,1); 

 

 

 

 

S=zeros(M-1,1); 

 

 

 

 

for ii=1:M-2 

 

 

 

 

 

 

 

 

B(ii)=1+2*r;A(ii)=-r;C(ii)=-r; 

 

 

 

 

 

 

 

 

S(ii)=u(ii+1,1); 

 

 

 

 

end 

 

 

 

 

B(M-1,1)=1+2*r;S(M-1,1)=u(M,1);u(1,2)=0;u(M+1,2)=0; 

 

 

 

 

S(1,1)=S(1,1)+r*u(1,2);S(M-1,1)=S(M-1,1)+r*u(M+1,2); 

 

 

 

 

%

杩借刀娉?/p>

 

 

 

 

 

S(1)=S(1)/B(1);T=B(1);k=2; 

 

 

 

 

while k~=M 

 

 

 

 

 

 

 

 

B(k-1)=C(k-1)/T; 

 

 

 

 

 

 

 

 

T=B(k)-A(k-1)*B(k-1); 

 

 

 

 

 

 

 

 

S(k)=(S(k)-A(k-1)*S(k-1))/T; 

 

 

 

 

 

 

 

 

k=k+1 

 

 

 

 

end 

 

 

 

 

k=1; 

 

 

 

 

while k~=M-1 

 

 

 

 

 

 

 

 

S(M-1-k)=S(M-1-k)-B(M-1-k)*S(M-k); 

 

 

 

 

 

 

 

 

k=k+1; 

 

 

 

 

end 

 

 

 

 

u(2:M,2)=S; 

 

 

 

 

u(:,1)=u(:,2); 

end 

%

璁＄畻绮剧‘瑙?/p>

 

for x=0:M