Code an equation in MATLAB

2 次查看（过去 30 天）

显示更早的评论

HN 2022-5-7

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1713505-code-an-equation-in-matlab

评论： HN 2022-5-18

在 MATLAB Online 中打开

In matlab, I write the following code to the above equation as follows

ri = sqrt(psi_coordinates.^2+theta_coordinates.^2);
A = sum((ri*delta_epsilon)/2);
V = delta_z*sum(A); 

But the result I obtained is not correct. Can anyone comment if I made some mistake please ?

Thank you .

27 个评论
显示 25更早的评论隐藏 25更早的评论

HN 2022-5-10

编辑：HN 2022-5-17

在 MATLAB Online 中打开

I didn’t get it. It is from a paper.

Below is the part of program @Dyuman Joshi

while point(3) <= z_end + 0.000001
    work_points = zeros(2,1);
    z_vector = [z_vector point(3)];
    delta_meridian = 2*pi/meridian_number;
    for meridian = 0 : delta_meridian : 2*pi
        point(4) = 0;				
        point(5) = 0;				
        
        	if point(3) < 0.2*leg_length					
        		 delta_angular = delta_start / 200;			
        	elseif point(3) >= 0.2*leg_length && point(3) <= 0.3*leg_length	
        		 delta_angular =  delta_start / 50;		
        	else								
        		 delta_angular = delta_start;
            end
        
       %delta_angular = delta_start;		
        
        ccontinue = 1;				
        point(4) = point(4) + delta_angular * cos(meridian);
        point(5) = point(5) + delta_angular * sin(meridian);
        while  ccontinue == 1
            ccontinue = workspace(constrain(point),max_cond);
            if ccontinue == 0 && delta_angular > tolerance		
                point(4) = point(4) - delta_angular * cos(meridian);
                point(5) = point(5) - delta_angular * sin(meridian);
                delta_angular = delta_angular / 2;			
                ccontinue = 1;					
                
            end
            psi_temp = point(4); theta_temp = point(5);
            point(4) = point(4) + delta_angular * cos(meridian);
            point(5) = point(5) + delta_angular * sin(meridian);
        end
        temp_vect = [psi_temp;theta_temp];
        work_points = [work_points temp_vect];
    end
    work_surf = [work_surf' work_points']';
    point(3) = point(3) + delta_z;
    
    
    
end
size_temp = size(work_surf);
work_surf = work_surf(:,2:size_temp(2));
size_temp = size(work_surf);
psi_coordinates = []; theta_coordinates = [];
for i = 1 :2:(z_slices*2+2)			
    psi_coordinates = [psi_coordinates work_surf(i,:)'];
    theta_coordinates = [theta_coordinates work_surf(i+1,:)'];
end
psi_coordinates = psi_coordinates';
theta_coordinates = theta_coordinates';
 psi = psi_coordinates'; 
 th = theta_coordinates';
z_position = ones(z_slices+1,size_temp(2));	
for i=1:z_slices+1
    z_position(i,:) = z_vector(i);
end
ri = sqrt(psi.^2 + th.^2); 
A = sum((ri*delta_meridian*delta_z)/2,2);
v = delta_z*sum(A)