当前位置: 首页 > news >正文

matlab不能盗版吗,matlab为了防止盗版,会不会篡改程序运行结果?这是明证

小弟不才,做了一个简单程序,如下:

x=[3.5  3.6  4.0  2.8  2.7  2.8  2.9  2.8  3.1  2.8 ....

3.0  4.3  3.1  2.3 3.0  2.7  2.5  3.8  2.7  4.1];

d=0.44;m=10;xx=sort(x,'descend');  %将数据由大到小排序

n=length(xx);

kk=[1:n];

p=100*(kk-d)./(n+1-2*d); %计算经验频率

e=mean(xx);  %均值

f=sqrt(sum((x./e-1).^2)/(n-1));   %Cv

X1=[1 2 4 10 20 50];  % 作业要求的重现期频率

X2=[0.01:0.01:99.99];  %  画图的横坐标间隔

for j=2:0.5:m

g=j*f;  %Cs

if  g>=2

disp('Cs不能大于2');

break;

end

h=4/(g^2);  %afa

a=e*(1-2*f/g); %a0

b=(e*f*g)/2; %bta

XX1=norminv(X1./100,0,1);

XX1=XX1-XX1(1);

pp=gaminv(X1./100,h,b)+a;  % 设计波高

XX2=norminv(p./100,0,1);

XX2=XX2-XX2(1);

pp2=gaminv(p./100,h,b)+a;   %原数据对应的累计率波高

figure;

X=[0.01 0.05 0.1 0.5 1 5 10 15 20 30 40 50 60 70 80 85 90 95 99 99.9 99.95 99.99]; %横坐标的值

Y=norminv(X/100,0,1);

Y=Y-Y(1);

for i=1:size(Y,2)

line([Y(i),Y(i)],[0 8]);

end

for i=0:0.5:8 % 10为纵坐标的格数,可以自已设置。

line([0 Y(end)],[i i]);

end

h=findobj('type','axes');set(h,'xtick',[],'ytick',[2:.5:8],'xlim',[0 Y(end)],'ylim',[0 ,8]);

XX=[0.01 1 2 10  40 50  70 90 ];%标横坐标的值,可以自己设置

YY=norminv(XX/100,0,1);YY=YY-YY(1);

for i=1:size(XX,2)

text('string',num2str(XX(i)),'HorizontalAlignment','center','pos',[YY(i),-0.5]);

end

xlabel('Probability');

ylabel('Date');

title('P-III Probability Plot');

hold on;

plot(XX2,pp2,'r*');

legend('data');

hold on;

XX3=norminv(X2./100,0,1);

XX3=XX3-XX3(1);

pp1=gaminv(X2./100,h,b)+a;

plot(XX3,pp1,':c');

legend('理论曲线');

hold on;

plot(XX1,pp,'bo');

legend('设计波高');

k=(j-2)/0.5+1;

RMSE(k)=sqrt(sum(pp2-xx).^2/n;

end

运行显示结果,函数已经超过了坐标范围。但实际上结果就在坐标范围内。原因错在h这个变量的计算上,一次循环只有一个h,但运行的时候有三个或者五个(已经在几个电脑运行过,不同电脑,结果不同),而且matlab给出的值与理论值差的太远。请问这算不算matlab的保护正版措施,还是我的程序有错?

相关文章:

  • matlab 复权数据,〖Matlab〗基于通达信股价数据的复权处理(fantuanxiaot版本)
  • 微信小程序向php传递数据,微信小程序 跳转传递数据的方法
  • mysql8默认存储,MySQL 8.0安装
  • php berkeleydb,berkeleydb-5.1
  • matlab vs 打包exe文件路径,Matlab中調用VS編譯的exe文件並傳遞變量 的方法
  • php 访问受保护的属性,php – 我们应该直接访问受保护的属性还是使用getter?
  • PHP限制请求类型,Laravel :API 请求频率限制(Throttle中间件),自定义返回JSON类型,自定义时间...
  • oracle归档增长快,归档日志增长过快的问题
  • oracle 行转列 顺序,Oracle PIVOT 行转列的单行小计如何实现?
  • oracle asm安装详情,oracle ASM安装记录
  • oracle in _all和any的区别,英语all和any的用法区别
  • linux用户隔离,CentOS下实现用户隔离模式的SFTP
  • linux kdump保存的目录,学习猿地-Linux中的Kdump服务
  • linux文件拆分,Linux文件分割和合并
  • linux 脚本遍历文件夹,Linux四种方法遍历执行目录下的全部脚本
  • hexo+github搭建个人博客
  • 【腾讯Bugly干货分享】从0到1打造直播 App
  • CSS实用技巧
  • Essential Studio for ASP.NET Web Forms 2017 v2,新增自定义树形网格工具栏
  • input的行数自动增减
  • java取消线程实例
  • JS基础之数据类型、对象、原型、原型链、继承
  • Python连接Oracle
  • React组件设计模式(一)
  • SpiderData 2019年2月16日 DApp数据排行榜
  • Vue.js 移动端适配之 vw 解决方案
  • vue从入门到进阶:计算属性computed与侦听器watch(三)
  • yii2中session跨域名的问题
  • 基于Vue2全家桶的移动端AppDEMO实现
  • 计算机在识别图像时“看到”了什么?
  • 前端知识点整理(待续)
  • 如何利用MongoDB打造TOP榜小程序
  • 深入 Nginx 之配置篇
  • 什么软件可以剪辑音乐?
  • 小程序开发之路(一)
  • (6)添加vue-cookie
  • (BFS)hdoj2377-Bus Pass
  • (DFS + 剪枝)【洛谷P1731】 [NOI1999] 生日蛋糕
  • (附源码)ssm旅游企业财务管理系统 毕业设计 102100
  • (附源码)计算机毕业设计SSM疫情社区管理系统
  • (四)【Jmeter】 JMeter的界面布局与组件概述
  • (转载)hibernate缓存
  • ./configure,make,make install的作用(转)
  • .NET 6 Mysql Canal (CDC 增量同步,捕获变更数据) 案例版
  • .Net环境下的缓存技术介绍
  • .NET连接数据库方式
  • .net项目IIS、VS 附加进程调试
  • .net中的Queue和Stack
  • @RequestBody与@ResponseBody的使用
  • @取消转义
  • [ Algorithm ] N次方算法 N Square 动态规划解决
  • [2009][note]构成理想导体超材料的有源THz欺骗表面等离子激元开关——
  • [30期] 我的学习方法
  • [383] 赎金信 js
  • [Android Pro] AndroidX重构和映射