应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 求助大侠帮忙，有关matlab模拟数据
北京石油化工学院2026年研究生招生接收调剂公告
301/3返回列表123下一页
查看: 2822  |  回复: 29
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

飞也如铃

新虫 (小有名气)

[求助] 求助大侠帮忙，有关matlab模拟数据

x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];

y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3,774];

上面是x,y的数据，公式是
请问哪位大侠帮忙给编辑一下公式代码啊，我自己用matlab模拟不出来，请哪位好心人帮帮忙！
回复此楼

» 收录本帖的淘帖专辑推荐

matlab

» 本帖已获得的红花（最新10朵）

xch09

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-11-22 20:34:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

飞也如铃(金币+2): 非常感谢您提供的这个信息，但是我找不到我需要的公式，请问我的公式适合fitting里面的哪个类型的公式呢，谢谢回复！ 2011-11-22 21:58:19
matlab里面有个cftool工具，拟合起来非常方便的
赞一下 回复此楼
2楼2011-11-22 21:15:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

拟合出来的结果
General model:
       f(x) = a+b*log(1-b*x)+(1-b)*log(1-b*x+x)
Coefficients (with 95% confidence bounds):
       a =      -1.309  (-9942, 9939)
       b =      -117.5  (-1.173e+006, 1.173e+006)

Goodness of fit:
  SSE: 0.4103
  R-square: 0.337
  Adjusted R-square: 0.2541
  RMSE: 0.2265
赞一下 回复此楼
3楼2011-11-22 21:25:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

用这个式子拟合貌似不好
赞一下 回复此楼
4楼2011-11-22 21:27:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

★ ★
飞也如铃(金币+2): 谢谢您的回复，可否把matlab里你输入的所有信息发到帖子了呢，我自己试了试，不知道是哪里出问题了，拟合不出来，谢谢！ 2011-11-22 22:15:30
dbb627(金币+2): 感谢应助 2011-11-22 22:57:51
重新拟合了一下
General model:
       f(x) = a+b*log(1-b*x)+(1-b)*log(1-b*x+x)
       where x is normalized by mean 0.6767 and std 0.03709
Coefficients (with 95% confidence bounds):
       a =       4.203  (4.157, 4.25)
       b =       0.294  (0.2726, 0.3154)

Goodness of fit:
  SSE: 0.03251
  R-square: 0.9475
  Adjusted R-square: 0.9409
  RMSE: 0.06375

这个好点
赞一下(1人) 回复此楼
5楼2011-11-22 21:32:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞也如铃

新虫 (小有名气)
引用回帖:
5楼: Originally posted by juloong at 2011-11-22 21:32:47:
重新拟合了一下
General model:
       f(x) = a+b*log(1-b*x)+(1-b)*log(1-b*x+x)
       where x is normalized by mean 0.6767 and std 0.03709
Coefficients (with 95% confidence bounds):
       a =  ...

x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];

y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];

plot(x,y,'*')

f = fittype('a+b*log(1-b*x)+(1-b)*log(1-b*x+x)');

%拟合函数
[c2,gof2] = fit(x,y,f)

% 显示拟合曲线与数据点的结果
plot(c2,'m')  % 绘制拟合曲线
hold on,plot(x,y,'*') %绘制原数据点
求大侠现身帮忙，我用上面的代码放到matlab里面，但是模拟出来的是一条直线，我的数据散点图是一曲线，可是我只会上面这个代码，请问大侠你用的代码是什么，能否发我一下呢？
赞一下 回复此楼
6楼2011-11-23 09:31:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

飞也如铃(金币+2): 谢谢您的回复，按照这个步骤拟合出来的效果不是很好，没有你第二次给的拟合效果好，请问你第二次的拟合是怎么弄的？ 2011-11-23 11:05:46
引用回帖:
6楼: Originally posted by 飞也如铃 at 2011-11-23 09:31:18:
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];

y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];

plot(x,y,'*')

f = ...

不需要这样写代码的，matlab有内建 的拟合工具cftool

过程如下
首先，在matlab中输入x,y的值
x = [0.768;0.711;0.685;0.671;0.659;0.657;0.655;0.654;0.653;0.654];
y = [4.584;4.375;4.234;4.128;4.041;3.972;3.913;3.861;3.815;3.774];

然后，键入命令cftool，这时会出现工具窗口，点data键，x data和y data分别选择x,y, 再点creat data set。然后关闭窗口，再回到工具窗口

然后fitting - new fit- type of fit - custmos equations
然后输入你的拟合公式，拟合就行了
赞一下 回复此楼
7楼2011-11-23 10:30:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

飞也如铃

新虫 (小有名气)
引用回帖:
7楼: Originally posted by juloong at 2011-11-23 10:30:35:
不需要这样写代码的，matlab有内建 的拟合工具cftool

过程如下
首先，在matlab中输入x,y的值
x = [0.768;0.711;0.685;0.671;0.659;0.657;0.655;0.654;0.653;0.654];
y = [4.584;4.375;4.234;4.128;4.041;3 ...

恩，非常感谢您的回复，我尝试了一下，给出的结果是
我自己拟合的效果不是很好，你拟合的这个效果比较好，请问你是调试过了吗，具体是怎么弄的呢？
这个是你拟合的数据
General model:
       f(x) = a+b*log(1-b*x)+(1-b)*log(1-b*x+x)
       where x is normalized by mean 0.6767 and std 0.03709
Coefficients (with 95% confidence bounds):
       a =       4.203  (4.157, 4.25)
       b =       0.294  (0.2726, 0.3154)

Goodness of fit:
  SSE: 0.03251
  R-square: 0.9475
  Adjusted R-square: 0.9409
  RMSE: 0.06375
赞一下 回复此楼
8楼2011-11-23 10:55:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juloong

木虫 (职业作家)

X射线衍/散射测试

【答案】应助回帖

多调一下参数选项，肯定能找到合适的值
赞一下 回复此楼
9楼2011-11-23 11:53:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dbb627

荣誉版主 (著名写手)

【答案】应助回帖

飞也如铃(金币+2): 谢谢您的回复，把你的代码放到matlab里能出来上面的结果，但是没有出来非线性拟合图，我想要这个非线性拟合图，请问怎么弄呢？ 2011-11-28 16:05:44
我给你写了一个
你的方程是lg应该是10为底对数
matlab里面log表示ln自然对数
CODE:
x = [0.768;
0.711;
0.685;
0.671;
0.659;
0.657;
0.655;
0.654;
0.653;
0.654];

y = [4.584;
4.375;
4.234;
4.128;
4.041;
3.972;
3.913;
3.861;
3.815;
3.774];
opts = fitoptions('Method','Nonlinear','Normalize','On');
opts.StartPoint =  [1 0.1];
ftype = fittype('a+b*log10(1-b*x)+(1-b)*log10(1-b*x+x)','options',opts);
[fresult,gof] = fit(x,y,ftype)

结果
fresult =

     General model:
     fresult(x) = a+b*log10(1-b*x)+(1-b)*log10(1-b*x+x)
       where x is normalized by mean 0.6767 and std 0.03709
     Coefficients (with 95% confidence bounds):
       a =       4.167  (4.126, 4.208)
       b =      0.1139  (0.05999, 0.1677)

gof =

           sse: 0.0207
       rsquare: 0.9665
           dfe: 8
    adjrsquare: 0.9623
          rmse: 0.0509
赞一下 回复此楼
The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.
10楼2011-11-23 23:48:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 飞也如铃 的主题更新
301/3返回列表123下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿9材料学硕297已过六级求调剂推荐 +11 adaie 2026-04-04 12/600 2026-04-05 19:04 by 蓝云思雨
[考研] 085500机械专硕初试288求调剂 +3 GZJguo666- 2026-04-05 3/150 2026-04-05 18:06 by jkddd
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 318求调剂 +11 ykyhsa 2026-04-05 13/650 2026-04-05 12:44 by aidndnjck
[考研] 302分求调剂 一志愿安徽大学085601 +5 zyx上岸！ 2026-04-04 5/250 2026-04-05 07:50 by 544594351
[考研] 26考研调剂0710 0860 +9 补补不补 2026-04-03 14/700 2026-04-04 23:32 by 果冻大王
[考研] 266分，一志愿电气工程，本科材料，求材料专业调剂 +11 哇呼哼呼哼 2026-04-01 12/600 2026-04-04 23:17 by 永字号
[考研] 材料383求调剂 +5 郭阳阳阳成 2026-04-04 5/250 2026-04-04 19:06 by dongzh2009
[考研] 求调剂：085600材料与化工，考材科基，总分319 +21 678lucky 2026-03-31 26/1300 2026-04-04 16:22 by dongzh2009
[考研] 一志愿东北大学085901土木专硕345求调剂 +3 zxt11111 2026-04-04 3/150 2026-04-04 14:21 by 土木硕士招生
[考研] 一志愿中国石油大学化学工程323分求调剂 +4 化工专硕323分 2026-04-03 6/300 2026-04-03 22:12 by dongzh2009
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 292求调剂 +21 是妍子也是研子 2026-03-30 22/1100 2026-04-03 21:44 by qlm5820
[考研] 考研调剂 +8 不爱喝饮料 2026-04-03 8/400 2026-04-03 16:40 by Mistake-J
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 0856材料与化工调剂，339 +14 10213207 2026-03-31 14/700 2026-04-02 21:01 by 1104338198
[考研] 260求调剂 +6 朱芷琳 2026-04-02 6/300 2026-04-02 20:27 by 6781022
[考研] 一志愿北京科技大学材料学硕328分求调剂 +6 1段时间 2026-03-31 7/350 2026-04-02 13:57 by 3041
[考研] 求调剂，一志愿南京师范大学计算机专硕，初试373，六级通过， +3 计算机追梦人 2026-04-01 3/150 2026-04-02 07:57 by fxue1114
[考研] 298求调剂 +4 什么是胖头鱼 2026-03-30 6/300 2026-04-01 22:06 by 客尔美德
信息提示
关闭
请填处理意见
关闭