24小时热门版块排行榜    

查看: 1955  |  回复: 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 lyz123lyz 2026-08-27 4/200 2026-08-27 13:08 by kulium
[基金申请] 基金不中,共勉 +8 eulota 2026-08-26 8/400 2026-08-27 12:34 by lice_1987
[基金申请] 基金未中,这种答复是模板吗? +4 zhaosm1982 2026-08-27 5/250 2026-08-27 12:17 by zhaosm1982
[基金申请] 投票:  有多少人是今天查系统知道结果的? +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 11:29 by xxniao123
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 让我中一个面上吧! +13 大萍1987 2026-08-20 17/850 2026-08-27 08:20 by 大萍1987
[基金申请] 国自然面上复盘~欢迎讨论 +7 晴天加油 2026-08-26 8/400 2026-08-27 07:53 by longfie172
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 明天应该可查了!? +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
[基金申请] 我面上完蛋了 +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
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 应该是下周三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
信息提示
请填处理意见