I'm getting a matrix as output when I'm expecting a vector

Question

Arthur Bergé 2024-8-6，9:57

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2143351-i-m-getting-a-matrix-as-output-when-i-m-expecting-a-vector

编辑： Torsten 2024-8-6，10:43

I'm trying to calculate the abscisses and weight values for the Newton-Cotes method. In the end I use a linsolve() expecting to get a vector as output but instead I get a matrix. Any ideas what I could have done wrong?

a = 0;
b = 1;
n = 20;
[x,H] = NewtonCotes(a,b,n)
Warning: Matrix is close to singular or badly scaled. Results may be inaccurate. RCOND =  5.742547e-17.
x = 1x21
         0    0.0500    0.1000    0.1500    0.2000    0.2500    0.3000    0.3500    0.4000    0.4500    0.5000    0.5500    0.6000    0.6500    0.7000    0.7500    0.8000    0.8500    0.9000    0.9500    1.0000
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
H = 21x21
1.0e+12 *

    7.0832   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026   -0.9026
    0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000
    5.6187   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159   -0.7159
    0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000
    3.0530   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890   -0.3890
    0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000   -0.0000
    0.0612   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078   -0.0078
    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000
   -2.5218    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213    0.3213
    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000    0.0000
<mw-icon class=""></mw-icon>
<mw-icon class=""></mw-icon>
function [x,H] = NewtonCotes(a,b,n)
if n ~= 0
    h = (b-a)/n;
else
    h = (b-a);
end
for i = 1:n+1
    x(i) = a + ((i-1)*h);
end
H = ComputeWeights(a,b,x);
end
function  [H] = ComputeWeights(a,b,x)
n = length(x)-1;
I = ones(n+1);
I(1) = b - a;
I(2) = 0;
for i = 3:n+1
    I(i) = ((b-a)*((-1)^i + 1))/(2*(1-i^2));
end
for i = 1:length(x)
    x(i) = ((2/(b-a))*(x(i)-a))-1;
end
T = zeros(n+1,n+1);
for i = 1:n+1
    for k = 1:n+1
        T(i,k) = cos(k*acos(x(i)));
    end
end
H = linsolve(T,I);
end

2 个评论
显示无隐藏无

Rahul 2024-8-6，10:21

Hi @Arthur Bergé,

What inputs are you passing to your function?

I have tried to use some dummy inputs using your function and it was returning a vector instead of a matrix.

Torsten 2024-8-6，10:27

I is an (n+1)x(n+1) matrix in the code. So H will always be a (n+1)x(n+1) matrix, not a vector.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Torsten 2024-8-6，10:19

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2143351-i-m-getting-a-matrix-as-output-when-i-m-expecting-a-vector#answer_1495291

在 MATLAB Online 中打开

Maybe you mean

I = ones(n+1,1);

instead of

I = ones(n+1);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

Shivansh 2024-8-6，10:14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2143351-i-m-getting-a-matrix-as-output-when-i-m-expecting-a-vector#answer_1495286

编辑：Shivansh 2024-8-6，10:15

Hello Arthur,

You can get a matrix as output while using linsolve to solve the linear systems of equation.

The function returns a vector or matrix that satisfies AX=B. The size of X depends on whether "opts.TRANSA" = true:

If A is m-by-n and B is m-by-k, then X is n-by-k and is the solution to AX = B.
If "opts.TRANSA" = true, then A is m-by-n and B is m-by-k. In this case, X is m-by-k and is the solution to A'X = B.

Here, "TRANSA" is Conjugate transpose which specifies whether the function solves A*X = B (when opts.TRANSA = false) or the transposed problem A'*X = B (when opts.TRANSA = true).

Please refer to the following documentation link for more information:

Output type for Linsolve: https://www.mathworks.com/help/matlab/ref/linsolve.html?searchHighlight=linsolve%20output%20matrix&s_tid=srchtitle_support_results_1_linsolve%20output%20matrix#mw_7666b08b-4d6b-4175-8f3b-c0c979bd9c6d.

I hope it resolves your query.