2014-09-25 108 views
0

我有这3个矢量,我需要它们的3D图形。
X和Y是坐标,Z是能量。 这个问题可能可以用griddata或interp2解决 - 但我真的不明白如何处理它们。3个矢量Matlab 3d图

x = [703 704 705 706 762 760 762 764 766 766 771 771 ... 
    828 831 502 836 891 892 895 318 320 379 379 382 383 ... 
    384 442 443 444 444 445 447 449 511 508 511 509 509 ... 
    510 509 510 574 573 573 573 573 573 574 576 575 578 ... 
    635 635 634 637 637 639 641 642 643 827 827 828 644 ... 
    697 697 891 698 699 892 704]; 

y = [498 563 629 692 179 241 305 371 435 498 561 624 ... 
    429 496 176 559 374 430 565 502 565 374 439 505 569 ... 
    631 183 312 375 440 501 567 631 247 309 373 438 502 ... 
    566 628 694 181 244 308 372 437 502 565 629 693 755 ... 
    181 245 310 372 436 500 564 629 693 244 308 371 755 ... 
    180 244 308 308 372 503 436]; 

z = [71172.63 59691.115 41209.71 45061.89 59223.98 ... 
    42330.585 41934.85 66172.105 45911.76 45078.01 68746.4 ... 
    51346.567 57932.58 48625.975 42202 44483.783 62927.595 ... 
    53761.958 43314.128 42754.26 31610.293 45376.275 34661.647 ... 
    35105.147 39571.948 38936.785 83547.518 39633.85 44048.802 ... 
    35590.97 49392.215 38008.345 42108.543 47486 26103 28599 ... 
    31419 29271.158 41591.745 32203.217 35669.092 50788.315 ... 
    29967.715 30581.138 30694.737 33782.143 40074.855 28531.265 ... 
    38337.658 31429.015 33893.448 58948.058 39229.93 39465.64 ... 
    37449.57 40265.768 43811.802 42839.647 50270.707 55402.198 ... 
    54075.66 43573.7325 38312.1075 57198.402 58962.967 50880.345 ... 
    52337.115 59699.333 43770.645 55210.7025 77166.775]; 
+1

嘿!如果它适用于您接受答案,那是多么StackOerflow工作! ;) – 2014-09-26 11:28:43

回答

4

您可以创建一个3D画面与插网:

x=[703 704 705 706 762 760 762 764 766 766 771 771 828 831 502 836 891 892 895 318 320 379 379 382 383 384 442 443 444 444 445 447 449 511 508 511 509 509 510 509 510 574 573 573 573 573 573 574 576 575 578 635 635 634 637 637 639 641 642 643 827 827 828 644 697 697 891 698 699 892 704] 
y=[498 563 629 692 179 241 305 371 435 498 561 624 429 496 176 559 374 430 565 502 565 374 439 505 569 631 183 312 375 440 501 567 631 247 309 373 438 502 566 628 694 181 244 308 372 437 502 565 629 693 755 181 245 310 372 436 500 564 629 693 244 308 371 755 180 244 308 308 372 503 436] 
z=[71172.63 59691.115 41209.71 45061.89 59223.98 42330.585 41934.85 66172.105 45911.76 45078.01 68746.4 51346.567 57932.58 48625.975 42202 44483.783 62927.595 53761.958 43314.128 42754.26 31610.293 45376.275 34661.647 35105.147 39571.948 38936.785 83547.518 39633.85 44048.802 35590.97 49392.215 38008.345 42108.543 47486 26103 28599 31419 29271.158 41591.745 32203.217 35669.092 50788.315 29967.715 30581.138 30694.737 33782.143 40074.855 28531.265 38337.658 31429.015 33893.448 58948.058 39229.93 39465.64 37449.57 40265.768 43811.802 42839.647 50270.707 55402.198 54075.66 43573.7325 38312.1075 57198.402 58962.967 50880.345 52337.115 59699.333 43770.645 55210.7025 77166.775] 

% interpolation on regular grid 
xlin=linspace(min(x),max(x),50); 
ylin=linspace(min(y),max(y),50); 
[X,Y]=meshgrid(xlin,ylin); 
Z=griddata(x,y,z,X,Y,'cubic'); 

% visualization 
mesh(X,Y,Z); 
axis tight; hold on 
plot3(x,y,z,'.', 'MarkerSize',15) 
xlabel('x') 
ylabel('y') 
%surf(X,Y,Z) 

结果是这样的: enter image description here

的代码进行详细道格对他的site解释。
当然,您可以用'线性'替换'立方'插值..

+0

非常感谢你 – 2014-09-25 20:05:18

+0

+1 - 这是我会做的方法。 – rayryeng 2014-09-25 21:26:39

+0

很酷,超级棒!我需要它,因为我正在进行一些数据分析,但我没有太多时间花时间研究matlab的工作原理,所以这是一项了不起的工作! Thnx再次:) – 2014-09-28 01:26:12