2014-09-25 154 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