Question

我在运行程序时收到此消息：

metallror出错（第23行） A（端，端-1）=（2 * R）; ％byter utdenförstaelementeninnan - （1 /（2 * h））

我的A矩阵是否有问题？

该计划：

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

Answer 1

通过看起来似乎的外观，您想要为您的矩阵添加一些边框。这样做得恰到好处：

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世界，已经有一段时间了，因为我看到这么多被滥用的括号。： - ）

如何解决“尺寸不匹配”？

1 个答案: