网页
资讯
视频
图片
知道
文库
贴吧
地图
采购
进入贴吧
全吧搜索
吧内搜索
搜贴
搜人
进吧
搜标签
日
一
二
三
四
五
六
签到排名:今日本吧第
个签到,
本吧因你更精彩,明天继续来努力!
本吧签到人数:0
一键签到
成为超级会员,使用一键签到
一键签到
本月漏签
0
次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行
补签
。
连续签到:
天 累计签到:
天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
02月09日
漏签
0
天
matlab吧
关注:
290,729
贴子:
1,655,822
看贴
图片
吧主推荐
游戏
7
回复贴,共
1
页
<<返回matlab吧
>0< 加载中...
求助
fmincon使用错误
只看楼主
收藏
回复
荒誕演绎
自带板凳
3
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
[lbk]y_optimized, fval, exitflag[rbk] = fmincon(objective_function, initial_guess, [lbk][rbk], [lbk][rbk], A_eq, b_eq, lb, ub, [lbk][rbk], options);
我有点不知道哪里有问题,这个影不影响运行结果啊大佬们帮我看看
lhmhz
内牛满面
13
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
这类错误可能出现在自定义函数中。你可以把代码完整的贴出来,才能发现你的问题。
荒誕演绎
自带板凳
3
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
% 定义原始曲线函数
x_func = @(t) t;
y_func = @(t) t.^3 + t;
% 定义参数范围和采样点数量
num_points = 10; % 采样点数量
t_values = linspace(0, 1, num_points); % 参数范围
% 计算采样点的坐标
x_samples = x_func(t_values);
y_samples = y_func(t_values);
% 计算采样点的切线向量和法向量
dx_dt = gradient(x_samples) ./ gradient(t_values);
dy_dt = gradient(y_samples) ./ gradient(t_values);
tangent_vectors = [dx_dt; dy_dt];
% 计算采样点的曲率
dx2_dt2 = gradient(dx_dt) ./ gradient(t_values);
dy2_dt2 = gradient(dy_dt) ./ gradient(t_values);
curvatures = abs(dx_dt .* dy2_dt2 - dy_dt .* dx2_dt2) ./ (dx_dt.^2 + dy_dt.^2).^(3/2);
x_values = x_samples;
%y_initial = []; % 最初的 y 值
k = curvatures; % 对应点的曲率
% 生成 10个随机数(继续优化)
random_numbers = -3 + 6 * rand(1, 10);
y_initial = random_numbers;
% 定义拟合的多项式次数
degree = 3; % 这里选择 3 次多项式
% 使用 polyfit 函数进行多项式拟合
coeffs = polyfit(x_values, y_initial, degree);
% 使用拟合参数来预测拟合曲线上的点
x_fit = linspace(min(x_values), max(x_values), 100); % 生成更密集的横坐标
y_fit_initial = polyval(coeffs, x_fit);
% 绘制原始数据和拟合曲线
figure;
plot(x_values, y_initial, 'o', x_fit, y_fit_initial, '-', 'LineWidth', 1.5);
hold on;
% 绘制原始曲线
x_original = linspace(0, 1, 100);
y_original = x_original.^3 + x_original;
plot(x_original, y_original, 'r--', 'LineWidth', 1.5);
legend('原始数据', '最初的拟合曲线', '原始曲线');
title('最初的拟合曲线与原始曲线');
% 假设已经拟合了数据并得到了多项式系数 coeffs
% 定义多项式函数
poly_func = @(x) polyval(coeffs, x);
% 定义二阶导数函数
poly_second_derivative = @(x) 2 * coeffs(1) + 6 * coeffs(2) * x + 12 * coeffs(3) * x.^2;
% 计算最初的曲率
curvature_values_initial = abs(poly_second_derivative(x_values)) ./ (1 + polyval(coeffs, x_values).^2).^(3/2);
% 定义目标函数,使得每个曲率的差值不大于1
%objective_function = @(y) sum((max(abs(abs(poly_second_derivative(x_values)) ./ (1 + polyval(polyfit(x, y, degree), x_values).^2).^(3/2) - k) - 1, 0)).^2);
% 定义目标函数,使得每个曲率的差值都不大于给定曲率的一定百分比
percentage_threshold = 0.01; % 设定百分比阈值为1%
objective_function = @(y) sum((max(abs(abs(poly_second_derivative(x_values)) ./ (1 + polyval(polyfit(x_values, y, degree), x_values).^2).^(3/2) - k) - percentage_threshold * k, 0)).^2);
% 定义目标函数
lhmhz
内牛满面
13
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
运行没有出现错误啊,
荣光金陵
内牛满面
13
该楼层疑似违规已被系统折叠
隐藏此楼
查看此楼
你是这个优化工具箱没有安装好吧,就用不了
登录百度账号
扫二维码下载贴吧客户端
下载贴吧APP
看高清直播、视频!
贴吧页面意见反馈
违规贴吧举报反馈通道
贴吧违规信息处理公示