24小时热门版块排行榜    

CyRhmU.jpeg
查看: 2043  |  回复: 1

zxiaofeizjnu

新虫 (初入文坛)

[求助] matlab中大的数值矩阵行列式等于0,如何处理

大侠:
  请问matlab中需要编写大的数值矩阵,且最后要求行列式等于0,这样的要求可以达到么?
下面是我编写的程序:
i=1:n;
  W(i)=(m+0.5)./(i.*a)+pi.*B.*i.*a/L;
  Z1=[w,0;-1+W(1),0];
  Z2=[0,-1;0,0];
  for j=1:n-1;
    Y=[-1,W(j)+1;0,w];
    X=[0,w;1,-1+W(j+1)];
    K(2.*j-1:2.*j,n+j:n+j+1)=  K(2.*j-1:2.*j,n+j:n+j+1)+Y;
    K(2.*j+1:2.*j+2,j:j+1)= K(2.*j+1:2.*j+2,j:j+1)+X;
  end

    K(1:2,1:2)= K(1:2,1:2)+Z1;

    K(n-1:n,n-1:n)= K(n-1:n,n-1:n)+Z2;

    y=abs(det(K))
    if  (y=0)
        z(t)=E;
        v(t)=B;
   
     t=t+1;
    end

   end
end
plot(v,z)
问下错那里了,如何修改呢?谢谢
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

我ai圈圈糖

银虫 (小有名气)

...............我想问 你怎么解决的,我现在也遇到相关问题了
2楼2018-12-07 15:03:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zxiaofeizjnu 的主题更新
信息提示
请填处理意见