主题:复数矩阵满足这样条件该怎么办
tianhy2010
[专家分:60] 发布于 2010-10-20 15:46:00
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
最后更新于:2010-10-20 15:47:00
回复列表 (共4个回复)
沙发
jstzhurj [专家分:4680] 发布于 2010-10-20 16:23:00
复数能比较大小吗?应该了解上面式子的涵义。
板凳
tianhy2010 [专家分:60] 发布于 2010-10-23 22:53:00
[quote]
复数能比较大小吗?应该了解上面式子的涵义。[/quote]
[em2][em2][em2][em8][em8][em8]
如果t是复矩阵,维数为(m,m),eps=1e-5,有没有对t内每个矩阵元求模,然后再对所有的模求和,然后判断这个和是否小于eps,有这样的函数吗?求模的函数
3 楼
dongyuanxun [专家分:7180] 发布于 2010-10-23 23:27:00
[quote][quote]
复数能比较大小吗?应该了解上面式子的涵义。[/quote]
[em2][em2][em2][em8][em8][em8]
如果t是复矩阵,维数为(m,m),eps=1e-5,有没有对t内每个矩阵元求模,然后再对所有的模求和,然后判断这个和是否小于eps,有这样的函数吗?求模的函数[/quote]
活用abs、矩阵操作和sum
4 楼
tianhy2010 [专家分:60] 发布于 2010-10-24 16:56:00
[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]
我来回复