|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
因为处理的数组比较大,小弟只会用循环来做:
clear
global Z % Z 是NxNx3-dimension matrix
global dataY % Y 是NxN-dimension matrix
N=308;
h=input('Input :'); % h 最佳值为4.8
RGBZ=Z;
Y=dataY;
ziel=zeros(N);
for i=1:N % for循环运行时间太慢
for j=1:N
ziel(i,j)=myBild([Z(i,j,1) Z(i,j,2) Z(i,j,3)]',h,N,RGBZ,Y); %
end
end
myBild()函数,文件名myBild.m:
function k=myBild(u,h,N,RGBZ,Y)
mymulti=zeros(N);
for i=1:N
for j=1:N
mymulti(i,j)=myK((u-RGBZ(i,j))/h);
end
end
Y3=mymulti.*Y;
k=(1/(N*h^3))*sum(sum(Y3));
因为循环太多,所以运行太慢(当N=25时,运行时间为十分钟),求各位大虾帮帮忙,实在不知道怎么改:@L :@L :@L |
|