数据处理

score=fix(rand(7,5)*40+60);sunmlie=sum(score);sumhang=sum(score‘);meanlie=mean(score);meanhang=mean(score‘);a=meanhang‘;stdlie=std(score);stdhang=std(score‘);youliang=score>=80;youlianghang=sum(youliang‘)/5;youlianglie=sum(youliang)/7;

 1 // Simulation.h: interface for the CSimulation class. 2 // 3 ////////////////////////////////////////////////////////////////////// 4  5 #if !defined(SIMULATION_H) 6 #define SIMULATION_H 7  8 #include <vector> 9 #include "Rocket.h"10 #include "Spark.h"11 using namespace std;12 13 class CSimulation  14 {15 public:16 //Member functions:17     //Constructor18     CSimulation();19     //Deconstructor20     virtual ~CSimulation();21     //Display all the particles22     void Display();23     //Calculate total number of particles24     int GetParticleNum();25     //One Euler step of the simulation26     void EulerStep();27     //When a rocket reaches its top height, it explodes and calls this function to generate sparks28     void Explode(float posx, float posy, float RocketSpeed, float* color);29     //When a key is pressed, Simulation calls this function to generate a rocket30     void FireRocket(float pos, float* color);31     //Set the speed of the wind32     void SetWindSpeed(float windSpeedValue);33     //Set Total Explosions34     void SetTotalExplosions(int explosions);35     36 //Member variables:37     //Vector of CRocket. Rockets in this vector are either still flying or generating sparks38     vector<CRocket*> rockets;39     //Vector of CSpark.40     vector<CSpark*> sparks;41     //Delta time for a Euler step.42     float deltaT;43     //Velocity of the wind in x direction44     float windSpeed;45     unsigned int totalExplosions;46 };47 48 #endif // !defined(SIMULATION_H)

t=0:100;y1=sin(t);y2=sin(3*t);y3=sin(t)+sin(3*t);y4=t.^2;y5=t.^3;y6=t.^-1;y7=log2(t);y8=log10(t);subplot(1,1,1),plot(t,y1,‘r‘);hold on;plot(t,y2,‘g‘);hold on;plot(t,y3,‘b‘);legend(‘sin(t)‘,‘sin(3t)‘,‘sin(t)+sin(3t)‘);title(‘y1/y2/y3‘);xlabel(‘t‘);ylabel(‘y‘);axis([0 20 -2 2]);grid on;figure(4);subplot(1,1,1),plot(t,y4,‘r‘);hold on;plot(t,y5,‘g‘);legend(‘t^2‘,‘t^3‘);title(‘y4/y5‘);xlabel(‘t‘);ylabel(‘y‘);figure(3);subplot(1,1,1),plot(t,y6,‘b‘);legend(‘1/t‘);title(‘y6‘);xlabel(‘t‘);ylabel(‘y‘);figure(2);subplot(1,1,1),plot(t,y7,‘r‘);hold on;plot(t,y8,‘g‘);legend(‘log2(t)‘,‘log10(t)‘);title(‘y7/y8‘);xlabel(‘t‘);ylabel(‘y‘);figure(1);