MATLAB仿真程序代码-模拟品牌传播的过程

MATLAB仿真程序代码--模拟品牌传播的过程附录二MATLAB仿真程序代码function[Et/Dt]=BrandSpread(uO/params)%这个程序用来模拟品牌传播的过程%n研究的区域的大小，为n*n的点阵%t时刻%u(x,y,t)t时刻(x,y)点的品牌信心值%Q(x,y,t)t时刻(x,y)点的品牌信号值%c(x,y)(x,y)点的品牌传播速率%k(x,y)(x,y)点对品牌信号的敏感程度%Initializationifnargin<2params=strcut;end;ifisfield(params/con')&&-isempty(params.con)con=params.con;elsecon=0;end;ifisfield(params/maxf')&&-isempty(params.maxf)maxf=params.maxf;elsemaxf=10;end;ifisfield(params/v')&&-isempty(params.v)v=params.v;elsev=1;end;ifisfield(params/n')&&-isempty(params.n)n=params.n;elsen=10;end;ifisfield(params/maxt')&&-isempty(params.maxt)maxt=params.maxt;elsemaxt=20;end;t=l；%循环次数tflag=l;%tflag=l时说明未达到平衡，需要做交换;tflag=0时则达到平衡，区域内不会出现信号交换，停止信号交换u=u0;%随机生成［L10］上的均匀分布k(:z:)=round(9*rand(n))+l;%根据u和k计算每一个点上的QQ(:,:)=u(:,:).*k(:,:);%顺序式传播,%先由点(ij)与其右边的点(ij+l)进行交换%再由点(ij)与其下方的点(i+Lj)进行交换E(t)=sum(sum(u))/nA2;%D(t)=sum(sum((u-E(t)).A2))/nA2;%随机获得一个需要观察的点(x_s,y_s)x_s=round(rand*n)y_s=round(rand*n)f⑴=u(x_s,y_s);fprintf('Thecoordinateoftheobservepointis\n');fprintf('(%d,%d)\n',x_s,y_s);fprintf('it\tmeanvalue\tVariance\t\n');fprintf('%d\t%,3e\t%,3e\t\n',t,E,D);ift<=maxffigure;x=l:n;y=x;[X,Y]=meshgrid(x,y);mesh(X,Y/u(:/:));end;while(tflag==l)t=t+l;tflag=O;%Phase1for(i=l:n-l)for(j=l:n)deltaQ=(k(i+Lj)*Q(i,j)-k(i,j)*Q(i+Lj))/(k(i+Lj)+k(i,j));ifdeltaQ>=vdeltaQ=v;endif(deltaQ〜=0)tflag=l;endQ(ij)=Q(ij)-deltaQ-abs(deltaQ)*con;Q(i+1J)=Q(i+lj)+deltaQ-abs(deltaQ)*con;u(iJ)=Q(i,j)/k(i,j);u(i+lj)=Q(i+lj)/k(i+lj);endend%Phase2for(i=l:n)for(j=l:n-l)deltaQ=(k(ij+l)*Q(i,j)-k(iJ)*Q(iJ+l))/(k(iJ+l)+k(i,j));ifdeltaQ>=vdeltaQ=v;endif(deltaQ-=0)tflag=l;endQ(ij);Q(i,j)-deltaQ-abs(deltaQ)*con;Q(ij+l)=Q(ij+l)+deltaQ-abs(deltaQ)*con;u(ij)=Q(ij)/k(ij);u(ij+l)=Q(ij+l)/l<(ij+l);endendE(t)=sum(sum(u))/nA2;D(t)=sum(sum((u-E(t)).A2))/nA2;f(t)=u(x_s/y_s);ift<=maxffigure;x=l:n;y=x;[X,Y]=meshgrid(x/y);mesh(X,Y,u(:,:));end;fprintf('%d\t%.3e\t%.3e\t\n)t,E(t),D(t));if(t>=maxt)break;end;iftflag==0break;end;end;figure;plot(l:t/E,'b-,/'linewidth,/l);xlabel('time');ylabel('meanvalue');legend('changeofmeanvalue');figure;plotaxD/r-'/linewidth'J);xlabel('time');ylabel('variance');