全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 SAS专版
匿名
843 28
2012-11-30
悬赏 200 个论坛币 已解决
请大家帮忙解决一下这个IML程序如何写?谢谢!!!

1.  TheGram-Schmidt process is used to find an orthogonal basis that spans a givensubset of Rn from an original spanning set:


Given a basis {x1, …, xp}for a subspace W of  Rn,define


v1 = x1

v2 = x2– projv2(x2)

v3 = x3– projv1(x3)  – projv2(x3)

vp = xp– projv1(xp) - projv2(xp)- … - projvp-1(xp)


where projv(x) is the projection operator (inprod(u,v)/inprod(u,u))*u and inprodis the inner product.


Then {v1, …, vp} is anorthogonal basis for W and span{x1, …, xp}= span{v1, …, vp}.


Write a SAS prgoram to find an orthogonal basis for Wusing IML.  Transform it into anorthonormal basis.


            1 2 6

          -1 0 2

W =     1 4 3

            0 1 0

            1 2 0

sec6_4.pdf
大小:(758.99 KB)

 马上下载


最佳答案

ziyenano 查看完整内容

这是按上面的要求编写的 proc iml; w={1 2 6, -1 0 2, 1 4 3, 0 1 0, 1 2 0}; v=j(nrow(w),ncol(w),0); do i=1 to ncol(w); if i=1 then do; v[,1]=w[,1]; end; else do; u=w[,i]; vi=v[,1:i-1]; v[,i]=u-(t(u)*vi/((vi#vi)[+,])*t(vi))`; /*v[,i]=u-vi*(t(vi)*u/((vi#vi)[+,])`);*/ end; end; print v; quit;
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2012-11-30 11:15:40
此帖仅作者可见
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-30 14:15:58
此帖仅作者可见
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-30 14:35:40
此帖仅作者可见
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-30 14:45:17
此帖仅作者可见
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-11-30 14:47:37
此帖仅作者可见
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群