声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: xjzuo

[编程技巧] 讨论: 关于如何画出自定义大小的网格

[复制链接]
发表于 2007-8-6 16:45 | 显示全部楼层
like this?

function myrectangularmesh
%
clear all
x=0:10;y=0:10;
X0=1;Y0=1;
n1=0;
z=zeros(11,11);
for i=0:10;
    for j=0:10
        n1=n1+1;
        x(n1)=i*X0;
        y(n1)=j*Y0;
        z(i+1,j+1)=n1;
    end
end
ndm=n1;
for i=1:11
    for j=1:11
        L=0;a=0;
        k5=z(i,j);
        if i-1>0
           L=L+1;
           k1=z(i-1,j);
           AA(k5,k1)=-1;
        end
        if j-1>0
           L=L+1;
           k2=z(i,j-1);
           AA(k5,k2)=-1;
        end
        if i+1<=11
           L=L+1;
           k3=z(i+1,j);
           AA(k5,k3)=-1;
        end
        if j+1<=11
           L=L+1;
           k4=z(i,j+1);
           AA(k5,k4)=-1;
        end
        AA(k5,k5)=L;
    end
end
hold on
AA=sparse(AA);
BB=full(AA);
   figure(1);
   xy=[x' y'];
   size(xy);
   gplot(AA,xy)
   axis equal,axis tight
   set(gca,'xtick',[0:10])
   set(gca,'xticklabel',{'0';'1';'2';'3';'4';'5';'6';'7';'8';'9';'10'})
回复 支持 反对
分享到:

使用道具 举报

发表于 2010-8-26 14:56 | 显示全部楼层
难道楼主是要这样?
x = (-2:.2:2);
y = (-1.5:.2:1.5)';
X = x(ones(size(y)),:);
Y = y(:,ones(size(x)));
F = X.*exp(-X.^2-Y.^2);
它与
[X,Y] = meshgrid(-2:.2:2, -1.5:.2:1.5);
F = X .* exp(-X.^2 - Y.^2);
是等价的
发表于 2010-8-26 15:13 | 显示全部楼层
若要用类C语言编写,楼主看这样行不?

function [X,Y]=mesh_grid(x,y)

x=1:3;y=2:4;
xLen=length(x);yLen=length(y);
X=zeros(yLen,xLen);Y=X;
for b=1:xLen
    X(:,b)=x(b);
end
for a=1:yLen
    Y(a,:)=y(a);
end
发表于 2011-4-11 15:50 | 显示全部楼层
如何在三角形区域划分三角形网格呢?请高手指点一下!
发表于 2014-4-4 22:10 | 显示全部楼层
学习了  感谢大家分享
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-25 03:03 , Processed in 0.051012 second(s), 16 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表