Matlab 如何解决;维度不匹配;?
我在运行程序时收到以下消息: MetalRor中的错误(第23行) A(结束,结束-1)=(2.*r);%byter ut den första elementen innan-(1/(2*h)) 我的Matlab 如何解决;维度不匹配;?,matlab,matrix,dimension,mismatch,Matlab,Matrix,Dimension,Mismatch,我在运行程序时收到以下消息: MetalRor中的错误(第23行) A(结束,结束-1)=(2.*r);%byter ut den första elementen innan-(1/(2*h)) 我的矩阵有什么问题吗 该方案: clc,clf,clear all %Givet Ti=450; r0=1;R=2;K=1;Te=20; x=[1 2 4 8 ]; %används för att fördubbla våra intervall g=0; %fo
矩阵有什么问题吗
该方案:
clc,clf,clear all
%Givet
Ti=450; r0=1;R=2;K=1;Te=20;
x=[1 2 4 8 ]; %används för att fördubbla våra intervall
g=0; %fortsätt läs koden!
% N=4; %Antal intervall
% n=N;
for i=1:4
N=25*x(i); %fördubblar delintervallena
h=(R-r0)/N; %steglängd
r=r0+h:h:R; %Inre punkter plus sista obekanta
sub=(2*r-h)'; %'sub diagonal
dia=(-4.*r)'; %'huvuddiagonal
sup=(2*r+h)'; %'super diagonal
%vänsterled------
A=spdiags([sub dia sup],-1:1,N,N); % bildar en gles matris
size(A)
A(end,end-1)=(2.*r); % byter ut den första elementen innan -(1/(2*h))
% med vårt handberäknade värde
A(end,end)=((-2.*r)-(2.*r.*h.*K)-((h.^2).*K)); % byter ut den sista elementen innan K
%högerled
B=zeros(N,1); % högerled
B(1,1)=-Ti*((2*r0)-h); % byter ut elementen med
% vårt handberäknade värde
B(end,1)=((-2*h*K*r)*Te-((h.^2)*K)*Te); % byter ut elementen innan KTe
%Beräknar
U=A\B
u=[Ti U']; %'
%visar resultatet
disp(['temperaturen vid r=2 med ',num2str(N),' delintervall är: ',num2str(u(end))])
%Ritar vår graf
g=g+1; %används för att rita olika diagram
subplot(2,2,g)
size(r)
size(u)
plot([1 r],u),grid on %[1 r],u)
xlabel('Radie')
ylabel('Temperatur')
title(['Diagram med ',num2str(N),' delintervall' ])
end
从你试图做的事情来看,你想给你的矩阵添加一些边框。这样做是正确的:
A(:,end-1) = 2*r;
A(:,end) = - 2*r - 2*h*K*r - h*h*K;
%// some more code
B(:,1) = (- 2*h*K*r - h*h*K)*Te;
而且,如果你们允许我开个玩笑的话:你们一定来自Lisp世界,我已经很久没有看到这么多被滥用的括号了。:-) 请至少努力提供英文评论。欢迎使用StackOverflow:)您正在将数组(r
)分配给a
中的单个矩阵项。那当然不合适。这就像把你所有的20辆车都停在一个车库里。正如JaBe所说,这里不是每个人都说瑞典语,英语网站也是如此。提供英文注释对您的事业是有益的,因为这样每个人(不仅仅是我和少数讲瑞典语的人)都能理解您的代码。