使用matlab编写一个4个输入、3个隐藏和1个输出的FNN，使用随机权重进行初始化，sigmoid激活函数，迭代次数3000次的代码

以下是符合您要求的MATLAB代码：

% 定义输入数据矩阵
inputData = [0.1 0.2 0.3 0.4;
             0.5 0.6 0.7 0.8;
             0.9 1.0 1.1 1.2;
             1.3 1.4 1.5 1.6];

% 定义输出数据矩阵
outputData = [0.3 0.4 0.5 0.6];

% 设置模糊神经网络的参数
inputNum = size(inputData, 1);  % 输入层神经元个数
hiddenNum = 3;  % 隐藏层神经元个数
outputNum = 1;  % 输出层神经元个数
epochNum = 3000;  % 训练次数
learnRate = 0.1;  % 学习率

% 初始化模糊神经网络的权重矩阵
w1 = randn(hiddenNum, inputNum);  % 输入层到隐藏层的权重矩阵
w2 = randn(outputNum, hiddenNum);  % 隐藏层到输出层的权重矩阵

% 循环训练模糊神经网络
for epoch = 1:epochNum
    % 前向传播
    hiddenOutput = 1 ./ (1 + exp(-(w1 * inputData)));  % 隐藏层输出
    networkOutput = w2 * hiddenOutput;  % 输出层输出

    % 计算误差
    error = outputData - networkOutput;

    % 反向传播
    delta2 = error;  % 输出层误差
    delta1 = (w2' * delta2) .* hiddenOutput .* (1 - hiddenOutput);  % 隐藏层误差

    % 更新权重矩阵
    w2 = w2 + learnRate * delta2 * hiddenOutput';  % 隐藏层到输出层的权重矩阵
    w1 = w1 + learnRate * delta1 * inputData';  % 输入层到隐藏层的权重矩阵
end

% 测试模糊神经网络
testData = [0.3 0.4 0.5 0.6];
testOutput = w2 * (1 ./ (1 + exp(-(w1 * testData'))));  % 输出层输出
disp(testOutput);

请注意，此代码中的“迭代次数”可能会使神经网络过度拟合，因此可能需要进行调整。此外，为了使测试数据成为一维行向量，我们需要使用`testData'`，而不是`testData`。