用MATLAB计算和绘制向量场、幅度场、拉格朗日张量场? 5C

Given:
▪(x, y) coordinates of the reference and deformed points
▪Ref–reference points
▪Def–deformed points

    Task:

▪Calculate and plot the displacement field
include:Vector field and Magnitude field
▪Calculate and plot the Lagrangian strain field

其中:
ref为
图片说明

def为
图片说明

需要MATLAB的编程代码,谢谢!

1个回答

qq_42026116
Lws63000 刚写完第一题、
一年多之前 回复
qq_42026116
Lws63000 回复weixin_42555055: 我觉得吧 还是自己写比较靠谱 hhh 竟然能遇上你。
一年多之前 回复
weixin_42555055
澄心伊兰 能否给出MATLAB的编程代码?谢谢!
一年多之前 回复
Csdn user default icon
上传中...
上传图片
插入图片
抄袭、复制答案,以达到刷声望分或其他目的的行为,在CSDN问答是严格禁止的,一经发现立刻封号。是时候展现真正的技术了!
其他相关推荐
用matlab进行拉格朗日插值还原曲线,和分段插值,详细代码包括测试结果

还原图1.4(绘制Lagrange 10次插值曲线和原函数曲线),并加上图1.5中的分段线性差值曲线。(三条曲线绘制在同一张Figure里面,注意用不同的线形区分) ![图片说明](https://img-ask.csdn.net/upload/201803/28/1522225085_470202.jpg) ![图片说明](https://img-ask.csdn.net/upload/201803/28/1522225106_459261.jpg)

【求帮助】拉格朗日插值法-补值

setwd("C:/Users/DELL1/Desktop") library(XLConnect) missing_data<-XLConnect::readWorksheetFromFile(file = "./data/missing_data.xls",sheet = 1,header=F) lagrange<-function(x,xi,yi){ n<-length(xi) lage<-0 for(i in 1:n){ li<-1 for(j in 1:n){ if(i!=j) li<-li*(x-xi[j])/(xi[i]-xi[j]) } lage<-li*yi[i]+lage } return(lage) } missdata = missing_data for(k in 1:3) { x = which(is.na(missing_data[,k])) x1=c(0,x) x2=c(x,nrow(missing_data)) x12 = x2-x1-1 xx1<-x12[1:(length(x12)-1)] #缺失值前面的行数 xx2<-x12[2:(length(x12))] #缺失值后面的行数 j = 1 for(m in x) { if(xx1[j]>=5) #空值前的判断 xi<- (m-5):(m-1) else xi<-(m-xx1[j]):(m-1) if(xx2[j]>=5) #空值后的判断 xi<-c(xi,(m+1):(m+5)) else xi<-c(xi,(m+1):(m+xx2[j])) yi=missing_data[xi,k] missdata[m,k]=lagrange(m,xi,yi) print(c(m,missdata[m,k])) j =j+1 } } 数据见下(无法上传附件,抱歉) 21*3的表格 235.8333 324.0343 478.3231 236.2708 325.6379 515.4564 238.0521 328.0897 517.0909 235.9063   514.89 236.7604 268.8324     404.048 486.0912 237.4167 391.2652 516.233 238.6563 380.8241   237.6042 388.023 435.3508 238.0313 206.4349 487.675 235.0729     235.5313 400.0787 660.2347   411.2069 621.2346 234.4688 395.2343 611.3408 235.5 344.8221 643.0863 235.6354 385.6432 642.3482 234.5521 401.6234   236 409.6489 602.9347 235.2396 416.8795 589.3457 235.4896   556.3452 236.9688   538.347 补充的值不对!!! 麻烦指出一下,谢谢 [1] 6.0000 237.1512 [1] 13.000 235.315 [1] 4.0000 203.4621 [1] 11.0000 237.3481 [1] 20 NA [1] 21 NaN [1] 5.0000 503.7444 [1] 8.0000 472.0948 [1] 11.0000 634.3942 [1] 17.0000 615.7713

如何用C++绘制插值得到的曲线并取值?

如何用C++用插值的方法得到曲线并取值? 有一些数据点(准确的,所以用插值,点不多,30个左右) 插值出一个曲线 取点(x是整数,y两位小数)(对任意x得到y)

基于双线性插值的拉格朗日像素修复算法

工程编译环境vc6.0+opencv1.0 ``` #include <stdlib.h> #include <stdio.h> #include <iostream.h> #include <math.h> #include <cv.h> #include <highgui.h> #ifdef _CH_ #pragma package <opencv> #endif #include <stdio.h> #include <stdlib.h> struct record//两个代表偏移 { int up,down,left,right;//差距 unsigned char dataup; unsigned char datadown; unsigned char dataleft; unsigned char dataright;// int **BT; };//负责记录空间位置,和边界像素值//换成指针结构体? IplImage* inpaint_mask = 0; IplImage* img0 = 0; IplImage* img = 0; IplImage* inpainted = 0; CvPoint prev_pt = {-1,-1}; unsigned char aver(unsigned char *data,unsigned char *bill,IplImage* masktt,int i,int j,int k,int step,int channels); void scan(unsigned char *data,unsigned char *bill,IplImage* masktt,record *re,int i,int j,int k,int step,int channels);//分别扫描出左右上下最短距离赋给t采样 unsigned char biline(record *re,int wh,int i,int j); double lagrange(double *x,double *y,double xx,int n); double xlag(record re[],int i); double ylag(record re[],int j); void cvrepair(IplImage* img,IplImage* mask,IplImage* masktt);//img破坏图像,mask轨迹母板,masktt2值图像 void on_mouse( int event, int x, int y, int flags, void* zhang) { if( !img ) return; if( event == CV_EVENT_LBUTTONUP || !(flags & CV_EVENT_FLAG_LBUTTON) ) prev_pt = cvPoint(-1,-1); else if( event == CV_EVENT_LBUTTONDOWN ) prev_pt = cvPoint(x,y); else if( event == CV_EVENT_MOUSEMOVE && (flags & CV_EVENT_FLAG_LBUTTON) ) { CvPoint pt = cvPoint(x,y); if( prev_pt.x < 0 ) prev_pt = pt; cvLine( inpaint_mask, prev_pt, pt, cvScalarAll(255), 15, 8); cvLine( img, prev_pt, pt, CV_RGB(0, 255, 0), 15, 8); prev_pt = pt; cvShowImage( "image", img ); } } int main() { char strIn[1024]; char strOut[1024]; cout<<"请输入原始图片路径:"; cin>>strIn; //cout<<endl; cout<<"请输入结果图片路径: "; cin>>strOut; if( (img0 = cvLoadImage(strIn,-1)) == 0 ) return 0; printf( "Hot keys: \n" "\tESC - quit the program\n" "\tr - restore the original image\n" "\ti or ENTER - run inpainting algorithm\n" "\t\t(before running it, paint something on the image)\n" ); cvNamedWindow( "image", 1 ); img = cvCloneImage( img0 ); inpainted = cvCloneImage( img0 ); inpaint_mask = cvCreateImage( cvGetSize(img), 8, 1 ); cvZero( inpaint_mask ); cvZero( inpainted ); cvShowImage( "image", img ); //cvShowImage( "watershed transform", inpainted ); cvSetMouseCallback( "image", on_mouse, 0 ); for(;;) { int c = cvWaitKey(0); if( (char)c == 27 ) break; if( (char)c == 'r' ) { cvZero( inpaint_mask ); cvCopy( img0, img,0 ); cvShowImage( "image", img ); } if( (char)c == 'i' || (char)c == '\n' ) { cvNamedWindow( "Mask image", 1 ); cvShowImage( "Mask image", inpaint_mask ); cvWaitKey(0); cvNamedWindow( "inpainted image", 1 ); cvrepair( img, inpaint_mask,inpainted); //cvInpaint( img, inpaint_mask, inpainted, 3, CV_INPAINT_TELEA ); cvShowImage( "inpainted image", inpainted ); cvSaveImage(strOut, inpainted); } } return 1; } void cvrepair(IplImage* img,IplImage* mask,IplImage* inpainted)//img破坏图像,mask轨迹母板,masktt2值图像 { int x=img->height; int y=img->width; int step=img->widthStep; int channels=img->nChannels; unsigned char *origin=(uchar *)inpainted->imageData; unsigned char *data=(uchar *)img->imageData; int i,j,k; unsigned char *bill=(uchar *)mask->imageData; for(i=0;i<x;i++) for(j=0;j<y;j++) { for( k=0;k<channels;k++) origin[i*step+j*channels+k]=data[i*step+j*channels+k];//原图 } for(i=0;i<x;i++)//边界 { for(j=0;j<y;j++) { if(bill[i*step+j*channels]==1) { for(k=0;k<channels;k++) { unsigned char buff=aver(origin,bill,mask,i,j,k,step,channels); origin[i*step+j*channels+k]=buff;//像素点覆 } } } } } unsigned char aver(unsigned char data[],unsigned char bill[],IplImage* masktt,int i,int j,int k,int step,int channels) { record re[2]; scan(data,bill,masktt,re,i,j,k,step,channels);//扫描采样->re int len=re[1].down-re[1].up; int lon=re[1].right-re[1].left; int wh=len/(lon+len);//高宽比 return biline(re,wh,i,j);//目标还原像素值 } void scan(uchar data[],uchar bill[],IplImage* masktt,record *re,int i,int j,int k,int step,int channels)//分别扫描出左右上下最短距离赋给t采样 { int m,p,n,q; int h=masktt->height; int w=masktt->width; for(m=0;m<i;m++) { if(bill[(i-m)*step+j*channels]==0) { if(m!=i)//上边界有像素 完成2个采样,边缘记0 { re[1].dataup=data[(i-m)*step+j*channels+k]; re[1].up=m;//上边界第1个最近 if((bill[(i-m-1)*step+j*channels]==0)&&((m+1)!=i)) { re[2].dataup=data[(i-m-1)*step+j*channels+k];re[2].up=m+1;}//上边界第2个最近 else{ re[2].dataup=0;re[2].up=i;} } else//上边界有像素破坏 { re[1].up=i;re[1].dataup=0; re[2].up=i;re[2].dataup=0; } printf("find up");break; } } for(p=i;p<h;p++) { if(bill[p*step+j*channels]==0) { if(p!=h)//下边界有像素 { re[1].datadown=data[p*step+j*channels+k];//下边界第1个最近 re[1].down=p; if((bill[(p+1)*step+j*channels]==0)&&((p+1)!=h)) { re[2].datadown=data[(p+1)*step+j*channels+k];re[2].down=p+1; }//下边界第2个最近 else{ re[2].datadown=0;re[2].down=h;} } else { re[1].down=h;re[1].datadown=0; re[2].down=h;re[2].datadown=0; } printf("find down");break; } } for(n=0;n<j;n++) { if(bill[i*step+(j-n)*channels]==0) { if(n!=j)//上边界有像素 完成2个采样,边缘记0 { re[1].dataleft=data[i*step+(j-n)*channels+k]; re[1].left=n;//上边界第1个最近 if((bill[i*step+(j-n-1)*channels]==0)&&((n+1)!=j)) { re[2].dataleft=data[i*step+(j-n-1)*channels+k];re[2].left=n+1;}//上边界第2个最近 else { re[2].dataleft=0;re[2].left=j;} } else//上边界有像素破坏 { re[1].left=j;re[1].dataleft=0; re[2].left=j;re[2].dataleft=0; } break; } } for(q=j;q<w;q++) { if(bill[i*step+q*channels]==0) { if(q!=w)//下边界有像素 { re[1].dataright=data[i*step+q*channels+k];//下边界第1个最近 re[1].right=q; if((bill[i*step+(q+1)*channels]==0)&&((q+1)!=w)) { re[2].dataright=data[i*step+(q+1)*channels+k];re[2].right=q+1; }//下边界第2个最近 else{ re[2].dataright=0;re[2].right=w;} } else { re[1].right=w;re[1].dataright=0; re[2].right=w;re[2].dataright=0; } break; } } } unsigned char biline(record *re,int wh,int i,int j)//双线性插值算法 { double temp=(double)wh*ylag(re,j)+(double)(1-wh)*xlag(re,i);//强制转型 return (uchar)temp; } double ylag(record re[],int j) { double x[3],y[3]; x[0]=re[2].left;y[0]=re[2].dataleft; x[1]=re[1].left;y[1]=re[1].dataleft; x[2]=re[1].left;y[2]=re[1].dataleft; x[3]=re[2].left;y[3]=re[2].dataleft; return lagrange(x,y,j,4); } double xlag(record re[],int i) { double x[3],y[3]; x[0]=re[2].up;y[0]=re[2].dataup; x[1]=re[1].up;y[1]=re[1].dataup; x[2]=re[1].down;y[2]=re[1].datadown; x[3]=re[2].down;y[3]=re[2].datadown; return lagrange(x,y,i,4); } double lagrange(double *x,double *y,double xx,int n) /*拉格朗日插值算法*/ { int i,j; double *a,yy=0.0; /*a作为临时变量,记录拉格朗日插值多项式*/ a=(double *)malloc(n*sizeof(double)); for(i=0;i<=n-1;i++) { a[i]=y[i]; for(j=0;j<=n-1;j++) if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i]; } free(a); return yy; } ``` 然后在调试的时候报错 unhandled exception in tss3.exe:0xC0000005:Access Violation. 停下的位置在scan()函数里的 ``` for(p=i;p<h;p++) { if(bill[p*step+j*channels]==0) { if(p!=h)//下边界有像素 { re[1].datadown=data[p*step+j*channels+k];//下边界第1个最近 re[1].down=p; if((bill[(p+1)*step+j*channels]==0)&&((p+1)!=h)) { re[2].datadown=data[(p+1)*step+j*channels+k];re[2].down=p+1; }//下边界第2个最近 else{ re[2].datadown=0;re[2].down=h;} } else { re[1].down=h;re[1].datadown=0; re[2].down=h;re[2].datadown=0; } printf("find down");break; } } ``` 的第二次循环 时间比较紧迫所以就不等了

matlab中怎么检测一幅二值图中的线是直线还是曲线?

一幅二值图中有几条线,怎么判断出它是直线还是曲线?![图片说明](http://forum.csdn.net/PointForum/ui/scripts/csdn/Plugin/001/face/54.gif)![图片说明](http://forum.csdn.net/PointForum/ui/scripts/csdn/Plugin/001/face/54.gif)![图片说明](http://forum.csdn.net/PointForum/ui/scripts/csdn/Plugin/001/face/54.gif)![图片说明](http://forum.csdn.net/PointForum/ui/scripts/csdn/Plugin/001/face/54.gif)![图片说明](http://forum.csdn.net/PointForum/ui/scripts/csdn/Plugin/001/face/54.gif)

拉格朗日插值查找

拉格朗日插值查找法:查找整型数据是非常快的 我有一个问题: 拉格朗日插值查找法:适合查找字符串吗? 求各位大神解释一下

帮帮忙吧,急哭了:python实现拉格朗日插值进行缺失值插补时报错,出现Keyerror:0问题?

![图片说明](https://img-ask.csdn.net/upload/201911/13/1573618586_795847.png) 数据截图 ``` import pandas as pd from scipy.interpolate import lagrange inputfile = './test1.xlsx' outputfile = './test1_1.xlsx' def ployinterp_column(s,n,k=6): y = s[list(range(n-k,n)) + list(range(n+1,n+1-k))] y = y[y.notnull()] return lagrange(y.index,list(y))(n) for i in data.columns: for j in range(len(data)): if (data[i].isnull())[j]: data[i][j] = ployinterp_column(data[i],j) ``` 报错 ``` KeyError Traceback (most recent call last) <ipython-input-28-02cf873cddff> in <module>() 1 for i in data.columns: 2 for j in range(len(data)): ----> 3 if (data[i].isnull())[j]: 4 if (data[i].isnull())[j] in data: 5 data[i][j] = ployinterp_column(data[i],j) C:\Anaconda3\lib\site-packages\pandas\core\series.py in __getitem__(self, key) 765 key = com._apply_if_callable(key, self) 766 try: --> 767 result = self.index.get_value(self, key) 768 769 if not is_scalar(result): C:\Anaconda3\lib\site-packages\pandas\core\indexes\base.py in get_value(self, series, key) 3116 try: 3117 return self._engine.get_value(s, k, -> 3118 tz=getattr(series.dtype, 'tz', None)) 3119 except KeyError as e1: 3120 if len(self) > 0 and self.inferred_type in ['integer', 'boolean']: pandas\_libs\index.pyx in pandas._libs.index.IndexEngine.get_value() pandas\_libs\index.pyx in pandas._libs.index.IndexEngine.get_value() pandas\_libs\index.pyx in pandas._libs.index.IndexEngine.get_loc() pandas\_libs\hashtable_class_helper.pxi in pandas._libs.hashtable.Int64HashTable.get_item() pandas\_libs\hashtable_class_helper.pxi in pandas._libs.hashtable.Int64HashTable.get_item() KeyError: 0 ```

拉格朗日插值查找的疑问

最近学到拉格朗日插值查找,发现了一个bug 网上找到的写法大多是这样 ```c void lagrangeSearch(int num, int *a, int length) { int count = 0; //记录查找几次,对比二分查找 int flag = 0; //标志是否找到 int max = length - 1; int min = 0; //int mid = (max+min)/2 //容易溢出 int mid; while (min <= max) { count++; //没有判断num-a[min]是否小于0,如果小于0不就出错了吗,mid有可能变成负数,造成不正确的访问 mid = min + (max - min) * 1.0 *(num - a[min]) / (a[max] - a[min]); if (num > a[mid]) { min = mid + 1; } else if (num < a[mid]) { max = mid - 1; } else { flag = 1; break; } } if (flag == 1) { printf("拉格朗日找到了,下标为:%d,查找了%d次\n", mid, count); } else { printf("不存在\n"); } } ``` 如代码所示(mid = min + (max - min) * 1.0 *(num - a[min]) / (a[max] - a[min]);) **如果不加以判断的话,num-a[min]是有可能小于0的,造成整个式子小于0,mid就变成了负数,如此访问的话不就访问到了错误的地址吗?** 请问我考虑的对吗? 自己改动后的: ```c //拉格朗日查找 void lagrangeSearch(int num, int *a, int length) { int count = 0; //记录查找几次,对比二分查找 int flag = 0; //标志是否找到 int max = length - 1; int min = 0; //int mid = (max+min)/2 //容易溢出 int mid; //当num<a[min]时,拉格朗日差值为负,无法正常判断 if (num >= a[min] && num <= a[max]) { while (min <= max) { count++; mid = min + (max - min) * 1.0 *(num - a[min]) / (a[max] - a[min]); if (num > a[mid]) { min = mid + 1; } else if (num < a[mid]) { max = mid - 1; } else { flag = 1; break; } } } if (flag == 1) { printf("拉格朗日找到了,下标为:%d,查找了%d次\n", mid,count); } else { printf("不存在\n"); } } ```

python3.8运行拉格朗日插值报错Passing list-likes to .loc or [] with any missing labels is no longer supported

在anaconda里python3.8运行拉格朗日插值代码时提示错误 请大神看看应该怎么改呢? 小白谢过! excel下载链接 链接:https://pan.baidu.com/s/1IYAPr-4UFnHZNZk_vLHhiw 提取码:ejwn #拉格朗日插值代码 import pandas as pd #导入数据分析库Pandas from scipy.interpolate import lagrange #导入拉格朗日插值函数 inputfile = '../data/catering_sale.xls' #销量数据路径 outputfile = '../tmp/sales.xls' #输出数据路径 data = pd.read_excel(inputfile) #读入数据 data[u'销量'][(data[u'销量'] < 400) | (data[u'销量'] > 5000)] = None #过滤异常值,将其变为空值 #自定义列向量插值函数 #s为列向量,n为被插值的位置,k为取前后的数据个数,默认为5 def ployinterp_column(s, n, k=5): y = s[list(range(n-k, n)) + list(range(n+1, n+1+k))] #取数 y = y[y.notnull()] #剔除空值 return lagrange(y.index, list(y))(n) #插值并返回插值结果 #逐个元素判断是否需要插值 for i in data.columns: for j in range(len(data)): if (data[i].isnull())[j]: #如果为空即插值。 data[i][j] = ployinterp_column(data[i], j) data.to_excel(outputfile) #输出结果,写入文件 运行错误提示: KeyError: 'Passing list-likes to .loc or [] with any missing labels is no longer supported, see https://pandas.pydata.org/pandas-docs/stable/user_guide/indexing.html#deprecate-loc-reindex-listlike'

哪位大神给写一个oracle中关于拉格朗日插值法的函数??大神 求指教啊

哪位大神给写一个oracle中关于拉格朗日插值法的函数??大神 求指教啊

能否将空间中曲线投影在三个坐标轴上来唯一确定曲线?

小弟在做一个空间中插值算法,参考资料都是一些关于曲面的,我想做曲线插值找不到具体算法。就像三个点拟合光滑曲线这种。和平面上一样!在平面上的算法都可以明白,像样条插值,拉格朗日插值等等…但是不知道也不会如何直接运用在空间中,因为空间曲线是由两个平面确定的,解析式都不一样,根本不知道如何写公式!请各位大佬帮帮忙![图片说明](http://forum.csdn.net/PointForum/ui/scripts/csdn/Plugin/001/face/1.gif)现在有一个想法,能否将空间中的点投影在三个坐标轴面上,使用平面上的插值拟合算法求出后还原回空间,不知道这样是否可取,大佬们帮忙解惑。万分感谢!

基于模糊身份的加密算法中拉格朗日系数要怎么理解,为什么是连乘,i,j要怎么理解,急求大神解答

![图片说明](https://img-ask.csdn.net/upload/201610/18/1476800370_352384.png)

如何练习java web应用的优化?

对于自己写好的java小型购物网站,数据库数据不多,应该从哪方面的性能优化入手,需要手动插入很多数据,以此来体现优化效果吗?

算法加密的优化 写出你的优化率

加密算法的优化: 任意加密算法,如ECC,RSA,IDEA,3DES,AES,SMS4等 请写出你的优化点及优化率

将Java改成Python中插值问题怎么解决?

Java代码为 PolynomialSplineFunction estimateFunc = lit.interpolate(sx, sy); double[] x = new double[num_cy]; for (int i = 0; i < num_cy; i++) { x[i] = Math.round(buchang) / (double)(num_cy+1) * i; if(x[i]>=Math.round(buchang)-1) x[i]=Math.round(buchang)-2; GLSJ144[i] = estimateFunc.value(x[i]); if (estimateFunc.isValidPoint(x[i])) { GLSJ144[i] = estimateFunc.value(x[i]); } else { GLSJ144[i] = 0.0; } } python改成这样不对为什么? estimateFunc = interp1d(sx, sy,kind='linear') x=list(range(self.num_cy)) for i in range(0, self.num_cy): x[i] = round(buchang)/(self.num_cy+1) * i if x[i]>=round(buchang)-1: x[i]=round(buchang)-2 GLSJ144[i] = estimateFunc(x[i]) print((x[i])) if estimateFunc(x[i]): GLSJ144[i] = estimateFunc(x[i]) else: GLSJ144[i] = 0.0 测试发现两个代码的estimateFunc(x[i])值不一样。求大佬帮忙!

利用群论的相关知识证明费马小定理

题目要求用代数结构里群的相关知识证明费马小定理,急!求帮!谢谢各位大神!

KKT条件 是否是对偶问题等价的条件

![图片说明](https://img-ask.csdn.net/upload/201801/19/1516371714_949476.jpg) 这是我在某PPT里面看到的 然而我看了一些文章和西瓜书后总结出KKT条件是 对于 x∈D, max L(x, lambda, mu) = f(x)的条件 ![图片说明](https://img-ask.csdn.net/upload/201801/19/1516372126_180896.jpg) 然后这两个问题的解在值上相等的条件是西瓜书里说的f和g为凸函数,h为仿射函数(Page 406)不知道值相等是否是PPT中所指的等价? 请教我的理解和第一张PPT截图里的内容哪一个是对的。

opencv窗口是透明的,调不了,留残影

``` #include <iostream> #include <process.h> #include "stdafx.h" #include <GL/glew.h> #include <GL/glut.h> #pragma comment(lib, "glew32.lib") using namespace std; void lagrange() { float p[10][2]; p[0][0] = 5, p[0][1] = 31.2; p[1][0] = 7, p[1][1] = 43.7; p[2][0] = 13, p[2][1] = 59.1; p[3][0] = 16, p[3][1] = 72; p[4][0] = 20, p[4][1] = 94.5; p[5][0] = 25, p[5][1] = 122; p[6][0] = 27, p[6][1] = 129.2; p[7][0] = 34, p[7][1] = 153; p[8][0] = 52, p[8][1] = 183; p[9][0] = 60, p[9][1] = 197.8; //得到数据点 int nmax = 10; int n = 4; float x, y; int i, j, t; float temp; glColor3f(0.0, 0.0, 1.0); for (i = 0; i <= (nmax - 1); i += (n - 1)) { glBegin(GL_LINE_STRIP); for (x = p[i][0]; x <= p[i + n - 1][0]; x += 1.0) { y = 0.0; for (j = 0; j<n; j++) { temp = 1.0; for (t = 0; t<n; t++) { if (t == j) continue; temp *= (x - p[i + t][0]) / (p[i + j][0] - p[i + t][0]); } y += temp*p[i + j][1]; } glVertex2f(x/nmax, y/200); } glEnd(); } glFlush(); } //拉格朗日插值法&画图 int main(int argc, char *argv[]) { glutInit(&argc, argv); glutInitDisplayMode(GLUT_RGB | GLUT_SINGLE); glutInitWindowSize(600, 600); glutCreateWindow("The first OpenGL Application"); glClearColor(0.0, 0.0, 0.0, 0.0); glutDisplayFunc(&lagrange); glutMainLoop(); return 0; //lagrange(p, 4, 10); } ``` ![图片说明](https://img-ask.csdn.net/upload/201905/07/1557194219_810530.png) 如图,透明的,我咋调glclearcolor的透明度都这样

大学四年自学走来,这些私藏的实用工具/学习网站我贡献出来了

大学四年,看课本是不可能一直看课本的了,对于学习,特别是自学,善于搜索网上的一些资源来辅助,还是非常有必要的,下面我就把这几年私藏的各种资源,网站贡献出来给你们。主要有:电子书搜索、实用工具、在线视频学习网站、非视频学习网站、软件下载、面试/求职必备网站。 注意:文中提到的所有资源,文末我都给你整理好了,你们只管拿去,如果觉得不错,转发、分享就是最大的支持了。 一、电子书搜索 对于大部分程序员...

在中国程序员是青春饭吗?

今年,我也32了 ,为了不给大家误导,咨询了猎头、圈内好友,以及年过35岁的几位老程序员……舍了老脸去揭人家伤疤……希望能给大家以帮助,记得帮我点赞哦。 目录: 你以为的人生 一次又一次的伤害 猎头界的真相 如何应对互联网行业的「中年危机」 一、你以为的人生 刚入行时,拿着傲人的工资,想着好好干,以为我们的人生是这样的: 等真到了那一天,你会发现,你的人生很可能是这样的: ...

程序员请照顾好自己,周末病魔差点一套带走我。

程序员在一个周末的时间,得了重病,差点当场去世,还好及时挽救回来了。

ArrayList源码分析(入门篇)

ArrayList源码分析 前言: 写这篇博客的主要原因是,在我上一次参加千牵科技Java实习生面试时,有被面试官问到ArrayList为什么查找的速度较快,插入和删除的速度较慢?当时我回答得不好,很大的一部分原因是因为我没有阅读过ArrayList源码,虽然最后收到Offer了,但我拒绝了,打算寒假学得再深入些再广泛些,下学期开学后再去投递其他更好的公司。为了更加深入理解ArrayList,也为

我以为我学懂了数据结构,直到看了这个导图才发现,我错了

数据结构与算法思维导图

String s = new String(" a ") 到底产生几个对象?

老生常谈的一个梗,到2020了还在争论,你们一天天的,哎哎哎,我不是针对你一个,我是说在座的各位都是人才! 上图红色的这3个箭头,对于通过new产生一个字符串(”宜春”)时,会先去常量池中查找是否已经有了”宜春”对象,如果没有则在常量池中创建一个此字符串对象,然后堆中再创建一个常量池中此”宜春”对象的拷贝对象。 也就是说准确答案是产生了一个或两个对象,如果常量池中原来没有 ”宜春” ,就是两个。...

技术大佬:我去,你写的 switch 语句也太老土了吧

昨天早上通过远程的方式 review 了两名新来同事的代码,大部分代码都写得很漂亮,严谨的同时注释也很到位,这令我非常满意。但当我看到他们当中有一个人写的 switch 语句时,还是忍不住破口大骂:“我擦,小王,你丫写的 switch 语句也太老土了吧!” 来看看小王写的代码吧,看完不要骂我装逼啊。 private static String createPlayer(PlayerTypes p...

和黑客斗争的 6 天!

互联网公司工作,很难避免不和黑客们打交道,我呆过的两家互联网公司,几乎每月每天每分钟都有黑客在公司网站上扫描。有的是寻找 Sql 注入的缺口,有的是寻找线上服务器可能存在的漏洞,大部分都...

讲一个程序员如何副业月赚三万的真实故事

loonggg读完需要3分钟速读仅需 1 分钟大家好,我是你们的校长。我之前讲过,这年头,只要肯动脑,肯行动,程序员凭借自己的技术,赚钱的方式还是有很多种的。仅仅靠在公司出卖自己的劳动时...

上班一个月,后悔当初着急入职的选择了

最近有个老铁,告诉我说,上班一个月,后悔当初着急入职现在公司了。他之前在美图做手机研发,今年美图那边今年也有一波组织优化调整,他是其中一个,在协商离职后,当时捉急找工作上班,因为有房贷供着,不能没有收入来源。所以匆忙选了一家公司,实际上是一个大型外包公司,主要派遣给其他手机厂商做外包项目。**当时承诺待遇还不错,所以就立马入职去上班了。但是后面入职后,发现薪酬待遇这块并不是HR所说那样,那个HR自...

女程序员,为什么比男程序员少???

昨天看到一档综艺节目,讨论了两个话题:(1)中国学生的数学成绩,平均下来看,会比国外好?为什么?(2)男生的数学成绩,平均下来看,会比女生好?为什么?同时,我又联想到了一个技术圈经常讨...

副业收入是我做程序媛的3倍,工作外的B面人生是怎样的?

提到“程序员”,多数人脑海里首先想到的大约是:为人木讷、薪水超高、工作枯燥…… 然而,当离开工作岗位,撕去层层标签,脱下“程序员”这身外套,有的人生动又有趣,马上展现出了完全不同的A/B面人生! 不论是简单的爱好,还是正经的副业,他们都干得同样出色。偶尔,还能和程序员的特质结合,产生奇妙的“化学反应”。 @Charlotte:平日素颜示人,周末美妆博主 大家都以为程序媛也个个不修边幅,但我们也许...

MySQL数据库面试题(2020最新版)

文章目录数据库基础知识为什么要使用数据库什么是SQL?什么是MySQL?数据库三大范式是什么mysql有关权限的表都有哪几个MySQL的binlog有有几种录入格式?分别有什么区别?数据类型mysql有哪些数据类型引擎MySQL存储引擎MyISAM与InnoDB区别MyISAM索引与InnoDB索引的区别?InnoDB引擎的4大特性存储引擎选择索引什么是索引?索引有哪些优缺点?索引使用场景(重点)...

如果你是老板,你会不会踢了这样的员工?

有个好朋友ZS,是技术总监,昨天问我:“有一个老下属,跟了我很多年,做事勤勤恳恳,主动性也很好。但随着公司的发展,他的进步速度,跟不上团队的步伐了,有点...

我入职阿里后,才知道原来简历这么写

私下里,有不少读者问我:“二哥,如何才能写出一份专业的技术简历呢?我总感觉自己写的简历太烂了,所以投了无数份,都石沉大海了。”说实话,我自己好多年没有写过简历了,但我认识的一个同行,他在阿里,给我说了一些他当年写简历的方法论,我感觉太牛逼了,实在是忍不住,就分享了出来,希望能够帮助到你。 01、简历的本质 作为简历的撰写者,你必须要搞清楚一点,简历的本质是什么,它就是为了来销售你的价值主张的。往深...

玩转springboot启动banner定义所得

最近接手了一个springboot项目,不是不熟悉这个框架,启动时打印的信息吸引了我。 这不是我熟悉的常用springboot的打印信息啊,我打开自己的项目: 还真是的,不用默认的感觉也挺高大上的。一时兴起,就去研究了一下源代码,还正是有些收获,稍后我会总结一下。正常情况下做为一个老程序员,是不会对这种小儿科感兴趣的,不就是一个控制台打印嘛。哈哈! 于是出于最初的好奇,研究了项目的源代码。看到

带了6个月的徒弟当了面试官,而身为高级工程师的我天天修Bug......

即将毕业的应届毕业生一枚,现在只拿到了两家offer,但最近听到一些消息,其中一个offer,我这个组据说客户很少,很有可能整组被裁掉。 想问大家: 如果我刚入职这个组就被裁了怎么办呢? 大家都是什么时候知道自己要被裁了的? 面试软技能指导: BQ/Project/Resume 试听内容: 除了刷题,还有哪些技能是拿到offer不可或缺的要素 如何提升面试软实力:简历, 行为面试,沟通能...

离职半年了,老东家又发 offer,回不回?

有小伙伴问松哥这个问题,他在上海某公司,在离职了几个月后,前公司的领导联系到他,希望他能够返聘回去,他很纠结要不要回去? 俗话说好马不吃回头草,但是这个小伙伴既然感到纠结了,我觉得至少说明了两个问题:1.曾经的公司还不错;2.现在的日子也不是很如意。否则应该就不会纠结了。 老实说,松哥之前也有过类似的经历,今天就来和小伙伴们聊聊回头草到底吃不吃。 首先一个基本观点,就是离职了也没必要和老东家弄的苦...

男生更看重女生的身材脸蛋,还是思想?

往往,我们看不进去大段大段的逻辑。深刻的哲理,往往短而精悍,一阵见血。问:产品经理挺漂亮的,有点心动,但不知道合不合得来。男生更看重女生的身材脸蛋,还是...

为什么程序员做外包会被瞧不起?

二哥,有个事想询问下您的意见,您觉得应届生值得去外包吗?公司虽然挺大的,中xx,但待遇感觉挺低,马上要报到,挺纠结的。

当HR压你价,说你只值7K,你该怎么回答?

当HR压你价,说你只值7K时,你可以流畅地回答,记住,是流畅,不能犹豫。 礼貌地说:“7K是吗?了解了。嗯~其实我对贵司的面试官印象很好。只不过,现在我的手头上已经有一份11K的offer。来面试,主要也是自己对贵司挺有兴趣的,所以过来看看……”(未完) 这段话主要是陪HR互诈的同时,从公司兴趣,公司职员印象上,都给予对方正面的肯定,既能提升HR的好感度,又能让谈判气氛融洽,为后面的发挥留足空间。...

面试:第十六章:Java中级开发(16k)

HashMap底层实现原理,红黑树,B+树,B树的结构原理 Spring的AOP和IOC是什么?它们常见的使用场景有哪些?Spring事务,事务的属性,传播行为,数据库隔离级别 Spring和SpringMVC,MyBatis以及SpringBoot的注解分别有哪些?SpringMVC的工作原理,SpringBoot框架的优点,MyBatis框架的优点 SpringCould组件有哪些,他们...

面试阿里p7,被按在地上摩擦,鬼知道我经历了什么?

面试阿里p7被问到的问题(当时我只知道第一个):@Conditional是做什么的?@Conditional多个条件是什么逻辑关系?条件判断在什么时候执...

终于懂了TCP和UDP协议区别

终于懂了TCP和UDP协议区别

你打算用Java 8一辈子都不打算升级到Java 14,真香

我们程序员应该抱着尝鲜、猎奇的心态,否则就容易固步自封,技术停滞不前。

无代码时代来临,程序员如何保住饭碗?

编程语言层出不穷,从最初的机器语言到如今2500种以上的高级语言,程序员们大呼“学到头秃”。程序员一边面临编程语言不断推陈出新,一边面临由于许多代码已存在,程序员编写新应用程序时存在重复“搬砖”的现象。 无代码/低代码编程应运而生。无代码/低代码是一种创建应用的方法,它可以让开发者使用最少的编码知识来快速开发应用程序。开发者通过图形界面中,可视化建模来组装和配置应用程序。这样一来,开发者直...

面试了一个 31 岁程序员,让我有所触动,30岁以上的程序员该何去何从?

最近面试了一个31岁8年经验的程序猿,让我有点感慨,大龄程序猿该何去何从。

大三实习生,字节跳动面经分享,已拿Offer

说实话,自己的算法,我一个不会,太难了吧

程序员垃圾简历长什么样?

已经连续五年参加大厂校招、社招的技术面试工作,简历看的不下于万份 这篇文章会用实例告诉你,什么是差的程序员简历! 疫情快要结束了,各个公司也都开始春招了,作为即将红遍大江南北的新晋UP主,那当然要为小伙伴们做点事(手动狗头)。 就在公众号里公开征简历,义务帮大家看,并一一点评。《启舰:春招在即,义务帮大家看看简历吧》 一石激起千层浪,三天收到两百多封简历。 花光了两个星期的所有空闲时...

《经典算法案例》01-08:如何使用质数设计扫雷(Minesweeper)游戏

我们都玩过Windows操作系统中的经典游戏扫雷(Minesweeper),如果把质数当作一颗雷,那么,表格中红色的数字哪些是雷(质数)?您能找出多少个呢?文中用列表的方式罗列了10000以内的自然数、质数(素数),6的倍数等,方便大家观察质数的分布规律及特性,以便对算法求解有指导意义。另外,判断质数是初学算法,理解算法重要性的一个非常好的案例。

《Oracle Java SE编程自学与面试指南》最佳学习路线图(2020最新版)

正确选择比瞎努力更重要!

立即提问
相关内容推荐