3.1作业

作业要求:

通过w(红色臂角度增大)s(红色臂角度减小)d(蓝色臂角度增大)a(蓝色臂角度减小)按键控制机械臂

程序代码:

#include<myhead.h>
#define SER_IP "192.168.126.12"      //服务器ip
#define SER_PORT 8888                //服务器端口号

#define CLI_IP "192.168.244.140"      //客户端IP
#define CLI_PORT 9999                //客户端端口号

int main(int argc, const char *argv[])
{
	//1.创建用于通信的套接字文件描述符
	
	int cfd = socket(AF_INET,SOCK_STREAM,0);
	if(cfd == -1)
	{
		perror("socket error");
		return -1;
	}
	printf("cfd = %d\n",cfd);

	//2.绑定(非必须)
	//填充地址信息结构体
	struct sockaddr_in cin;
	cin.sin_family = AF_INET;
	cin.sin_port = htons(CLI_PORT);
	cin.sin_addr.s_addr = inet_addr(CLI_IP);
	//绑定
	if(bind(cfd,(struct sockaddr*)&cin , sizeof(cin)) == -1)
	{
		perror("bind error");
		return -1;
	}

	printf("bind success\n");
	
	//3.连接服务器
	//3.1填充要连接的服务器地址信息结构体
	
	struct sockaddr_in sin;
	sin.sin_family = AF_INET;               //地址族
	sin.sin_port = htons(SER_PORT);         //端口号
	sin.sin_addr.s_addr = inet_addr(SER_IP); //ip地址


	//3.2连接服务器
	if(connect(cfd ,(struct sockaddr*)&sin,sizeof(sin)) == -1)
	{
		perror("connect error");
		return -1;
	}
	printf("connect success\n");


	//4.数据收发
	char rbuf[5] = {0xff,0x02,0x00,0x00,0xff};
	unsigned char bbuf[5] = {0xff,0x02,0x01,0x00,0xff};

	//发送给服务器当做初始值
	send(cfd,rbuf,sizeof(rbuf),0);
	sleep(1);                 //防止沾包
	send(cfd,bbuf,sizeof(bbuf),0);

	char key = 0;       //接收键盘输入的字符

	while(1)
	{
		system("stty -icanon");
		key = getchar();
		
		//scanf("%c",&key);                //键盘输入一个字符
		//getchar();                     //吸收回车

		switch(key)
		{
			//操作红色臂角度增大,范围不大于90
		case 'w':case 'W':
			{
				rbuf[3]+=1;                //每次操作的角度偏移一度
				if(rbuf[3]>90)
				{
					rbuf[3] = 90;
				}
				send(cfd,rbuf,sizeof(rbuf),0);
			}
			break;
			//操作红色臂角度减小,范围不小于-90
		case 's':case 'S':
			{
				rbuf[3]-=1;                
				if(rbuf[3]<-90)
				{
					rbuf[3] = -90;
				}
				send(cfd,rbuf,sizeof(rbuf),0);
			}
			break;
			//操作蓝色臂角度增大,范围不大于180
		case 'd':case 'D':
			{
				bbuf[3]+=1;               			
				if(bbuf[3]>180)
				{
					bbuf[3] = 180;
				}
				send(cfd,bbuf,sizeof(bbuf),0);
			}
			break;
			//操作蓝色臂角度减小,范围不小于0
		case 'a':case 'A':
			{
				bbuf[3]-=1;              
				if(bbuf[3]<0)
				{
					bbuf[3] = 0;
				}
				send(cfd,bbuf,sizeof(bbuf),0);
			}
			break;
		}


	}
	//5.关闭套接字
	close(cfd);
	return 0;
}

运行结果:

基于UDP的TFTP文件传输

#include <stdio.h>
#include <sys/socket.h>
#include <sys/types.h>
#include <sys/stat.h>
#include <fcntl.h>
#include <netinet/in.h>
#include <string.h>
#include <unistd.h>
#include <stdlib.h>
#include <errno.h>
#include <arpa/inet.h>

#define ERR_LOG(msg) do{\
	perror(msg);\
	printf("%d %s %s\n", __LINE__, __func__, __FILE__);\
}while(0)


#define PORT 	69
#define N 		516

int do_download(int , struct sockaddr_in );
int do_upload(int sfd, struct sockaddr_in sin);


int main(int argc, const char *argv[])
{
	if(argc < 2)
	{
		printf("请输入IP\n");
		return -1;
	}
	//1.创建报式套接字
	int sfd = socket(AF_INET, SOCK_DGRAM, 0);
	if(sfd < 0)
	{
		perror("socket");
		return -1;
	}

	//2.填充服务器ip和端口号
	struct sockaddr_in sin;
	sin.sin_family 		= AF_INET;
	sin.sin_port 		= htons(PORT);
	sin.sin_addr.s_addr = inet_addr(argv[1]);


	char choose;
	while(1)
	{
		system("clear");
		printf("******************\n");
		printf("******1.下载******\n");
		printf("******2.上传******\n");
		printf("******3.退出******\n");
		printf("******************\n");

		scanf("%c", &choose);
		while(getchar()!=10);

		switch(choose)
		{
		case '1':
			//下载功能, todo
			do_download(sfd, sin);
			break;
		case '2':
			//上传功能,todo
			do_upload(sfd, sin);
			break;
		case '3':
			goto END;
			break;
		default:
			printf("输入错误\n");
		}
		printf("请输入任意字符清屏>>>");
		while(getchar()!=10);		
	}
	
	//关闭文件描述符
END:
	close(sfd);
	return 0;
}


//上传功能
int do_upload(int sfd, struct sockaddr_in sin)
{
	char filename[20] = "";
	printf("请输入要上传的文件名>>>");
	fgets(filename, 20, stdin);
	filename[strlen(filename)-1] = 0;

	//判断该文件是否存在
	int fd = open(filename, O_RDONLY);
	if(fd < 0)
	{
		if(errno == ENOENT)
		{
			printf(">>>文件不存在,请重新输入<<<\n");
			return -2;
		}
		else
		{
			ERR_LOG("open");
			return -1;
		}
	}
	
	//发送上传请求
	//上传协议
	char buf[N] = "";
	int size = sprintf(buf, "%c%c%s%c%s%c", 0, 2, filename, 0, "octet", 0);

	if(sendto(sfd, buf, size, 0, (struct sockaddr*)&sin, sizeof(sin))<0)
	{
		ERR_LOG("sendto");
		return -1;
	}

	//循环接收发送数据包
	int recv_len;
	unsigned short num = 0;
	socklen_t addrlen = sizeof(sin);
	while(1)
	{
		bzero(buf, N);
		recv_len = recvfrom(sfd, buf, sizeof(buf), 0, (struct sockaddr*)&sin, &addrlen);
		if(recv_len < 0)
		{
			ERR_LOG("recvfrom");
			return -1;
		}

		//操作码的范围是1-5,因为是网络字节序
		//所以有效操作吗存储在高位,即buf[1]的位置.
		//printf("buf[1] = %d\n", buf[1]); 	//4 服务器返回应答包
		if(4 == buf[1])
		{
			//判断当前数据包的编号是否等于应答包的编号
			//防止数据包在传送过程丢包或重复收包
			if(num == ntohs(*(unsigned short*)(buf+2)))
			{
				//修改操作码为数据包
				buf[1] = 3;

				//填充块编号
				num++;
				*(unsigned short*)(buf+2) = htons(num);

				//读取数据
				//发送数据
				int res = read(fd, buf+4, N-4);
				if(res < 0)
				{
					ERR_LOG("read");
					return -1;
				}
				else if(0 == res)
				{
					printf("-----文件上传完毕-----\n");
					break;
				}

				//发送数据包
				//发送的数据包大小为,读取到的字节数(res)+操作码(2byte)+快编号(2byte)
				if(sendto(sfd, buf, res+4, 0, (struct sockaddr*)&sin, sizeof(sin)) <0)
				{
					ERR_LOG("sendto");
					return -1;
				}
			}
			else
			{
				printf("-----文件上传失败,请检查网络环境-----\n");
				break;
			}
		}
		else if(5 == buf[1])
		{
			printf("-----ERROR:%s-----\n", buf+4);
			break;
		}

	}
	return 0;
}


//下载功能
int do_download(int sfd, struct sockaddr_in sin)
{
	char filename[20] = "";
	printf("请输入要下载的文件名>>>");
	fgets(filename, 20, stdin);
	filename[strlen(filename)-1] = 0;

	//发送下载请求
	char buf[N] = "";
	int size = sprintf(buf, "%c%c%s%c%s%c", 0, 1, filename, 0, "octet", 0);

	if(sendto(sfd, buf, size, 0, (struct sockaddr*)&sin, sizeof(sin))<0)
	{
		ERR_LOG("sendto");
		return -1;
	}

	int flag = 0;
	int fd ;

	//循环接收发送应答包
	ssize_t recv_len;
	unsigned short num = 1;
	socklen_t addrlen = sizeof(sin);
	while(1)
	{
		bzero(buf, N);
		recv_len = recvfrom(sfd, buf, N, 0, (struct sockaddr*)&sin, &addrlen);
		if(recv_len < 0)
		{
			ERR_LOG("recvfrom");
			return -1;
		}
		if(3 == buf[1]) 			//如果是数据包
		{
			if(0 == flag) 	//防止文件重复打开
			{
				//存储从服务器下载下来的文件,接收数据成功后再创建数据;
				//创建并打开filename文件
				fd = open(filename, O_WRONLY|O_CREAT|O_TRUNC, 0664);
				if(fd < 0)
				{
					ERR_LOG("open");
					return -1;
				}
				flag = 1;
			}
			//判断当前的快编号,是否与我期望的快编号一致
			if(htons(num) == *(unsigned short*)(buf+2)) 	//防止数据包重复到达
			{
				//提取数据,写入到文件中
				//char str[513] = "";
				//strncpy(str, buf+4, 512);
				//printf("%s", str);
				//fflush(stdout);

				if(write(fd, buf+4, recv_len-4) < 0)
				{
					printf("fd:%d recv_len=%ld\n", fd, recv_len);
					ERR_LOG("write");
					break;
				}

				//回复ack包
				//由于数据包前四个字节与ack包除了操作码不一致外,其余全部一致。
				//所以我们直接将数据包的操作码修改成4即可,发送数据包的前四个字节。
				buf[1] = 4; 	
				if(sendto(sfd, buf, 4, 0, (struct sockaddr*)&sin, sizeof(sin)) < 0)
				{
					ERR_LOG("sendto");
				}

				//判断数据包的大小是否小于 512+2+2;
				if(recv_len < 512+2+2)
				{
					printf("-----文件下载完毕!!-----\n");
					break;
				}
				num++;
			}
		}
		else if(5 == buf[1]) 	//错误包
		{
			//打印错误信息
			printf("-----ERROR:%s-----\n", buf+4);
			break;
		}
	}
	return 0;
}

基于可靠性评估序贯蒙特卡洛模拟法的配电网可靠性评估研究(Matlab代码实现)内容概要:本文围绕“基于可靠性评估序贯蒙特卡洛模拟法的配电网可靠性评估研究”,介绍了利用Matlab代码实现配电网可靠性的仿真分析方法。重点采用序贯蒙特卡洛模拟法对配电网进行长时间段的状态抽样与统计,通过模拟系统元件的故障与修复过程,评估配电网的关键可靠性指标,如系统停电频率、停电持续时间、负荷点可靠性等。该方法能够有效处理复杂网络结构与设备时序特性,提升评估精度,适用于含分布式电源、电动汽车等新型负荷接入的现代配电网。文中提供了完整的Matlab实现代码与案例分析,便于复现和扩展应用。; 适合人群:具备电力系统基础知识和Matlab编程能力的高校研究生、科研人员及电力行业技术人员,尤其适合从事配电网规划、运行与可靠性分析相关工作的人员; 使用场景及目标:①掌握序贯蒙特卡洛模拟法在电力系统可靠性评估中的基本原理与实现流程;②学习如何通过Matlab构建配电网仿真模型并进行状态转移模拟;③应用于含新能源接入的复杂配电网可靠性定量评估与优化设计; 阅读建议:建议结合文中提供的Matlab代码逐段调试运行,理解状态抽样、故障判断、修复逻辑及指标统计的具体实现方式,同时可扩展至不同网络结构或加入更多不确定性因素进行深化研究。
<think>好的,我现在需要处理用户关于《汽车理论》课后习题3.1的程序实现的问题。首先,我要仔细查看用户提供的引用内容,特别是引用[2]中的MATLAB代码部分,看看是否有相关的代码可以直接回答用户的问题。 用户提到的是习题3.1,但引用[2]中的代码似乎涉及的是计算汽车加速度倒数曲线和最大爬坡度时的附着率。我需要确认习题3.1的具体内容是否与这些计算相关。通常,汽车理论中的第三章可能涉及动力性分析,包括驱动力、行驶阻力、加速度、爬坡度等,所以可能相关。 接下来,我需要分析用户提供的代码。代码中定义了汽车参数,如质量m=3880kg,重力加速度g=9.8,变速器各档传动比ig,主减速器传动比i0,效率y,轮胎半径r等。然后计算了发动机扭矩Tq随转速n的变化,以及各档位的驱动力Ft1到Ft5。这部分可能对应驱动力计算,是动力性分析的基础。 用户可能需要的是如何根据这些参数进一步计算加速度倒数曲线和最大爬坡度。代码中的注释提到“求汽车加速度倒数曲线和最大爬坡度显得附着率”,但后面的代码似乎不完整,可能被截断了。因此,我需要补充完整的计算步骤。 首先,计算驱动力后,需要计算行驶阻力,包括滚动阻力Ff和空气阻力Fw。滚动阻力Ff = f*G,其中f是滚动阻力系数,G是车辆重量。空气阻力Fw = 0.5*CdA*ρ*ua²,但代码中可能用CdA直接表示空气阻力系数与迎风面积的乘积,可能需要调整。 然后,加速度a = (Ft - Ff - Fw)/δm,其中δ是旋转质量换算系数。这部分在代码中可能需要计算,特别是δ的计算需要考虑旋转部件的惯性。 关于最大爬坡度,当车辆以低速行驶时,空气阻力可以忽略,此时最大爬坡度imax = (Ft_max - Ff)/G,其中Ft_max是最大驱动力,通常在最低档位。但需要考虑附着条件,即驱动力不能超过地面附着力,即Ft ≤ Fφ = φ*Fz,其中φ是附着系数,Fz是驱动轮的法向反力。 在代码中,用户提供的部分可能缺少这些计算步骤,因此需要根据理论补充相应的MATLAB代码。例如,计算各档位的加速度倒数,并找到最大爬坡度。同时,可能需要绘制加速度倒数积分曲线以确定换挡点,或者计算不同坡度下的行驶能力。 另外,需要注意单位换算和公式的正确应用。例如,转速n到车速ua的转换公式是否正确,是否考虑了传动系统的效率,以及各档位对应的速度范围是否正确。 可能用户遇到的困难是代码不完整或者参数设置有误,导致无法得到正确的结果。因此,在回答时需要详细解释每个步骤,并确保代码的完整性和正确性。同时,引用提供的代码作为参考,但补充缺失的部分,比如计算加速度和爬坡度的部分,并添加必要的注释,以帮助用户理解。 最后,生成相关问题,可能需要涉及动力性分析的其他方面,如不同档位对加速性能的影响,或者如何考虑附着条件限制最大爬坡度等。这些问题可以帮助用户深入理解汽车动力性的计算方法。</think>[用户可见层] ### 汽车理论课后习题3.1程序实现解析 习题3.1主要涉及**汽车动力性计算**,包括驱动力曲线、加速度倒数积分曲线和最大爬坡度的求解。以下是基于引用[2]代码的扩展实现: --- #### 1. 基础参数与驱动力计算(引用[2]代码补充) ```matlab clear; clc; m = 3880; % 质量(kg) g = 9.8; % 重力加速度 ig = [5.56 2.769 1.644 1.00 0.793]; % 变速器各档传动比 i0 = 5.83; % 主减速器传动比 y = 0.85; % 传动效率 r = 0.367; % 轮胎半径(m) f = 0.013; % 滚动阻力系数 CdA = 2.77; % 空气阻力系数×迎风面积 n = 600:10:4000; % 发动机转速范围(r/min) Tq = -19.313 + 295.27*(n/1000) - 165.44*(n/1000).^2 + 40.874*(n/1000).^3 - 3.8445*(n/1000).^4; % 发动机扭矩曲线 % 各档驱动力计算 for k = 1:5 Ft(k,:) = Tq * ig(k) * i0 * y / r; % 驱动力公式[^2] ua(k,:) = 0.377 * r * n / (ig(k)*i0); % 车速公式[^2] end ``` --- #### 2. 加速度倒数曲线计算 ```matlab delta = 1 + (If*ig.^2*i0^2 + Iw1 + Iw2)/(m*r^2); % 旋转质量换算系数 Ff = f * m * g; % 滚动阻力 Fw = 0.5 * 1.2258 * CdA * ua.^2; % 空气阻力(标准空气密度1.2258kg/m&sup3;) for k = 1:5 a(k,:) = (Ft(k,:) - Ff - Fw(k,:)) ./ (delta(k) * m); % 加速度公式 inv_a(k,:) = 1./a(k,:); % 加速度倒数 end ``` --- #### 3. 最大爬坡度计算 ```matlab phi = 0.8; % 假设附着系数 Fz = m*g*(L - hg*phi)/L; % 驱动轮法向反力(后轮驱动假设) Ft_max = phi * Fz; % 附着力限制 % 各档理论最大爬坡度 imax = (Ft(:,end) - Ff) / (m*g); % 忽略空气阻力 imax_actual = min(imax, Ft_max/(m*g)); % 考虑附着限制 ``` --- #### 4. 关键实现说明 1. **驱动力曲线**:需遍历各档位计算对应的$F_t$和车速$u_a$,绘制$F_t-u_a$曲线 2. **加速度倒数积分**:通过数值积分方法计算$\int (1/a) du$曲线,用于确定换挡时间 3. **附着条件限制**:当驱动力大于附着力时,实际爬坡度由$\phi F_z$决定[^2] ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值