帮我写一段matlab代码，利用sigmoid实现与门感知器

clear all;
clc;

%定义输入矩阵
X = [0 0; 0 1; 1 0; 1 1];

%定义目标输出
Y = [0; 0; 0; 1];

%初始化权重矩阵和偏置
W = rand(2,1);
b = rand();

%定义学习率和迭代次数
alpha = 0.1;
epochs = 10000;

%使用sigmoid函数实现与门感知器
for i = 1:epochs
%前向传播
Z = X*W + b;
A = 1./(1+exp(-Z));

%计算误差
error = Y - A;

%反向传播
dW = X'*error.*A.*(1-A);
db = sum(error.*A.*(1-A));

%更新权重和偏置
W = W + alpha*dW;
b = b + alpha*db;
end

%输出结果
disp('权重矩阵：');
disp(W);
disp('偏置：');
disp(b);

%测试模型
test_input = [0 0; 0 1; 1 0; 1 1];
test_output = 1./(1+exp(-test_input*W-b));
disp('测试结果：');
disp(test_output);