Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/matlab/15.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
Matlab 如何解决;维度不匹配;?_Matlab_Matrix_Dimension_Mismatch - Fatal编程技术网

Matlab 如何解决;维度不匹配;?

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

我在运行程序时收到以下消息:

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;               %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所说,这里不是每个人都说瑞典语,英语网站也是如此。提供英文注释对您的事业是有益的,因为这样每个人(不仅仅是我和少数讲瑞典语的人)都能理解您的代码。