求助：关于图像显示


s1=zeros(511*10,1);        　%每组511个数据，一共十组
s1=textread('data1.txt')   　%  数据范围：0=<s1()<=256
s2=zeros(800,800);
s2=textread('data2.txt')      %  数据范围：0=<s2()<=900
要求:
     在800*800像素阵中填色，
当0=<s2()<510时，颜色大小由s1的大小控制,比如：s2(100,150)=211,则在(100,150)这个位置填上大小为s1(211)的颜色;
当s2()>510的区域，用白色填充
　　我用几次循环，plot 画，发现速度太慢，慢的结果出不来，这样做行不通；
　　现在想用像素直接来显示图像，不知怎么做？
　　请大家教教我应该怎么做。
　　谢谢！