MATLAB仿真程序代碼--模擬品牌傳播的過程

1、附錄MATLAB 仿真程序代碼function Et,Dt=BrandSpread(u0,params)% 這個程序用來模擬品牌傳播的過程% n 研究的區(qū)域的大小，為 n*n 的點(diǎn)陣% t 時刻% u(x,y,t) t時刻(x,y)點(diǎn)的品牌信心值% Q(x,y,t) t時刻(x,y)點(diǎn)的品牌信號值% c(x,y) (x,y)點(diǎn)的品牌傳播速率% k(x,y) (x,y)點(diǎn)對品牌信號的敏感程度% Initializationif nargin2 params=strcut;end;if isfield(params,con) & isempty(params.con)con = params.co

2、n; 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 = 1O;end;if isfield(params,maxt) & isempty(params.maxt) maxt = para

3、ms.maxt;elsemaxt = 2O;end;t=1; % 循環(huán)次數(shù)tflag=1; % tflag=1時說明未達(dá)到平衡，需要做交換；tflag=O時則達(dá)到平衡，區(qū)域內(nèi)不會出現(xiàn)信號交換，停止信號交換u=uO;% 隨機(jī)生成 1,10 上的均勻分布 k(:,:)=round(9*rand(n)+1;% 根據(jù) u 和 k 計(jì)算每一個點(diǎn)上的 QQ(:,:)=u(:,:).*k(:,:);% 順序式傳播，% 先由點(diǎn) (i,j) 與其右邊的點(diǎn) (i,j+1) 進(jìn)行交換% 再由點(diǎn) (i,j) 與其下方的點(diǎn) (i+1,j) 進(jìn)行交換E(t)=sum(sum(u)/nA2; %D(t)=sum(sum(u

4、-E(t)A2)/nA2;% 隨機(jī)獲得一個需要觀察的點(diǎn) (x_s,y_s) x_s=round(rand*n) y_s=round(rand*n) f(1)=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=1;endQ(i,j)=Q(i,j)-del

5、taQ-abs(deltaQ)*con;Q(i+1,j)=Q(i+1,j)+deltaQ-abs(deltaQ)*con;u(i,j)=Q(i,j)/k(i,j);u(i+1,j)=Q(i+1,j)/k(i+1,j);endend % Phase 2 for (i=1:n)for (j=1:n-1) deltaQ=(k(i,j+1)*Q(i,j)-k(i,j)*Q(i,j+1)/(k(i,j+1)+k(i,j); if deltaQ=vdeltaQ=v;endif (deltaQ=0) tflag=1;endQ(i,j)=Q(i,j)-deltaQ-abs(deltaQ)*con;Q(i,j+

6、1)=Q(i,j+1)+deltaQ-abs(deltaQ)*con;u(i,j)=Q(i,j)/k(i,j);u(i,j+1)=Q(i,j+1)/k(i,j+1);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=0 break;end;end;figure;plot(1:t,E,b-,linewidth,1); xlabel(time);ylabel(mean value); legend(change of mean value);figure; plot(1:t,D,r-,linewidth,1); xlabel(time); ylabel(variance); legend(