看这个构造有错误的数据对不对

prince00840

我打算构造一组数，再往里面增加错误数据。来检测程序的稳定性。

思路是利用两个矩阵相加的方式：
首先生成一个正态分布的矩阵，1000个数 : a =random('Normal',30,0.05,1,1000);
再生成一个随机的100个数 : b=random(1,100)？？？？？？？
两者加起来c=a+b

问题出在b的处理上面。a是一个1000个数，b就100个数
是不是两者不能相加啊。

[ 本帖最后由 ChaChing 于 2009-2-13 21:18 编辑 ]

ChaChing

可能所学背景不同, 不太清楚楼主要的东东!
c(=a+b)是1*1000 or 1*1100的矩阵? 请LZ说明确定一下!

prince00840

我感觉第二个矩阵b 应该是1000个数，其中随机100个数非零，另外的都是零，才能做到和a相加。
目的打算构造一组数，一共1000个数，其中包涵着100个错误



:@( 表达老不明确呜呜

ChaChing

依个人理解的, 猜测一下LZ要的!
原始数据a(1*1000), 欲加上一错误数据b(1*1000), 来检测程序的稳定性!
其中b含100个随机数, 其它为零(即维持a原本之值)
是这样吗?

prince00840

是啊是啊:loveliness:
终于能让别人理解出来了，感动ing:'(
现在不明白这个矩阵b怎么写？？

[ 本帖最后由 prince00840 于 2009-2-15 13:00 编辑 ]

ChaChing

具体100个随机数大小及摆放位置的准则, 个人并不清楚应该为何? LZ有规定吗?
个人先随意假设, 并利用另一随机数列, 摆放位置为这数列前100最小数! 具体应如何, LZ根据你的需求适当调整!
clc; clear;
a =random('Normal',30,0.05,1,1000);
cc=rand(1,1000); [cc,ind]=sort(cc);
b=zeros(1,1000); b(ind(1:100))=random('Normal',0,0.5,1,100);
c=a+b;

prince00840

没有要求随机的。
就是想把
矩阵b原来是全为零的1000个数，现在使其中的100个不为零:@)
那么还用sort函数吗？

ChaChing

看LZ的需求吧!
只是"其中的100个不为零", 如何决定那100个(前100个,後100个,中间100个...)? 所以我才想到使用随机给!
刚得知被送个礼物(鸡蛋), 也不晓得原因, 心情有些不好!

prince00840

我给你朵鲜花嘿嘿:lol
下面是您的程序画的图，看着是随机产生的:@) 十分感谢，我就是矩阵b那里不会来

[ 本帖最后由 prince00840 于 2009-2-16 08:49 编辑 ]