Gimme梦
Gimme梦
关注数: 16
粉丝数: 78
发帖数: 453
关注贴吧数: 12
业余爱好——Tropical house[br]哈哈 图一乐 就贼觉得有成就感
安徽往届高考生成绩也能查询了 安徽老考生突然想重温高考查分的那一天了 虽然分不高 但在我心里已经达到预期了 当年这分数 只能上边远地区的211 但对于农村出身的我 确是相当不容易的了 怀念那时的激动心情 也相信自己一直无畏地走下去
FL-Studio-Arp-Bass [br]业余爱好,弄了一个主要以Arp和Bass为特色的东西
溪涧草堂pt站是不是没了? 最近打不开了,不会是无限关停了吧?
◇07-30◆〖互帮互助〗安徽除了这个吧,还有什么吧吗? 看吧里的帖子,了解到中西部主要是河南人的天下,中华城市主要是川渝人,安徽是不是没有地区吧呢
Ios15.6.1怎样 手机iPhone13,去年9月官网抢的,还是出厂版本ios15.0 目前电池健康还剩90% 感觉想通过升级提升一点续航 有没有铁汁介绍下ios15.6.1续航跟15.0比怎样
今天的更新 2022-KB5013628 刚刚更新 May 10, 2022-KB5013628 Cumulative Update for .NET Framework 3.5 and 4.8 for Windows 11 This security update addresses an issue where a local user opening a specially crafted file could cause a denial of service condition on an affected system. For more information please see CVE-2022-30130. 最重要的更新:好像是终于对文件操作假死问题的更新,不知道有没有用。先试着用一两天看看。
Windows11到底为什么打不开文件夹 3月末更新了win11,好家伙,最大的问题就是文件夹有很大几率打不开,还有有时打开文件也卡死,最无语的是昨天终于更新了新的版本,上面说的问题非但没解决,连桌面都有几率卡死了。谁来帮帮我,退回win10最大的问题是很多专业软件非常难安装
吧里的🐭🐭们怎么不整活了 逛吧,怎么没有小故事让👨🏻🦳开心呢
西安有发货的吗? 等一个,到现在还是准备发货,看看有没有人发了
果 粉 过 年 现 状 由于换机来这里看看,在这里,我见识到了很多。1.很多人吐槽抢机的困难和漫长的交货日期,并都热心地展示了自己的订单截图,其中以13p蓝尤甚,吓得我这土鳖下了13的单瑟瑟发抖不敢吱声;2.偶有寻求购机意见贴,必然推荐13pm,宣称性价比极高,可是我看了跑分图,貌似几个差别不大,一定是我这个土狗仅看评分下结论了,因为很多果粉告诉我,跑分不能反映真实的性能。从这两点可以看出,很多果粉都是追求性价比的高端人群。不说话的都是和我一样的土鳖
有没有吧里的大神投过Physica A啊? 如图,这是我2月7日投的,2月8号with editor,然后系统状态一直都停在2月7号了,今天都已经13号了,这是什么情况?这家杂志的效率还是比较快的,是拒是送审很快就有结果的,前些日子投过另一篇文章,从with editor后第三天主编就直接拒稿了。。。
大家觉得我的小8需要换电池吗 兄嘚们,我是18年9月买的手机,到现在什么事都没有,唯一的缺点就是感觉电池不耐用了,现在至少一天三冲,主要是我有个习惯,低于50的电量就情不自禁想冲,最近感觉掉电快,于是想换个大容量(2180mA)的电池,大家觉得我有必要换吗
非线性常微分方程组的解失真是怎么回事? 如上图,要计算一个非线性常微分方程组,采用的是四阶龙格库塔算法,得到的结果在0-0.6左右计算还是可以的,但到后来的计算结果就完全失真发散了。网上查找了很多资料也没找到什么原因,有人能帮忙看一下吗?附带代码如下: #include <stdio.h> #include <math.h> #define M 3 // 维数,方程组中方程的个数// 算斜率的右端函数 void RightSlopeFun( double *pRightK, //算出来的右端斜率项[out] double iT, //计算点的微分自变量值 double *iY //计算点的函数初始值 ) { pRightK[0] = -(iY[0]-iY[1])-(2*iY[0]*iY[2]+iY[1]*iY[0]+iY[1]*iY[2]-iY[2]*iY[2])+3; pRightK[1] = -(iY[1]+iY[2])-(2*iY[0]*iY[0]/3-iY[0]*iY[1]/3+iY[1]*iY[0]-2*iY[1]*iY[2]/3+2*iY[2]*iY[0]/3+4*iY[2]*iY[1]/3+iY[2]*iY[2]/3)+1; pRightK[2] = -iY[1]-(-iY[0]*iY[0]/3+2*iY[0]*iY[1]/3-iY[0]*iY[2]-iY[1]*iY[0]+iY[1]*iY[2]/3-iY[2]*iY[0]/3-2*iY[2]*iY[1]/3+4*iY[2]*iY[2]/3)-1; } // 龙格库塔法解算程序 double LgktSolution( double *pOY, // 返回函数值[out] double iT, // 积分自变量输入值 double *iY, // 初始值 double h, // 单步步长 int Count=1 // 解算步数,默认形参为1,调用时不给此 // 参数赋值,则只解算一步 ) { double K1[M]; // 第一个斜率 double K2[M]; // 第二个斜率 double K3[M]; // 第三个斜率 double K4[M]; // 第四个斜率 double tempY[M]; // 保存Y的中间值的临时空间 int i; // 计数器 double tempT; // 临时保存x的空间 tempT=iT; // 载入初始的x值 for (i=0;i<M;i++) { tempY[i]=iY[i]; // 载入初始的函数值 pOY[i]=iY[i]; } while (Count>0) { RightSlopeFun(K1,tempT,tempY); // 解算斜率K1 tempT += h/2; for (i=0;i<M;i++) { pOY[i]=tempY[i]+K1[i]*h/2; } RightSlopeFun(K2,tempT,pOY); // 解算斜率K2 for (i=0;i<M;i++) { pOY[i]=tempY[i]+K2[i]*h/2; } RightSlopeFun(K3,tempT,pOY); // 解算斜率K3 tempT += h/2; for (i=0;i<M;i++) { pOY[i]=tempY[i]+K3[i]*h; } RightSlopeFun(K4,tempT,pOY); // 解算斜率K4 for (i=0;i<M;i++) { // 得到函数值 pOY[i]=tempY[i]+h*(K1[i]+2*K2[i]+2*K3[i]+K4[i])/6; // 推算了一步的函数值,此处可以对该数据进行相应的操作 } // 递推数据更新 Count--; if (Count<=0) { break; } for (i=0;i<M;i++) { tempY[i]=pOY[i]; } } return tempT; } void main() { FILE *fp=fopen("g:\\随机振动相关论文及程序代码\\试算\\非平稳解\\测试组\\无耦合项.txt","a"); double iy[M]; double it; int i; iy[0]=1; iy[1]=0.00001; iy[2]=0.00001; it=0; for(i=0;i<1000;i++) { it=LgktSolution(iy,it,iy,0.01); fprintf(fp,"%-8.2f %-8.6f %-8.6f %-8.6f %-8.6f %-8.6f %-8.6f\n",it,iy[0],iy[1],iy[2],iy[3],iy[4],iy[5],iy[6],iy[7],iy[8],iy[9],iy[10],iy[11],iy[12],iy[13],iy[14],iy[15]); } fclose(fp); }
1
下一页