24小时热门版块排行榜    

查看: 1368  |  回复: 3
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

sololee2010

木虫 (正式写手)


[交流] 虫友问的一个关于压缩感知问题的解答:如何将Y=A*S等价变形为Y(:)=B*S(:)

.............貌似在有些应用场合中会碰到如下的数学模型:
                                               Y = A*S ,                  (1)
其中Y, A, S均为矩阵。提出的问题是:已知矩阵Y和矩阵A,求解矩阵S,且已知S是一个稀疏矩阵。
       一看到稀疏,可能马上有人想起是不是可以利用压缩感知CS重构算法寻找稀疏矩阵S呢?但是又发现不对,CS的数学模型和(1)式不一样, CS的数学模型是:
                                               y = Phi * x,                 (2)
其中y和x是列向量,Phi是矩阵。如何能够将(1)式等价变形为(2)式,按照我的理解,就是这位虫友提出的问题。
        实际上,当然是可以的,只是个简单的线性代数的基本问题,以2维矩阵为例,我们看下如何变形,见附件的例图(图片上传不了,没办法,只能以附件的形式了)。这样就可以找到Y=A*S的等价变形:
                                                                                   Y(: ) = B * S(: ),            (3)
当然S(: )是稀疏的,那么就有可能利用CS重构算法求解S,但需要满足的条件就是B具有受限等距性(RIP)。
        如果这位虫友遇到的是低维(代数维)问题,就是S(: )的长度比较短,那么矩阵B就可以在计算机直接构造出来,但是对于高维问题,B是无法构造的,那么实际就应该找到(3)式的等价算子,即输入是S(: ),输出是Y(: )的函数。目前大多数流行的CS重构算法都是支持算子的,如果不支持算法的CS重构算法实际上就不是成熟的算法,因为该算法无法处理高维信号(如图像信号)。
        保证CS重构算法收敛的条件是B满足RIP性,那么接下来的任务就是去数学证明B是否满足RIP,这个像我这种工程出身的虫子,就非常难了。可以先随便弄个什么CS算法测试一下,看看可不可以有效复原S(: ),如果效果很好,就应该可以证明方法是可行的,再来思考B的RIP性,一旦可以用数学证明出来,那么这应该就是个不错的Idear了。[ Last edited by sololee2010 on 2013-5-31 at 11:15 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : IMG_20130531_100845.jpg
  • 2013-05-31 11:02:08, 562.26 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

taoyu9528

金虫 (小有名气)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
实际上做多维信号的稀疏表示的时候,用到就是类似的思想吧。
4楼2013-06-03 18:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

三木_1987

金虫 (初入文坛)


★
小木虫: 金币+0.5, 给个红包,谢谢回帖
去年 我试过 仿真可以得到合理可解释的结果  但是没去证明

[ Last edited by 三木_1987 on 2013-5-31 at 21:44 ]
2楼2013-05-31 21:41:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 12/600 2026-10-13 20:40 by Merryya
[硕博家园] 极限拉扯,如何抉择 +9 Isak0320 2026-10-12 13/650 2026-10-13 11:57 by Isak0320
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见