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

《MATLAB仿真程序代码--模拟品牌传播的过程.docx》由会员分享，可在线阅读，更多相关《MATLAB仿真程序代码--模拟品牌传播的过程.docx（5页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、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)点对品牌信号的敏感程度% Initializationif nargin2params=strcut;end;if isfield(params/con) & -isempty(params.

2、con) con = params.con;elsecon = 0;end;if isfield(params/maxf) & -isempty(params.maxf) maxf = params.maxf;elsemaxf = 10;end;if isfield(params/v) & -isempty(params.v)v = params.v;elsev = 1;end;if isfield(params/n) & -isempty(params.n)n = params.n;elsen = 10;end;if isfield(params/maxt) & -isempty(param

3、s.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计算每一个点上的Q Q(:,:)=u(:,:).*k(:,:);%顺序式传播,%先由点(ij)与其右边的点(ij+l)进行交换%再由点(ij)与其下方的点(i+Lj)进行交换 E(t)=sum(sum(u)/nA2; % D(t)=sum(sum(u-E

4、(t).A2)/nA2;%随机获得一个需要观察的点(x_s,y_s) x_s=round(rand*n) y_s=round(rand*n) f=u(x_s,y_s);fprintf(The coordinate of the observe point isn);fprintf( %d , %d )n,x_s,y_s);fprintf(itt mean valuet Variancetn);fprintf(%dt %,3et %,3et n,t,E,D);if t=vdeltaQ=v;endif (deltaQ=0)tflag = l;endQ(ij)=Q(ij)-deltaQ-abs(de

5、ltaQ)*con;Q(i+1 J)=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% Phase 2for (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);if deltaQ=vdeltaQ=v;endif (deltaQ-=0)tflag = l;endQ(ij); Q(i,j)-deltaQ-abs(deltaQ)*con;Q(ij + l)=Q(ij + l)+de

6、ltaQ-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);if t = maxt)break;end;if tflag =0break;end;end;figure;plot(l:t/E,b-,/linewidth,/l);xlabel(time);ylabel(mean value);legend(change of mean value);figure;plotaxD/r-/linewidthJ);xlabel(time);ylabel(variance);legend(change of variance);figure;plot(l:t,f,g-/linewidth,fl);xlabel(time);ylabel(U);legend(change of U);