回 帖 发 新 帖 刷新版面

主题:复数矩阵满足这样条件该怎么办

t(n+1),t(n+1)~都是复数矩阵,大小(m,m),s是个小量,比如s=1e-6,当
t(n+1),t(n+1)~都小于s时,执行某个操作。这个条件应该实部和虚报都满足小于s吗?

比如
if(all(real(t(n+1))<eps) .and. all(imag(t(n+1))<eps)) exit

回复列表 (共4个回复)

沙发


复数能比较大小吗?应该了解上面式子的涵义。

板凳

[quote]
复数能比较大小吗?应该了解上面式子的涵义。[/quote]
[em2][em2][em2][em8][em8][em8]

如果t是复矩阵,维数为(m,m),eps=1e-5,有没有对t内每个矩阵元求模,然后再对所有的模求和,然后判断这个和是否小于eps,有这样的函数吗?求模的函数

3 楼

[quote][quote]
复数能比较大小吗?应该了解上面式子的涵义。[/quote]
[em2][em2][em2][em8][em8][em8]

如果t是复矩阵,维数为(m,m),eps=1e-5,有没有对t内每个矩阵元求模,然后再对所有的模求和,然后判断这个和是否小于eps,有这样的函数吗?求模的函数[/quote]

活用abs、矩阵操作和sum

4 楼

[quote]
复数能比较大小吗?应该了解上面式子的涵义。[/quote]
紧接上次求助你的问题,定义type类型那个,
http://bbs.pfan.cn/post-350875.html
如果t(i).a表示一个矩阵,那么可以用t(i).a(j,k)表示这个矩阵里的元素吗?
eps=1.0e-3
pp=0.0
do j=1,m
  do k=1,m
  pp=pp+abs(t(i).a(j,k))
  if(pp<eps) exit
  end do
end do
[em2][em2][em2][em2][em2][em2][em2][em2][em2][em2][em2][em2]

我来回复

您尚未登录,请登录后再回复。点此登录或注册