在matlab中绘制一个砂轮在磨削加工中的平面受力图，其中包含磨削力，热应力，夹紧力和离心力，并且凸显出它们在砂轮中的位置，请给我代码

在MATLAB中，绘制砂轮在磨削加工中受力图通常涉及到使用`patch`函数创建矩形区域代表砂轮，然后分别添加各个力作为箭头或者其他形状表示。下面是一个基本的示例，展示如何创建这样的图形，假设我们已经定义了每个力的大小、方向以及砂轮的位置：

% 定义砂轮的中心点和半径
r = 0.1; % 砂轮半径
xc = 0;   % 横向中心点
yc = 0;   % 纵向中心点

% 定义四个力的数据 (这里假设数值为随机数)
force1_mag = rand(1);  % 磨削力
force1_dir_x = r / 4;  % 磨削力的方向偏移量（沿x轴）
force1_dir_y = 0;       % 磨削力垂直于砂轮表面
force2_mag = rand(1);  % 热应力
force2_dir_x = -r / 6; % 热应力反方向偏移
force2_dir_y = 0;
force3_mag = rand(1);  % 夹紧力
force3_dir_x = -r / 8;
force3_dir_y = 0;
force4_mag = rand(1);  % 离心力
force4_dir_x = 0;
force4_dir_y = r * 0.5; % 远离砂轮中心

% 创建矩形区域代表砂轮
[x, y] = meshgrid(-r:r, -r:r);
figure;
hold on;

% 绘制砂轮
patch(x(:), y(:), 'FaceColor', 'w', 'EdgeColor', 'k');
text(xc, yc, '砂轮');

% 添加受力箭头
% 使用quiver函数绘制矢量场
quiver([xc force1_dir_x + xc], [yc force1_dir_y + yc], [0 force1_mag], 'color', 'red', 'LineWidth', 1.5);
text(xc + force1_dir_x + 0.05, yc + force1_dir_y + 0.1, '磨削力');

% 其他力同理，你可以添加更多quiver指令，对应各自的力方向和大小

hold off;
xlabel('X方向');
ylabel('Y方向');
title('砂轮磨削受力图');