UC知道在线急求matlab解方程
来源：百度知道 编辑：UC知道 时间：2024/05/22 00:28:30
syms w;
md=1700
g=8.27*10^9
n=10
l=10
%给矩阵赋值
for i=1:n-1
for j=1:n-1
if (j==i-1)
s(i,j)=n^2/l^2;
elseif (j==i)
s(i,j)=md*(w^2)/g-(2*i-1)*(n^2)/i*l^2;
elseif (j==i+1)
s(i,j)=((n/l)^2)/i;
end
end
end
s(n-1,n-1)=md*w^2/g-(n/l)^2;
s
%%解行列式为0的方程，y即未知数w的值
y=det(s)
其中的w是未知数，也就是要求的，只就正w
麻烦大侠把上面的程序编完整
谢谢
能不能z给的详细些，我没有学过matlab，现在急着用，没有时间学了，谢谢！！！
md=1700
g=8.27*10^9
n=10
l=10
%给矩阵赋值
for i=1:n-1
for j=1:n-1
if (j==i-1)
s(i,j)=n^2/l^2;
elseif (j==i)
s(i,j)=md*(w^2)/g-(2*i-1)*(n^2)/i*l^2;
elseif (j==i+1)
s(i,j)=((n/l)^2)/i;
end
end
end
s(n-1,n-1)=md*w^2/g-(n/l)^2;
s
%%解行列式为0的方程，y即未知数w的值
y=det(s)
其中的w是未知数，也就是要求的，只就正w
麻烦大侠把上面的程序编完整
谢谢
能不能z给的详细些，我没有学过matlab，现在急着用，没有时间学了，谢谢！！！
325665
没时间仔细做，给你指出几个问题。
s(i,j)应该用function函数定义一下吧，你既然说是矩阵了，s(i,j)=n^2/l^2;应改为s(i,j)=n.^2/l.^2;
s(i,j)=md*(w^2)/g-(2*i-1)*(n^2)/i*l^2;应改为s(i,j)=md.*(w.^2)/g-(2*i-1)*(n.^2)/i*l.^2;
下面的几处也要改的………………
呵呵，多多看看书。