matlab实验报告

来源:百度知道 编辑:UC知道 时间:2024/05/07 03:49:06
空心导体内,外半径分别是a和b,若在球心放置一点电荷Q,求各处的点位和电势。。。。求MTALAB仿真程序。

clc;
  clear;
  %对空心球体的空间电位的仿真
  %首先绘制两个球面,作为空心球体的模型
  %假设内半径a=5m,外半径b=10m
  %绘制内球面
  a=5;b=10;
  x=-5:0.5:5;
  y=-5:0.5:5;
  [X1,Y1]=meshgrid(x,y);
  Z11=sqrt(a^2-X1.^2-Y1.^2);
  index=find(X1.^2+Y1.^2>a^2);
  Z11(index)=NaN;
  plot3(X1,Y1,Z11,'LineStyle','.','Color','g');
  hold on;
  Z12=-sqrt(a^2-X1.^2-Y1.^2);
  index=find(X1.^2+Y1.^2>a^2);
  Z12(index)=NaN;
  plot3(X1,Y1,Z12,'LineStyle','.','Color','g');
  hold on;
  %绘制外球面
  x=-10:1:10;
  y=-10:1:10;
  [X2,Y2]=meshgrid(x,y);
  Z21=sqrt(b^2-X2.^2-Y2.^2);
  index=find(X2.^2+Y2.^2>b^2);
  Z21(index)=NaN;
  plot3(X2,Y2,Z21,'LineStyle','.','Color','r');
  hold on;
  index=find(X2.^2+Y2.^2>b^2);
  Z22=-sqrt(b^2-X2.^2-Y2.^2);
  Z22(index)=NaN;
  plot3(X2,Y2,Z22,'Line