Cómo utilizar MATLAB para dibujar una superficie 3D en cualquier dato 3D

Método: aplique el método de ajuste a los dos conjuntos de datos para obtener la ecuación de regresión para el tercer conjunto de datos y expanda los puntos de datos a través de la ecuación de regresión y la cuadrícula.

Herramienta: Matlab2011a

Ejemplo:

(1) Columna de datos:

K/NM/° C/NMS/° At /NM

10 0,1 128

20 0,2 260

30 0,3 393

40 0,4 531

50 0,5 672

60 0,6 815

70 0,7 954

80,27 0,8 1102

90 0,9 1251

100 1 1406

(2) Dibujar la superficie tridimensional entre K_C_At

(3) Pasos:

1) Importar datos directamente en la ventana de comandos de matlab

2) Escribe m. Archivos

Los archivos fuente son los siguientes:

************ 3D*************** *

Datos de entrada

Superficie tridimensional entre rigidez K, amortiguación C y diferencia de pico del resorte de torsión Pt

x1: rigidez de entrada K, unidad NM/°

x2: Amortiguación de entrada C, unidad NMS/°

y: Diferencia de amplitud del resorte de torsión Pt

******** ******** **************** ****************

x1=datos(: ,1);c columna de datos

x2=data(：,2);k columna de datos

y=data(：,3);En la columna de datos

****** ************************* *************

Regresión lineal dividida múltiple Convertir a ecuación de regresión lineal para resolver

X=[ ones(size(x1)) x1.*x1 x1 x2.*x2 x2 x1.*x2];

[b, bint, r, rint, stats]=regress(y, X, 0.05);

b

Distribución del diagrama de dispersión 3D

scatter3(x1, 'lleno')

espera;

*************************** ********* ****************

Diferencia de datos

x1fit=10:1:109;

x2fit=0.1 :0.009:1;

****************************** ************ **********

[xi, yi]= meshgrid(x1fit, x2fit);

Ajuste no lineal ecuación

zi=b(1 ) b(2)*xi.*xi b(3)*xi b(4)*yi.*yi b(5)*yi b(6)*xi .*yi;

Dibujar superficie tridimensional

meshc (xi, yi, zi)

Mantenga presionado

; Ajustar la vista

vista(140, 30)

*************************** ********************* ****

Agregar coordenadas

xlabel('rigidez/(

NM/°)');ylabel('amortiguación/(NMS/°)');

zlabel('amplitud de aberración con par/(NM)'); (2)

Análisis de residuos

rcoplot(r, rint

);