24小时热门版块排行榜    

查看: 1918  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 娱乐 +3 Tide man 2026-08-03 3/150 2026-08-04 09:58 by gltch
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 5/250 2026-08-04 09:27 by v1TJC4qJauAB
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 cu9Nq1xK233Z 2026-08-03 10/500 2026-08-04 09:15 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 10/500 2026-08-04 08:05 by 88817753
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +9 锦衣卫寒战 2026-07-28 13/650 2026-08-04 07:21 by Ermito
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 cu9Nq1xK233Z 2026-08-03 6/300 2026-08-04 06:10 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 8/400 2026-08-04 03:21 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[考博] 2027年申博 50+3 射雕英雄胜 2026-07-30 3/150 2026-08-02 09:36 by lfy8008
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
信息提示
请填处理意见