Matlab跨类修改属性：磁粒子交互模拟技术问询

我正在开展具有磁特性的粒子间交互作用的模拟工作，需要实现从一个类中更新另一个类的属性——因为目标属性仅在存在至少一个其他粒子时才会发生变化。

我设计了Particle基类，以及Proton、Electron等子类（当前聚焦Proton类），同时定义了GroupParticles类，其初始化输入为Particle对象数组。现有代码如下：

classdef Particle < handle
    properties
        location
        magneticmoment
    end

    methods
        % Get-Methoden definieren
        function res = GetLocation(obj)
            if nargin ~= 0
                res = obj.location;
            else
                Error('location is not set yet')
            end
        end

        function res = GetMagneticMoment(obj)
            if nargin ~= 0
                res = obj.magneticmoment;
            else
                Error('Magnetic Moment is not set yet')
            end
        end
    end
end

classdef Proton < Particle
    properties (Access = private)
        gamma = 2.6752218744 * 10^8             % 1/ (s*T)
        length_mu = 1.41060679736 * 10^(-26)    % J/T
    end
    properties
        ActingBField
    end

    methods
        % Constructor 
        function obj = Proton(location, direction)
            obj.location = location;
            obj.magneticmoment = direction / abs(direction) * obj.length_mu;
            obj.ActingBFeld = zeros(1,3);
        end
    end
end

function res = abs(a)
    res = sqrt(a(1)^2 + a(2)^2 + a(3)^2);
end

classdef GroupParticles
    properties
        ParticleArray
        DistanceMatrix
    end

    methods
        function obj = GroupParticles(Particles)
            obj.ParticleArray = Particles;
        end
    end
end

我的需求是在GroupParticles类中定义方法，修改Proton对象的magneticmoment和ActingBField属性。请问Matlab中是否可行？若不可行，还有哪些对象关联的实现方式？

一、完全可行，核心依托Matlab的handle类特性 #

你的Particle基类继承自handle，这意味着所有子类（比如Proton）都是引用类型对象。当GroupParticles持有ParticleArray时，操作数组里的对象就是直接操作原对象实例，修改它们的属性会直接生效。

二、具体实现步骤 #

1. 修正代码小问题
   首先注意Proton构造函数里的拼写错误：obj.ActingBFeld应该改为obj.ActingBField，和属性名保持一致。

2. 在GroupParticles中添加更新属性的方法
   示例实现一个updateParticleProperties方法，遍历粒子数组，计算并更新每个质子的磁场和磁矩：

   classdef GroupParticles
       properties
           ParticleArray
           DistanceMatrix
       end
   
       methods
           function obj = GroupParticles(Particles)
               obj.ParticleArray = Particles;
               % 可选：初始化距离矩阵
               obj.DistanceMatrix = obj.calculateDistanceMatrix();
           end
   
           % 计算粒子间距离矩阵
           function mat = calculateDistanceMatrix(obj)
               n = length(obj.ParticleArray);
               mat = zeros(n);
               for i = 1:n
                   loc_i = obj.ParticleArray(i).location;
                   for j = i+1:n
                       loc_j = obj.ParticleArray(j).location;
                       dist = sqrt(sum((loc_i - loc_j).^2));
                       mat(i,j) = dist;
                       mat(j,i) = dist;
                   end
               end
           end
   
           % 更新所有质子的磁场和磁矩
           function updateParticleProperties(obj)
               n = length(obj.ParticleArray);
               if n < 2
                   % 少于2个粒子时无需更新
                   return;
               end
   
               for i = 1:n
                   current_particle = obj.ParticleArray(i);
                   if ~isa(current_particle, 'Proton')
                       continue; % 仅处理质子
                   end
   
                   % 计算当前粒子受到的总磁场ActingBField
                   total_B = zeros(1,3);
                   for j = 1:n
                       if i == j
                           continue;
                       end
                       other_particle = obj.ParticleArray(j);
                       % 替换为你实际的磁场计算公式（示例为磁偶极子磁场简化版）
                       r = other_particle.location - current_particle.location;
                       r_mag = sqrt(sum(r.^2));
                       B = (3*(dot(other_particle.magneticmoment, r)*r)/r_mag^5) - (other_particle.magneticmoment)/r_mag^3;
                       total_B = total_B + B;
                   end
                   current_particle.ActingBField = total_B;
   
                   % 根据磁场更新磁矩（示例逻辑，替换为你的物理模型）
                   current_particle.magneticmoment = current_particle.magneticmoment + current_particle.gamma * cross(current_particle.magneticmoment, total_B) * 1e-6;
               end
           end
       end
   end
   

3. 调用示例

   % 创建两个质子
   p1 = Proton([0,0,0], [1,0,0]);
   p2 = Proton([0,0,1], [0,1,0]);
   % 创建粒子组
   group = GroupParticles([p1,p2]);
   % 更新属性
   group.updateParticleProperties();
   % 查看更新后的属性
   disp(p1.ActingBField);
   disp(p1.magneticmoment);
   

三、其他可选的对象关联实现方式 #

如果不想直接在GroupParticles中修改粒子属性，还可以采用以下方式：

• 粒子主动请求数据：在Proton类中添加方法，接收GroupParticles对象作为参数，自行计算并更新自身属性。
• 事件驱动机制：给GroupParticles添加事件（比如ParticlesUpdated），粒子注册该事件，当组内粒子变化时自动触发更新逻辑。
• 中间调度类：创建单独的SimulationManager类，负责协调GroupParticles和Particle之间的属性更新，解耦两者的直接依赖。

内容的提问来源于stack exchange，提问作者Meck3l