数学建模案例三

目录

ps：如果你觉得对你有帮助的话，麻烦点个赞，加个关注

• 绘制三维空间曲线：
• 绘制三维空间曲面：
• 绘制等高线：

绘制三维空间曲线：

做出三维螺旋曲线

t=0:pi/50:10*pi;
plot3(sin(t),cos(t),t);

绘制三维空间曲面：

做出二元函数f=x^2 + y^2的图形

clf
x=-4:4; y=x;
[X,Y] = meshgrid(x,y);
Z = X.^2+Y.^2;
surf(X,Y,Z);

用来表现在格点上计算函数值

clf
x=-4:4; y=x;
[X,Y] = meshgrid(x,y);
Z = X.^2+Y.^2;
surf(X,Y,Z);
hold on,
colormap(hot)
stem3(X,Y,Z,'bo') %用来表现在格点上计算函数值
hold off

例题

clc;
clear;
x=0: 400: 5600;
y=0: 400: 4800;
z=[370 470 550 600 670 690 670 620 580 450 400 300 100 150 250;
510 620 730 800 850 870 850 780 720 650 500 200 300 350 320; 
650 760 880 970 1020 1050 1020 830 900 700 300 500 550 480 350;
740 880 1080 1130 1250 1280 1230 1040 900 500 700 780 750 650 550;
830 980 1180 1320 1450 1420 1400 1300 700 900 850 840 380 780 750;
880 1060 1230 1390 1500 1500 1400 900 1100 1060 950 870 900 930 950;
910 1090 1270 1500 1200 1100 1350 1450 1200 1150 1010 880 1000 1050 1100;
950 1190 1370 1500 1200 1100 1550 1600 1550 1380 1070 900 1050 1150 1200;
1430 1430 1460 1500 1550 1600 1550 1600 1600 1600 1550 1500 1500 1550 1550;
1420 1430 1450 1480 1500 1550 1510 1430 1300 1200 980 850 750 550 500;
1380 1410 1430 1450 1470 1320 1280 1200 1080 940 780 620 460 370 350;
1370 1390 1410 1430 1440 1140 1110 1050 950 820 690 540 380 300 210;
1350 1370 1390 1400 1410 960 940 880 800 690 570 430 290 210 150];
meshz(x,y,z); rotate3d,
figure(2)
surf(x,y,z), rotate3d;

冷色调

clc;
clear;
x=0: 400: 5600;
y=0: 400: 4800;
z=[370 470 550 600 670 690 670 620 580 450 400 300 100 150 250;
510 620 730 800 850 870 850 780 720 650 500 200 300 350 320; 
650 760 880 970 1020 1050 1020 830 900 700 300 500 550 480 350;
740 880 1080 1130 1250 1280 1230 1040 900 500 700 780 750 650 550;
830 980 1180 1320 1450 1420 1400 1300 700 900 850 840 380 780 750;
880 1060 1230 1390 1500 1500 1400 900 1100 1060 950 870 900 930 950;
910 1090 1270 1500 1200 1100 1350 1450 1200 1150 1010 880 1000 1050 1100;
950 1190 1370 1500 1200 1100 1550 1600 1550 1380 1070 900 1050 1150 1200;
1430 1430 1460 1500 1550 1600 1550 1600 1600 1600 1550 1500 1500 1550 1550;
1420 1430 1450 1480 1500 1550 1510 1430 1300 1200 980 850 750 550 500;
1380 1410 1430 1450 1470 1320 1280 1200 1080 940 780 620 460 370 350;
1370 1390 1410 1430 1440 1140 1110 1050 950 820 690 540 380 300 210;
1350 1370 1390 1400 1410 960 940 880 800 690 570 430 290 210 150];
meshz(x,y,z); rotate3d,
figure(2)
shading interp;
colormap(cool)
colorbar
surf(x,y,z), rotate3d;

绘制等高线：

[X,Y,Z]=peaks(30);surf(X,Y,Z)
figure(2);
[C,h] = contour(X,Y,Z,16);
clabel(C,h)
figure(3);
contour(X,Y,Z,16);