24小时热门版块排行榜    

查看: 1299  |  回复: 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的回帖

三木_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的回帖

taoyu9528

金虫 (小有名气)



小木虫: 金币+0.5, 给个红包,谢谢回帖
实际上做多维信号的稀疏表示的时候,用到就是类似的思想吧。
4楼2013-06-03 18:49:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
2013-05-31 22:58   回复  
相关版块跳转 我要订阅楼主 sololee2010 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 这样的filecode谁见过 +3 布布和一二 2026-08-08 5/250 2026-08-08 12:22 by 仁砚薪传
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] fileCode有新解读? +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见