声源定位和GPS模拟实验编程

上传人:h****6 文档编号:210598286 上传时间:2023-05-17 格式:DOCX 页数:4 大小:200.95KB
收藏 版权申诉 举报 下载
声源定位和GPS模拟实验编程_第1页
第1页 / 共4页
声源定位和GPS模拟实验编程_第2页
第2页 / 共4页
声源定位和GPS模拟实验编程_第3页
第3页 / 共4页
资源描述:

《声源定位和GPS模拟实验编程》由会员分享,可在线阅读,更多相关《声源定位和GPS模拟实验编程(4页珍藏版)》请在装配图网上搜索。

1、细心整理1、 声源定位:#include #include using namespace std;#define c 2982#define pai 3.1415926535int main()double x0,x1,x2,y0,y1,y2,t0,t1,t2,t3,A,B,D,angle1,angle2,dt1,dt2,r;x0=0,x1=0,x2=0.3;y0=0,y1=0.45,y2=0.45;int i = 8;while(i-)cout 请输入第8-i组时间数据 t0 t1 t2 t3;t0 = t0 / 1000000;t1 = t1 / 1000000;t2 = t2 / 10

2、00000;dt1 = t0 - t2;dt2 = t1 - t2;A = x2*(pow(x1,2)+pow(y1,2)-pow(c*dt1,2) - x1*(pow(x2,2)+pow(y2,2)-pow(c*dt2,2);B = y2*(pow(x1,2)+pow(y1,2)-pow(c*dt1,2) - y1*(pow(x2,2)+pow(y2,2)-pow(c*dt2,2);D = c*dt1*(pow(x2,2)+pow(y2,2)-pow(c*dt2,2) - c*dt2*(pow(x1,2)+pow(y1,2)-pow(c*dt1,2);angle2 = atan(B/A);i

3、f(acos(D/sqrt(pow(A,2)+pow(B,2)+angle2 pai/2)angle1 = acos(D/sqrt(pow(A,2)+pow(B,2)+angle2;elseangle1 = angle2 - acos(D/sqrt(pow(A,2)+pow(B,2);r = (pow(x1,2)+pow(y1,2)-pow(c*dt1,2)/(2*(x1*cos(angle1)+y1*sin(angle1)+c*dt1);cout 声源坐标为:r*cos(angle1)*1000mm , r*sin(angle1)*1000mmendl; / end while;return

4、 0; / end main;2、 GPS模拟#include #include #include using namespace std;#define c 2982#define min 0.0000000001void GPS(double x10 , double y10 , double t10 , double x1 , double y1 , int num)double a1,a2,a3,a4,b1,b2,mult,D=0,E=0;int i;while(1)a1 = a2 = a3 = a4 = b1 = b2 = 0;for(i = 0; i 10; i +)a1 = a1

5、-(3*pow(xi-x1),2)+pow(yi/1000-y1),2)-pow(c*ti/1000000,2);a2 = a2+2*(xi/1000-x1)*(y1-yi/1000);a3 = a2;a4 = a4-(3*pow(yi-y1),2)+pow(xi/1000-x1),2)-pow(c*ti/1000000,2);b1 = b1+(pow(xi/1000-x1),2)+pow(yi/1000-y1),2)-pow(c*ti/1000000,2)*(x1-xi/1000);b2 = b2+(pow(xi/1000-x1),2)+pow(yi/1000-y1),2)-pow(c*ti

6、/1000000,2)*(y1-yi/1000); /end for;mult = a3/a1;a4 -= a2*mult;b2 -= b1*mult;E = b2/a4;D = (b1-a2*E)/a1;x1 += D;y1 += E;if(pow(D,2)+pow(E,2) min)break; / end while;x1 = x1*1000;y1 = y1*1000;cout 点P num 定位于: X = fixed setprecision(8) x1;cout (mm) , Y = fixed setprecision(8) y1 (mm)endl; / end GPS;int main()int yy;double t110 , t210 , t310 , x10 , y10;cout 请输入时间数据(us)与坐标数据(mm): endl;int i;for(i = 0; i = 9; i +)cout 第 i+1 组: t1i t2i t3i xi yi; / end for;GPS(x,y,t1,0,0.45,1);GPS(x,y,t2,0.3,0.45,2);GPS(x,y,t3,0,0,3);cinyy;return 0; / end main;

展开阅读全文
温馨提示:
1: 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
2: 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
3.本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
5. 装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2023-2025  zhuangpeitu.com 装配图网版权所有   联系电话:18123376007

备案号:ICP2024067431-1 川公网安备51140202000466号


本站为文档C2C交易模式,即用户上传的文档直接被用户下载,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有。装配图网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对上载内容本身不做任何修改或编辑。若文档所含内容侵犯了您的版权或隐私,请立即通知装配图网,我们立即给予删除!