200论坛币求助一个IML程序如何写 [推广有奖]

匿名网友

楼主

匿名网友 发表于 2012-11-30 11:15:39 |AI写论文

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;

分享0 收藏0 回帖

关键词：200论坛币 0论坛币如何写论坛币 Orthogonal 程序论坛如何 orthogonal original

本帖被以下文库推荐

· SAS精彩问答|主题: 2530, 订阅: 30

加关注串个门加好友发消息 0关注 52粉丝院士 ziyenano 当前离线阅读权限 255 威望 2 级论坛币 6137 个通用积分 27.3565 学术水平 338 点热心指数 326 点信用等级 298 点经验 57352 点帖子 879 精华 1 在线时间 4132 小时注册时间 2012-1-4 最后登录 2025-5-7 雷达卡	沙发 ziyenano 发表于 2012-11-30 11:15:40 此帖仅作者可见

	举报

加关注串个门加好友发消息 0关注 52粉丝院士 ziyenano 当前离线阅读权限 255 威望 2 级论坛币 6137 个通用积分 27.3565 学术水平 338 点热心指数 326 点信用等级 298 点经验 57352 点帖子 879 精华 1 在线时间 4132 小时注册时间 2012-1-4 最后登录 2025-5-7 雷达卡	藤椅 ziyenano 发表于 2012-11-30 14:15:58 此帖仅作者可见

	举报

加关注串个门加好友发消息 4关注 16粉丝教授 wqf_cufe 当前离线阅读权限 255 威望 0 级论坛币 4283 个通用积分 6.5000 学术水平 10 点热心指数 17 点信用等级 13 点经验 565 点帖子 730 精华 0 在线时间 2236 小时注册时间 2009-5-26 最后登录 2021-12-29 雷达卡	板凳 wqf_cufe 发表于 2012-11-30 14:35:40 此帖仅作者可见

	举报

加关注串个门加好友发消息 4关注 16粉丝教授 wqf_cufe 当前离线阅读权限 255 威望 0 级论坛币 4283 个通用积分 6.5000 学术水平 10 点热心指数 17 点信用等级 13 点经验 565 点帖子 730 精华 0 在线时间 2236 小时注册时间 2009-5-26 最后登录 2021-12-29 雷达卡	报纸 wqf_cufe 发表于 2012-11-30 14:45:17 此帖仅作者可见

	举报

加关注串个门加好友发消息 0关注 52粉丝院士 ziyenano 当前离线阅读权限 255 威望 2 级论坛币 6137 个通用积分 27.3565 学术水平 338 点热心指数 326 点信用等级 298 点经验 57352 点帖子 879 精华 1 在线时间 4132 小时注册时间 2012-1-4 最后登录 2025-5-7 雷达卡	地板 ziyenano 发表于 2012-11-30 14:47:37 此帖仅作者可见

	举报

加关注串个门加好友发消息 4关注 16粉丝教授 wqf_cufe 当前离线阅读权限 255 威望 0 级论坛币 4283 个通用积分 6.5000 学术水平 10 点热心指数 17 点信用等级 13 点经验 565 点帖子 730 精华 0 在线时间 2236 小时注册时间 2009-5-26 最后登录 2021-12-29 雷达卡	7楼 wqf_cufe 发表于 2012-11-30 14:49:00 此帖仅作者可见

	举报

加关注串个门加好友发消息 0关注 52粉丝院士 ziyenano 当前离线阅读权限 255 威望 2 级论坛币 6137 个通用积分 27.3565 学术水平 338 点热心指数 326 点信用等级 298 点经验 57352 点帖子 879 精华 1 在线时间 4132 小时注册时间 2012-1-4 最后登录 2025-5-7 雷达卡	8楼 ziyenano 发表于 2012-11-30 14:58:39 此帖仅作者可见

	举报

加关注串个门加好友发消息 4关注 16粉丝教授 wqf_cufe 当前离线阅读权限 255 威望 0 级论坛币 4283 个通用积分 6.5000 学术水平 10 点热心指数 17 点信用等级 13 点经验 565 点帖子 730 精华 0 在线时间 2236 小时注册时间 2009-5-26 最后登录 2021-12-29 雷达卡	9楼 wqf_cufe 发表于 2012-11-30 15:02:51 此帖仅作者可见

	举报

加关注串个门加好友发消息 0关注 52粉丝院士 ziyenano 当前离线阅读权限 255 威望 2 级论坛币 6137 个通用积分 27.3565 学术水平 338 点热心指数 326 点信用等级 298 点经验 57352 点帖子 879 精华 1 在线时间 4132 小时注册时间 2012-1-4 最后登录 2025-5-7 雷达卡	10楼 ziyenano 发表于 2012-11-30 15:09:42 此帖仅作者可见

	举报