Matlab 求解电场分析和总结.docxVIP

function Vaule_Arry=Get_Value(Maxvaule,Divide,Num) % function: Vaule_Arry=Get_Value(Maxvaule,Divide,Num) % Maxvaule为定义的最大电势值，Divide为Maxvaule的等分数，Num为递归的次数 % Vaule_Arry返回的是递归后的矩阵 % 此函数也可以把电势分布图画出 clc n=Maxvaule:(-(Maxvaule/(Divide-1))):0; %产生Divide等分temparry=ones(1,size(n,2)); %得到n的列数 Temp_Memory =n*temparry; %构建数据存储矩阵 for loopnum=0:Num-1 %循化递归Temp_Memory (2:size(n,2),[1 size(n,2)])=0; %构建边界Temp_Memory (size(n,2),:)=0; %构建边界 %********************************************************************* %** circshift(Temp_Memory ,1)把Temp_Memory 矩阵循环往下移一位， ** %** circshift(Temp_Memory ,-1)把Temp_Memory 矩阵循环往上移一位。 ** %** circshift(Temp_Memory ,[0 1])把Temp_Memory 矩阵循环往右移一位， ** %** circshift(Temp_Memory ,[0 -1])把Temp_Memory 矩阵循环往左移一位。 ** %********************************************************************* Temp_Memory =circshift(Temp_Memory ,1)+circshift(Temp_Memory ,-1)+..., circshift(Temp_Memory ,[0 -1])+circshift(Temp_Memory ,[0 1]); %Temp_Memory得到每一点的上下左右矩阵元素相加的和。 Temp_Memory =Temp_Memory /4; %把Temp_Memory除4得平均 Temp_Memory (2:size(n,2),[1 size(n,2)])=0; %构建边界 Temp_Memory (size(n,2),:)=0; %构建边界 Temp_Memory end (1,:)=Maxvaule; %构建边界 disp([定义的最大电势值 Maxvaule=, num2str(Maxvaule)]); disp([Maxvaule的等分数 Divide=, num2str(Divide)]); disp([递归的次数 Num,=num2str(Num)]); xaxi=size(n,2)-1:-1:0; xtemp=ones(1,size(xaxi,2)); xaxi=xtemp*xaxi; %得到画图的X坐标 hold on for plotcount=1:size(n,2) plot(xaxi(plotcount,:),Temp_Memory (plotcount,:)); %画得到的最终数据的图end xlabel(X) ylabel(电势) title([电势分布图 最大电势: ,num2str(Maxvaule), 等分数: ,..., num2str(Divide), hold off 递归次数: ,num2str(Num)]); Vaule_Arry=Temp_Memory ; %返回最终结果Temp_Memory grid on 因为定义的函数有Maxvaule；Divide；Num三个变量（见图 1） 所以在仿真的时候要定义好变量值。 图 1 以书上164页的例 3.7.1为例子： 例 3.7.1 定义的电势是100V，共有25个点也就是分成了5等份我们假设只递归1次所以就有了如下的函数初始化（见图 2）。 T=Get_Value(100,5,1) 然后按下回车键就可以执行代码。 图 2 运行后的到的只递归1次后的结果（图 4,表格 1）和书上所给的结果一致说明是正确的。 图 3 表格 1 图 4 图 5是递归2次后的结果和书上的结果一致（见书164页数据），（详见 表格2） 图 5 表格 2 图 6 是电势为100V等分5份递归28次，函数的运行格式设置，按回车执行程序得到 图 7 的运行结果和书上的结果一致（见书165页数据）（详见 表格 3）。 图 6 图 7