matlab画三维曲面出现黑图,Matlab画三维图（一个面）为什么会出现线呢？

close all

clear all

clc

%% 离散曲率值矩阵k，弧长矩阵ss，(包含边界点k=0,s=0)求出对应弧段的弧长-曲率线性方程

%Strain=[0 1.75044E-04 1.53376E-04 1.32575E-04 1.12096E-04 9.16930E-05 7.12955E-05 5.08464E-05 3.01552E-05]%应变

%Dspl=[0 -8.57915E-04 -2.36228E-03 -4.51358E-03 -7.22375E-03 -1.04064E-02 -1.39754E-02 -1.78449E-02 -2.19286E-02]%位移

Data= load('F:\ANSYS working directory\modal2\simulation.txt'); %从数据采集文件中提取

Datas=Data(:,1);

Strain=reshape(Datas,11,4);

i=1:11;s(i)=0.1*(i-1); %弧长：单位m

ss=repmat(s',1,4);

h=0.005;%厚度

k=2*Strain/h;%曲率

l=length(s);

for i=1:1:(l-1)

for j=1:4

a(i,j)=(k(i+1,j)-k(i,j))/(ss(i+1,j)-ss(i,j)); %线性方程的参数

b(i,j)=(k(i,j)*ss(i+1,j)-k(i+1,j)*ss(i,j))/(ss(i+1,j)-ss(i,j));

end

end