询问MATLAB程序cleark=3;x=xlsread('2000.xls')[n,m]=size(x);%-----------弦长参数化--------------------------------------u(k+n)=0;for i=1:n-1u(k+i+1)=u(k+i)+sqrt((x(i+1,1)-x(i,1))^2+(x(i+1,2)-x(i,2))^2);end;L=u(n+k);for i=1:nu(k+i)=u(k+i)/L;end
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/03 23:12:58
![询问MATLAB程序cleark=3;x=xlsread('2000.xls')[n,m]=size(x);%-----------弦长参数化--------------------------------------u(k+n)=0;for i=1:n-1u(k+i+1)=u(k+i)+sqrt((x(i+1,1)-x(i,1))^2+(x(i+1,2)-x(i,2))^2);end;L=u(n+k);for i=1:nu(k+i)=u(k+i)/L;end](/uploads/image/z/10212908-68-8.jpg?t=%E8%AF%A2%E9%97%AEMATLAB%E7%A8%8B%E5%BA%8Fcleark%3D3%3Bx%3Dxlsread%28%272000.xls%27%29%5Bn%2Cm%5D%3Dsize%28x%29%3B%25-----------%E5%BC%A6%E9%95%BF%E5%8F%82%E6%95%B0%E5%8C%96--------------------------------------u%28k%2Bn%29%3D0%3Bfor+i%3D1%3An-1u%28k%2Bi%2B1%29%3Du%28k%2Bi%29%2Bsqrt%28%28x%28i%2B1%2C1%29-x%28i%2C1%29%29%5E2%2B%28x%28i%2B1%2C2%29-x%28i%2C2%29%29%5E2%29%3Bend%3BL%3Du%28n%2Bk%29%3Bfor+i%3D1%3Anu%28k%2Bi%29%3Du%28k%2Bi%29%2FL%3Bend)
询问MATLAB程序cleark=3;x=xlsread('2000.xls')[n,m]=size(x);%-----------弦长参数化--------------------------------------u(k+n)=0;for i=1:n-1u(k+i+1)=u(k+i)+sqrt((x(i+1,1)-x(i,1))^2+(x(i+1,2)-x(i,2))^2);end;L=u(n+k);for i=1:nu(k+i)=u(k+i)/L;end
询问MATLAB程序
clear
k=3;
x=xlsread('2000.xls')
[n,m]=size(x);
%-----------弦长参数化--------------------------------------
u(k+n)=0;
for i=1:n-1
u(k+i+1)=u(k+i)+sqrt((x(i+1,1)-x(i,1))^2+(x(i+1,2)-x(i,2))^2);
end;
L=u(n+k);
for i=1:n
u(k+i)=u(k+i)/L;
end;
for i=1:3
u(k+i+n)=1;
end
%控制多边线
plot(x(:,1),x(:,2),'o');
hold on
A=zeros(n+2);
A(1,1)=1;A(1,2)=-1;
A(2,2)=1;
A(n+2,n+1)=-1;A(n+2,n+2)=1;
A(n+1,n+1)=1;
for i=3:n
for j=0:2
A(i,i+j-1)=Base(i+j-1,k,u,u(i+2));
end
end
%e:方程右边.
e=0;
for i=1:m
e(n+2,i)=0;
end
for i=1:n
e(i+1,:)=x(i,:);
end
%求出控制点d
d=inv(A)*e;
plot(d(:,1),d(:,2),'g');
hold on
错误提示为:
Warning:Matrix is singular to working precision.
xls中的数据为110*109的矩阵
询问MATLAB程序cleark=3;x=xlsread('2000.xls')[n,m]=size(x);%-----------弦长参数化--------------------------------------u(k+n)=0;for i=1:n-1u(k+i+1)=u(k+i)+sqrt((x(i+1,1)-x(i,1))^2+(x(i+1,2)-x(i,2))^2);end;L=u(n+k);for i=1:nu(k+i)=u(k+i)/L;end
在求矩阵的逆时,必须保证矩阵非奇异,即det(A)不等于零.现在提示的错误是矩阵是奇异的,即det(A)等于零.这不是程序问题,请检查你的数据.