版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

zuo19880318

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 322.5
帖子: 112
在线: 8.1小时
虫号: 1395321
注册: 2011-09-08
性别: GG
专业: 计算数学与科学工程计算

[求助] 求高手用MATLAB编写一个画张量积曲面的程序

希望给出16个控制顶点，即4行4列的双三次Bezier网格，画出双三次Bezier曲面。

回复此楼

» 猜你喜欢

职称评审没过，求安慰已经有49人回复
26申博自荐已经有3人回复
A期刊撤稿已经有4人回复
垃圾破二本职称评审标准已经有17人回复
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿已经有22人回复
EST投稿状态问题已经有7人回复
毕业后当辅导员了，天天各种学生超烦已经有4人回复
三无产品还有机会吗已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

关于Matlab画三维曲面图的问题，求助高手，感激感激已经有3人回复

1楼 2013-03-29 18:55:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

wurongjun

专家顾问 (职业作家)

专家经验: +831
应助: 791 (博后)
贵宾: 0.308
金币: 24609
散金: 310
红花: 75
帖子: 3004
在线: 881.2小时
虫号: 1368482
注册: 2011-08-14
性别: GG
专业: 计算数学与科学工程计算
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zuo19880318: 金币+15, ★★★很有帮助 2013-03-30 19:53:46

function bzqum33
clear;
N=30;M=30;
x=[1 2 3 4];y=[1 2 3 4];
n=length(x)-1;
[x y]=meshgrid(x,y);
z=[1 2 2 1;
2 2.5 2.5 2;
2 2.5 2.5 2;
1 2 2 1 ];
t=linspace(0,1,N);
s=linspace(0,1,M);
xst=zeros(M,N);yst=zeros(M,N);zst=zeros(M,N);
xst=T(n,s)'*x*T(n,t);
yst=T(n,s)'*y*T(n,t);
zst=T(n,s)'*z*T(n,t);
axis equal;
axis([0 6 0 6 1 3]);
for i=1:4
plot3(x(i,

,y(i,

,z(i,

,'-or','LineWidth',2);
hold on;
plot3(x(:,i),y(:,i),z(:,i),'-ob','LineWidth',2);
end
axis([0 6 0 6 1 3]);
grid on
mesh(xst,yst,zst);
colormap([1,0,0]);
hold on;
xlabel('X-axis');ylabel('Y-axis');zlabel('Z-axis');
hold off;
function y=T(n,t)
y=zeros(n,length(t));
for i=1:n+1
y(i,

=nchoosek(n,i-1)*(1-t).^(n+1-i).*t.^(i-1);
end

qm.jpg

赞一下(2人)

回复此楼

善恶到头终有报,人间正道是沧桑.

2楼2013-03-30 18:12:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zuo19880318

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 322.5
帖子: 112
在线: 8.1小时
虫号: 1395321
注册: 2011-09-08
性别: GG
专业: 计算数学与科学工程计算

引用回帖:

2楼: Originally posted by wurongjun at 2013-03-30 18:12:10
function bzqum33
clear;
N=30;M=30;
x=;y=;
n=length(x)-1;
=meshgrid(x,y);
z=;
t=linspace(0,1,N);
s=linspace(0,1,M);
xst=zeros(M,N);yst=zeros(M,N);zst=zeros(M,N);
xst=T(n,s)'*x*T(n,t);
yst=T ...

请问在此基础上做渐进迭代逼近（PIA），你会嘛？

赞一下(1人)

回复此楼

3楼2013-03-30 19:53:24

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

wurongjun

专家顾问 (职业作家)

专家经验: +831
应助: 791 (博后)
贵宾: 0.308
金币: 24609
散金: 310
红花: 75
帖子: 3004
在线: 881.2小时
虫号: 1368482
注册: 2011-08-14
性别: GG
专业: 计算数学与科学工程计算
管辖: 数学

好吧!金币要给的!!
好的!
function bzqumnn
clear;
N=30;M=30;%网格大小
x=[1 2 3 4;
1 2 3 4;
1 2 3 4;
1 2 3 4];
y=[1 1 1 1;
2 2 2 2;
3 3 3 3;
4 4 4 5];
z=[1 2 2 1;
2 2.5 2.5 2;
2 2.5 2.5 2;
1 2 2 1 ];%十六个点的坐标,注意排成网格状
[xst,yst,zst]=bezqum(x,y,z,M,N,3);%调用函数
axis equal;
mesh(xst,yst,zst);
xlabel('X-axis');ylabel('Y-axis');zlabel('Z-axis');
function [xst,yst,zst]=bezqum(x,y,z,M,N,n)
t=linspace(0,1,N);
s=linspace(0,1,M);
xst=zeros(M,N);yst=zeros(M,N);zst=zeros(M,N);
xst=T(n,s)'*x*T(n,t);
yst=T(n,s)'*y*T(n,t);
zst=T(n,s)'*z*T(n,t);
function y=T(n,t)
y=zeros(n,length(t));
for i=1:n+1
y(i,

=nchoosek(n,i-1)*(1-t).^(n+1-i).*t.^(i-1);
end

qm.jpg

赞一下(1人)

回复此楼

善恶到头终有报,人间正道是沧桑.

8楼2013-03-30 22:35:36

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

zuo19880318

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 322.5
帖子: 112
在线: 8.1小时
虫号: 1395321
注册: 2011-09-08
性别: GG
专业: 计算数学与科学工程计算

引用回帖:

8楼: Originally posted by wurongjun at 2013-03-30 22:35:36
好吧!金币要给的!!
好的!
function bzqumnn
clear;
N=30;M=30;%网格大小
x=;
y=;
z=;%十六个点的坐标,注意排成网格状
=bezqum(x,y,z,M,N,3);%调用函数
axis equal;
mesh(xst,yst,zst);
xlabel('X-axis') ...

昨天金币全给你了。已经40个给完了，现在给不了。。。

赞一下(1人)

回复此楼

9楼2013-03-31 12:07:04

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖