clc;
clear all
x=[0 10 20];
y=[0 10 20 30];
z=[1 2 5 3;4 1 0 6;6 3 4 3];
x0=x';
y0=y;
R=x0(:,ones(1,length(y)));
G=y0(ones(1,length(x)),:);
B=z;
for i=1:2
N=length(x);
M=length(y);
alpha=0.3*ones(N,M);
xx=0;
yy=0;
s=1;
s1=1;
s2=1;
s3=0;
s4=0;
for i=1:N-1
for j=1:N
xx=x(i)+(x(i+1)-x(i))*(x(j)-x(1))/(x(N)-x(1));
xp(s1)=xx;s1=s1+1;
end
end
for i=1:M-1
for j=1:M
yy=y(i)+(y(i+1)-y(i))*(y(j)-y(1))/(y(M)-y(1));
yp(s2)=yy;s2=s2+1;
end
end
KK=[];
for i=2:N
for j=2:M
c=(z(i-1,j-1)-z(i-1,j)-z(i,j-1)+z(i,j)-alpha(i,j)*(z(1,1)-z(N,1)-z(1,M)+z(N,M)))/(x(1)*y(1)-x(N)*y(1)-x(1)*y(M)+x(N)*y(M));
b=(z(i-1,j-1)-z(i,j-1)-alpha(i,j)*(z(1,1)-z(N,1))-c*(x(1)*y(1)-x(N)*y(1)))/(x(1)-x(N));
d=(z(i-1,j-1)-z(i-1,j)-alpha(i,j)*(z(1,1)-z(1,M))-c*(x(1)*y(1)-x(1)*y(M)))/(y(1)-y(M));
k=z(i,j)-b*x(N)-d*y(M)-alpha(i,j)*z(N,M)-c*x(N)*y(M);
for t1=1:N
for t2=1:M
KK=[KK;[(i-2)*M+t1,(j-2)*N+t2]];
zz((i-2)*N+t1,(j-2)*M+t2)=b*x(t1)+d*y(t2)+c*x(t1)*y(t2)+alpha(i,j)*z(t1,t2)+k;
end
end
end
end
x=xp;
y=yp;
z=zz;
end
[xx1 yy1]=meshgrid(yp,xp);
plot3(yy1,xx1,zz,'k.');
hold on
plot3(R, G, B, 'kO');
axis tight
没有合适的资源?快使用搜索试试~ 我知道了~
fractalsurfacegenerater.zip_matlab 分形地形_三维表面特征_分形三维_地形 matlab_地形
共1个文件
m:1个
1.该资源内容由用户上传,如若侵权请联系客服进行举报
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
2.虚拟产品一经售出概不退款(资源遇到问题,请及时私信上传者)
版权申诉
5星 · 超过95%的资源 2 下载量 18 浏览量
2022-07-15
04:22:42
上传
评论 1
收藏 735B ZIP 举报
温馨提示
用于产生三维的、具有分形特征的地形表面或者具有相关性质的三维模型。
资源推荐
资源详情
资源评论
收起资源包目录
fractalsurfacegenerater.zip (1个子文件)
fractalsurfacegenerater.m 1KB
共 1 条
- 1
资源评论
- 楓淩2023-07-19资源质量不错,和资源描述一致,内容详细,对我很有用。
- m0_750287012023-02-25这个资源对我启发很大,受益匪浅,学到了很多,谢谢分享~
小波思基
- 粉丝: 72
- 资源: 1万+
下载权益
C知道特权
VIP文章
课程特权
开通VIP
上传资源 快速赚钱
- 我的内容管理 展开
- 我的资源 快来上传第一个资源
- 我的收益 登录查看自己的收益
- 我的积分 登录查看自己的积分
- 我的C币 登录后查看C币余额
- 我的收藏
- 我的下载
- 下载帮助
安全验证
文档复制为VIP权益,开通VIP直接复制
信息提交成功