ceshi = [ceshi delta]; %向其中连续追加数据,得到1xn的矩阵
% ceshi = [ceshi;delta];%得到nx1的矩阵
end
end
%加入三维曲面绘制程序
[X,Y]=meshgrid(xita_array,fin_array);
v1 = reshape(ceshi, 41, 41);
meshz(X,Y,v1); view(3);
73、两个等直径圆管的交线
clear
m=30;
z=1.2*(0:m)/m;
r=ones(size(z));
theta=(0:m)*2*pi/m;
x1=r'*cos(theta);y1=r'*sin(theta);
z1=z'*ones(1,m+1);
x=(-m:2:m)/m;
x2=x'*ones(1,m+1);y2=r'*cos(theta);
z2=r'*sin(theta);
surf(x1,y1,z1);
axis equal,axis off
hold on
surf(x2,y2,z2);
axis equal ,axis off
title('两个等直径圆管的交线');
hold off
74、Matlab的csvread读取数据的问题
u0=csvread('data.csv')不好用,因为我的文件第一行是文本,好像csvread要求所读取的csv文件必须都是数据。什么语句能删除csv文件里面的第一行,然后读取单元格A3到A10的数据?(假设csv里面有20*20的数据)
方法:
用如下方式读取:
csvread('D:\chengji.txt',1,0); %把第一行的内容排除
代码: txt文件的内容
liu1,liu2,liu3
80,75,65
60,65,35
代码:
>> csvread('D:\chengji.txt',1,0)
ans =
80 75 65
60 65 35
75、在对数坐标系插值的问题
我在双对数坐标系(loglog)中有一组坐标点[x0 y0],请问如何在些点之间做直线插值,让插值得到的点落在坐标之间的直线上。
由于对数坐标系的原因,我用 yi = interp1(x0, y0, xi) 得到的点 (xi yi) 都不在直线上,如下图所示,应如何让这些插值点落在直线上?