Large-scale constrained linear least-squares optimization


Slide 3







We can simulate the effect of vertical motion blurring by averaging each pixel with the 5 above and below. We construct a sparse matrix D, shown in part above, that will do this with a single matrix multiply.
>> blur=5; mindex=1:mn; nindex=1:mn;
>> for i=1:blur,
>> mindex=[mindex i+1:mn 1:mn-i]; nindex=[nindex 1:mn-i i+1:mn];
>> end
>> D=sparse(mindex,nindex,1/(2*blur+1));