西北大学《MATLAB基础与应用》课件-第8章matlab在理工课程中的应用.pptVIP

结果2：输入电荷数目N=:3 ----- 对电荷 # 1/n输入电荷位置[x y](米):[1 1] 输入电荷量(库仑):0.5 ----- 对电荷 # 2/n输入电荷位置[x y](米):[2 3] 输入电荷量(库仑):0.7 ----- 对电荷 # 3/n输入电荷位置[x y](米):[3 4] 输入电荷量(库仑):0.9 C0 = 8.9918e+009 其它电荷作用在电荷 # 1上的合力为: x-分量:-4.54141e+008 N y-分量:-8.21956e+008 N 其它电荷作用在电荷 # 2上的合力为: x-分量:-1.72133e+009 N y-分量:-1.43985e+009 N 其它电荷作用在电荷 # 3上的合力为: x-分量:2.17548e+009 N y-分量:2.2618e+009 N [例] 由电位的表示式计算电场并画出等电位线和电场方向 建模 如果已知空间的电位分布 则空间的电场等于电位场的负梯度 其中 分别为x,y,z三个方向的单位向量 MATLAB设有gradient函数，它是靠数值微分，因此空间观测点应该取得密一些，已获得较高的精度 fprintf(输入电位分布方程 V(x,y) \n); fprintf(例如: log(x.^2 + y.^2) \n); V = input(: ,s); % 读入字符串 V(x,y) NGrid = 20; % 绘图的网格线数 xMax = 5; % 绘图区从 x= -xMax 到 x= xMax yMax=5; xPlot = linspace( -xMax, xMax, NGrid); % 绘图取的x值 [x,y]=meshgrid(xPlot); VPlot=eval(V); [ExPlot, EyPlot] = gradient(-VPlot); % 电场是电位的负梯度 clf; subplot(1,2,1),meshc(VPlot); % 画含等高线的三维曲面 set(gcf,color,‘w) % 置图形背景色为白色 xlabel(x); ylabel(y); zlabel(电位); % 规定等高线图的范围及比例 subplot(1,2,2), axis([-xMax xMax -yMax yMax]); % 建立第二子图 cs = contour(x,y,VPlot); % 画等高线 clabel(cs); hold on; % 在等高线图上加上编号等高线值 %@ 在等高线图上加上电场方向 quiver(x,y,ExPlot,EyPlot); % 画电场 E 的箭头图 xlabel(x); ylabel(y);hold off; 练习目的： 输入，输出 产生等差向量 由向量生成网格点坐标 表达式取值 清理图形窗口 图形窗口分割 三维曲面绘图 设置句柄属性 坐标轴的添加 坐标轴的区间设置 箭头图形的绘制 电位三维立体图 等位线及电场分布图 2.5 稳恒磁场 [例]用毕奥－－萨伐定律计算电流环产生的磁场 z ? Idl dB dB P ? ? r I ?0真空的磁导率?0＝4??10-7H/m。 任一电流元Idl在空间任一点P处所产生的磁感应强度为 clear all; % 清工作空间及变量初始化 mu0 = 4*pi*1e-7; % 真空导磁率 (T*m/A) I0 = 5.0; % 环中电流(A) Rh = input(输入环半径Rh(m):(取2) ); C0 = mu0/(4*pi) * I0; % 归并常数 xMax = 3; yMax = 3; % 规定图的范围 NGx = 21; NGy = 21; % 规定观测点网格线数 x=linspace(-xMax, xMax, NG