24小时热门版块排行榜    

查看: 1961  |  回复: 14

abcxubu

新虫 (正式写手)

[求助] 线性方程组求解 已有2人参与

已知AX=b,A的行数小于列数,且X是一个稀疏的,能不能只求出X某几个位置的值?可以用广义逆么?十分感谢您的答复
回复此楼

» 猜你喜欢

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

好男儿,胸怀天下
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

askuyue

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
这是一个优化问题:
     min || x ||_{0}
         s.t. Ax = b
See:
       1) Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency information
       2) Stable signal recovery from incomplete and inaccurate measurements
        3) Near-optimal signal recovery from random projections: Universal encoding strategies?
HEHE
11楼2016-04-08 13:28:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

askuyue

铁杆木虫 (著名写手)

【答案】应助回帖

还有一个点就是:
     广义逆
的意义,要弄明白是几个意思,很自然地回牵涉到一正则化的问题。
HEHE
12楼2016-04-08 13:30:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

abcxubu

新虫 (正式写手)

或者有什么相关方面的书籍和资料能推荐一下么?谢谢~
好男儿,胸怀天下
2楼2016-04-06 21:02:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
Not possible. More constraints are needed

[ 发自手机版 http://muchong.com/3g ]
小木虫给我温暖,给我希望,爱就要爱小木虫。
3楼2016-04-07 00:46:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abcxubu

新虫 (正式写手)

引用回帖:
3楼: Originally posted by zaq123321 at 2016-04-07 00:46:42
Not possible. More constraints are needed

麻烦问一下,先对矩阵A求伪逆,然后选出所要求解位置的行,与b相乘不可以么?(x除某几个位置外都是零),或者先把要求位置的A的列组成一个矩阵,求伪逆也不行么?谢谢啦
好男儿,胸怀天下
4楼2016-04-07 09:17:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zaq123321

专家顾问 (著名写手)

【答案】应助回帖

You need to put further constraint on the solution since the solution is not unique. Usually put like 2-norm minimal etc.
小木虫给我温暖,给我希望,爱就要爱小木虫。
5楼2016-04-07 10:21:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abcxubu

新虫 (正式写手)

引用回帖:
5楼: Originally posted by zaq123321 at 2016-04-07 10:21:39
You need to put further constraint on the solution since the solution is not unique. Usually put like 2-norm minimal etc.

您的意思是,如果将b-Ax'的二范数作为目标函数,并求min,就可以用伪逆形式?在实验中我用了选取A中要求的列,作为矩阵,求其伪逆的方法,求x某些位置的值。当选取的列数小于A行数的时候,结果比较好;选的列数大于A的行数的时候,结果较差,不知为什么(自己已将min||b-Ax'||作为目标函数,x'为近似解),百思不得其解,希望大神帮忙解释一下。(列<行,列满秩,用的(T'T)^(-1)T';反之,行满秩,用的T'(TT')^(-1))
好男儿,胸怀天下
6楼2016-04-07 14:25:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abcxubu

新虫 (正式写手)

引用回帖:
5楼: Originally posted by zaq123321 at 2016-04-07 10:21:39
You need to put further constraint on the solution since the solution is not unique. Usually put like 2-norm minimal etc.

刚刚那个A是用的高斯随机矩阵,当A是哈达玛矩阵的时候,便可以先求A整体的伪逆,然后用A的某些行与b相乘,求出x对应位置的值。这种情况下,无论是行大于列还是列大于行,都能得到很好的结果,是因为哈达玛矩阵很好的正交性么?(第二种情况只有A为哈达玛的时候才成立),谢谢您的回复!!!
好男儿,胸怀天下
7楼2016-04-07 14:40:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abcxubu

新虫 (正式写手)

8楼2016-04-07 21:01:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

200971014

禁虫 (小有名气)

本帖内容被屏蔽

9楼2016-04-08 00:22:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

abcxubu

新虫 (正式写手)

引用回帖:
9楼: Originally posted by 200971014 at 2016-04-08 00:22:26
方程数目必须等于未知数数目,也就是说A必须是方阵

十分感谢您的解答
好男儿,胸怀天下
10楼2016-04-08 08:03:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 abcxubu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[基金申请] 基金不中,共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[基金申请] 国自然面上复盘~欢迎讨论 (金币+5) +10 晴天加油 2026-08-26 11/550 2026-08-27 16:35 by msuhlx
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧? +4 哈哈蛤? 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见