1234567891011121314151617181920212223 |
- function modelmat=creatcontmodel(euclidean)
- th=linspace(-pi,pi,17);
- th=th(1:end-1);
- [x,y]=pol2cart(th',1); % draw unit circle in Cartesian coords
- xy=[x y]; % cat for easier indexing below
- for a=1:length(th)
- for b=1:length(th)
-
- if euclidean==0
- modelmat(a,b)=(abs(angdiff(th(a),th(b)))/pi); %angular
- elseif euclidean==1
- modelmat(a,b)=norm(xy(a,:) - xy(b,:)); %euclidean
- end
-
- end
- end
- % figure;imagesc(modelmat);
- end
|