Unexpected complex coefficients in a matrix

2 次查看（过去 30 天）
Enterprixe
Enterprixe 2017-3-30
回答： David Goodmanson 2017-3-30
Hello, im trying to create a panel method for the study of an airfoil. Inside this problem, i have to obtain a matrix by combination of points created which define the panel distributions and their middle points. So the code is:
for n=1:201
x(n)= 1/2-1/2*cos((n-1)*pi/200);
end
for n=1:200
Pc(n)=(x(n)+x(n+1))/2;
end
for i=1:200
for j=1:200
B(i,j)=(-((Pc(i)-x(j+1))/(x(j+1)-x(j)))*log((Pc(i)-x(j+1))/(Pc(i)-x(j)))-1);
end
end
So, in theory, the matrix B should be not complex at all, with its main diagonal full of just '-1' in each tearm, and that main diagonal should be dominant compared to the rest of coefficients of the matrix.

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

采纳的回答

David Goodmanson
David Goodmanson 2017-3-30
Hello Javier, The only opportunity here to get a complex number is if the argument of the log term is negative. For all the nondiagonal terms of B you are all right, but for i=j the code produces log(-1) and you get an extra i pi/2 on the diagonal.

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Airfoil tools 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by