[2dPIC调试笔记]初始化变量1014(2)

%**********初始化变量*********%
xlen = nx;         %XÍø¸ñÊý
ylen = ny;         %YÍø¸ñÊý
nxny = nx*ny;      %×ÜÍø¸ñÊý
g = 1.0;
npt =sum(np(1:ns));%×Ü
b0 = wc/qm(1);costh = cos(theta*pi/180);
sinth = sin(theta*pi/180);
cosph = cos(phi  *pi/180);
sinph = sin(phi  *pi/180);bx0 = b0*sinth*cosph;
by0 = b0*sinth*sinph;
bz0 = b0*costh;
csq = cv*cv;
tcs = 2.0*csq;
it = 0;nxp1 = nx+1;nyp1 = ny+1;nxp2 = nx+2;nyp2 = ny+2;nxp3 = nx+3;nyp3 = ny+3;nxc = nx/2+1;nyc = ny/2+1;nxny = nx*ny;bx = ones(nxp2, nyp2)*bx0;by = ones(nxp2, nyp2)*by0;bz = ones(nxp2, nyp2)*bz0;ex = zeros(nxp2, nyp2);ey = zeros(nxp2, nyp2);ez = zeros(nxp2, nyp2);rho = zeros(nxp2, nyp2);psi = zeros(nxp2, nyp2);%  definition of inductive electric field
%  exind = zeros(nxp2, nyp2);
%  exind(nx/4+1:nx/4*3, ny/8*3+1:ny/8*5) = ones(nx/4*2, ny/8*2)*exinductive;% Indiex Arry for Matrix OperationX1 = 1:nx;Y1 = 1:ny;X2 = 2:nxp1;Y2 = 2:nyp1;X3 = 3:nxp2;Y3 = 3:nyp2;V1 = 1:nxp1;V2 = 2:nxp2;VA = 1:nxp3;W1 = 1:nyp1;W2 = 2:nyp2;WA = 1:nyp3;A2 = 2:nx;B2 = 2:ny;A3 = 3:nxp1; B3 = 3:nyp1;
% Paritlce Initializationx = zeros(1,npt);y = zeros(1,npt);vx = zeros(1,npt);vy = zeros(1,npt);vz = zeros(1,npt);n2 = 0;
for  k=1:nsn1 = n2 + 1;n2 = n2 + np(k);pchr = pi/180.0*pch(k);vdpa = vd(k)*cos(pchr);vdpe = vd(k)*sin(pchr);x(n1:n2) = xlen*rand(1,np(k));% y(n1:n2) = ylen*rand(1,np(k));y(n1:n2) = ylen/np(k)*(1:np(k));%x(n1:n2) = xlen*0.5+0.1*xlen*randn(1,np(k));%y(n1:n2) = ylen*0.5+0.1*ylen*randn(1,np(k));for i = n1:n2phase = 2*pi*rand;uxi = vpe(k)*randn + vdpe*cos(phase);uyi = vpe(k)*randn + vdpe*sin(phase);uzi = vpa(k)*randn + vdpa; % rotation to the direction of the magnetic fieldux = uxi*costh*cosph - uyi*sinph + uzi*sinth*cosph;uy = uxi*costh*sinph + uyi*cosph + uzi*sinth*sinph;uz = -uxi*sinth                  + uzi*costh;%g = cv /sqrt(csq +ux*ux +uy*uy +uz*uz);vx(i) = ux*g;vy(i) = uy*g;vz(i) = uz*g;end
end
xcen= xlen*0.5;
ycen= ylen*0.5;
%for m = 1:np(1)
%  r = sqrt((x(m)-xcen)^2+(y(m)-ycen)^2);
%   if (r < 4.0) & (r >0.1)
%      vx(m) = (x(m) - xcen)*vd(1)/r;
%      vy(m) = (y(m) - ycen)*vd(1)/r;
%   end
%end
%
% Coefficiendts for Poisson Solver
rkfact = zeros(nx,ny);
rkxmin = pi/nx;
rkymin = pi/ny;
nxh = nx/2;
nyh = ny/2;
fftx = 1.0/nxh;
ffty = 1.0/nyh;
for j=3:nyrky  = sin(rkymin*floor((j-1)/2))*2.0;for i=1:nx rkx  = sin(rkxmin*floor((i-1)/2))*2.0;rkfact(i,j) = 1.0/(rkx^2+rky^2);endrkx = 2.0;rkfact(2,j) = 1.0/(rkx^2+rky^2);
end
for i=3:nx rkx  = sin(rkxmin*floor((i-1)/2))*2.0;rkfact(i,1) = 1.0/(rkx^2);
end
rky = 2.0;
for i=3:nx rkx  = sin(rkxmin*floor((i-1)/2))*2.0;rkfact(i,2) = 1.0/(rkx^2+rky^2);
end
rkfact(1,2) = 1.0/(rky^2);
rkx = 2.0;
rkfact(2,2) = 1.0/(rkx^2+rky^2);
rkfact(2,1) = 1.0/(rkx^2);
rkfact(1,1) = 0.0;
%
q = (xlen*ylen)./np(1:ns).*(wp(1:ns).^2)./qm(1:ns);
mass = q./qm(1:ns);
rho0 = - sum(q(1:ns).*np(1:ns))/nxny*ones(nxp2,nyp2);color = get(gca,'ColorOrder');
fign = 1;
itime = 0;
jtime = 0;
rmax = dx*nx;
rmin = 0.0;
%ke=zeros(ntime,ns);

[2dPIC调试笔记]初始化变量1014(2)相关推荐

调试笔记--jlink 变量转实时波形小技巧
调试笔记–jlink 变量转实时波形小技巧上篇讲了jlink RTT组件在RAM中开辟一段内存,将printf字符串存入然后通过jlink调试接口读取RAM,将printf字符串在jlink-RTT ...
[2dPIC调试笔记]输入参数归一化1014(1)
%**********输入参数归一化*********% rnx = dx; % distance rnt = dt/2; % time rnv = rnx/rnt; % velocity rne = ...
MCP2515调试笔记----初始化注意事项
在调试基于F103单片机的MCP2515时,怎么调试都不正常,用的还是以前调试过的代码,隔了一段时间后再重新运行就不行了: 于是单步运行调试,没想到,单步运行调试竟然成功了:于是想到是延时的问题,在初 ...
[C++调试笔记]初始化3种粒子数据initmaxw
#include <stdio.h> #include "define.h" /* 坐标变换 Michael Komm,Interakce plazmatu se st ...
[C++调试笔记]初始化
/* 初始化 *///static const int numb_e = 1.0e6; //模拟电子的超粒子数//static const int numb_D = 1.0e6; //模拟D的超粒子数 ...
[2DPIC调试笔记]parameter_antenna_radiation1013(3)
% Parameters for KEMPO2 % For a continuous job, jobnumber > 1 jobnumber = 1; if jobnumber == 1cle ...
OMAPL138调试笔记
title: Linux 嵌入式使用 tags: 新建,模板,小书匠 grammar_cjkRuby: true Copyright (C) @2018 Widic Filename:File-fun ...
android 调试笔记
android 调试笔记 adb使用: 连接不上手机 C:\Users\Administrator\.android\adb_usb.ini 增加一行内容: 0x1782 安装第三方应用 $ a ...
调试笔记--keil 测量周期小技巧
调试笔记–keil 测量周期小技巧本文参考安富莱专题教程第7期 http://www.armbbs.cn/forum.php?mod=viewthread&tid=87176&ext ...

[2dPIC调试笔记]初始化变量1014(2)

[2dPIC调试笔记]初始化变量1014(2)相关推荐

最新文章

热门文章