曲线程序部分已经编好,但是导出为txt部分一直导出失败,图形如下
自己摸索修改的转化为txt文件的程序:
c = (xe-xc)/0.01+1;
c1 = (xd-xb)/0.01+1;
fix(c)
fix(c1)
z1 = zeros(1,fix(c));
z2 = zeros(1,fix(c1));
%figure
%plot(s1,y1,s2,y2);
a = [xt1;yt1;z1]';
fid1=fopen(['d:','curve2.txt'],'w');
[l,m,n]=size(a);%利用size函数读取a矩阵的行,列,维
for i=1:l
for j=1:m
for k=1:n
fprintf(fid1,'%f\t',a(i,j,k));%fprintf?将a矩阵的数据以指定的格式写入文本中
end
end
fprintf(fid1,'\r\n');
end
fclose(fid1);
b = [xt2;yt2;z2]';
fid2=fopen(['d:','curve3.txt'],'w');
[l,m,n]=size(b);
for i=1:l
for j=1:m
for k=1:n
fprintf(fid2,'%f\t',b(i,j,k));
end
end
fprintf(fid2,'\r\n');
end
format long
两个弧线为贝塞尔曲线,参数方程。
两弧线分方程为
xt1=@(t) xc*(1-t)^3+3xm1t*(1-t)^2+3xm2t^2*(1-t)+xet^3
yt1=@(t) yc(1-t)^3+3ym1t*(1-t)^2+3ym2t^2*(1-t)+yet^3
xt2=@(t) xb(1-t)^3+3xn1t*(1-t)^2+3xn2t^2*(1-t)+xdt^3
yt2=@(t) yb(1-t)^3+3yn1t*(1-t)^2+3yn2t^2*(1-t)+yd*t^3
运行结果报错如下
真诚求帮助 求这一段的程序